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Abstract 

This article reports world averages for measurements on 6-hadron properties obtained 
by the Heavy Flavor Averaging Group (HFAG) using the available results as of summer 
2004 conferences. In the averaging, the input parameters used in the various analyses are 
adjusted (rescaled) to common values, and all known correlations are taken into account. 
The averages include fe-hadron lifetimes, S-oscillation (mixing) parameters, semileptonic 
decay parameters, rare decay branching fractions, and CP violation measurements. 
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1 Introduction 



The flavor dynamics is one of tlie important elements in understanding tlie nature of particle 
physics. The accurate knowledge of properties of heavy flavor hadrons, especially b hadrons, 
play an essential role for determination of the Cabbibo-Kobayashi-Maskawa (CKM) matrix PJ. 
Since asymmetric B factories started their operation, available amounts of B meson samples 
has been dramatically increased and the accuracies of measurements have been improved. 

The Heavy Flavor Averaging Group (HFAG) has been formed, continuing the activities of 
LEP Heavy Flavor Steering group j2], to provide the averages for measurements dedicated to 
the 6-flavor related quantities. The HFAG consists of representatives and contacts from the 
experimental groups: BABAR, Belle, CDF, CLEO, D0, and LEP. 

The HFAG is currently organized into four subgroups. 

• the "Lifetime and mixing" group provides averages for 6-hadron lifetimes, 6-hadron frac- 
tions in T[4:S) decay and high energy collisions, and various parameters in B^ and B^ 
oscillation (mixing). 

• the "Semileptonic B decays" group provides averages for inclusive and exclusive -B-decay 
branching fractions, and best values of the CKM matrix elements \Vcb\ and \Vub\- 

• the ^^CP{t) and Unitarity Triangle angles" group provides averages for time-dependent 
CP asymmetry parameters and angles of the unitarity triangles. 

• the "Rare decays" group provides averages of branching fractions and their asymmetries 
between B and B for charmless mesonic, radiative, leptonic, and baryonic B decays. 

The first two subgroups continue the activities from LEP working groups with some reor- 
ganization (merging four groups into two groups). The latter two groups are newly formed to 
take care of new results which are available from asymmetric B factory experiments. 

In this article, we report the world averages using the available results as of summer 2004 
conferences (ICHEP04 and FPCP04). All results that are publicly available, including recent 
preliminary results, are used in averages. We do not use preliminary results which remain 
unpublished for a long time or for which no publication is planned. Close contacts have been 
established between representatives from the experiments and members of different subgroups in 
charge of the averages, to ensure that the data are prepared in a form suitable for combinations. 

We do not scale the error of an average (as is presently done by the Particle Data Group [3]) 
in case x^/dof > 1, where dof is the number of degrees of freedom in the average calculation. 
In this case, we examine the systematics of each measurement and try to understand them. 
Unless we find possible systematic discrepancies between the measurements, we do not make 
any special treatment for the calculated error. We provide the confidence level of the fit so 
that one can know the consistency of the measurements included in the average. We attach a 
warning message in case that some special treatment is done or the approximation used in the 
average calculation may not be good enough {e.g., Gaussian error is used in averaging though 
the likelihood indicates non-Gaussian behavior). 

Section 121 describes the methodology for averaging various quantities in the HFAG. In the 
averaging, the input parameters used in the various analyses are adjusted (rescaled) to common 
values, and, where possible, known correlations are taken into account. The general philosophy 
and tools for calculations of averages are presented. 
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Sections OHHl describe the averaging of the quantities from each subgroup mentioned above. 
A summary of the averages described in this article is given in Sec. [7| 
The complete listing of averages and plots described in this article are also available on the 
HFAG Web page: 

http : //www. slac . Stanford. edu/xorg/hfag and 
http://belle.kek.jp/mirror/hfag (KEK mirror site). 

2 Methodology 

The general averaging problem that HFAG faces is to combine the information provided by 
different measurements of the same parameter, to obtain our best estimate of the parameter's 
value and uncertainty. The methodology described here focuses on the problems of combining 
measurements performed with different systematic assumptions and with potentially-correlated 
systematic uncertainties. Our methodology relies on the close involvement of the people per- 
forming the measurements in the averaging process. 

Consider two hypothetical measurements of a parameter x, which might be summarized as 

X = xi± 6xi ± Aa^i^i ± Ax2^i ■ ■ ■ 

X = X2 ± 6X2 ± Axi,2 ± AX2,2 • • • , 

where the 6xk are statistical uncertainties, and the Axi^k are contributions to the systematic 
uncertainty. One popular approach is to combine statistical and systematic uncertainties in 
quadrature 

x = xi± {6x1 © Axi,i © Ax2,i © . . .) 

X = X2± {6X2 © Axi,2 © AX2,2 © • • •) 

and then perform a weighted average of Xi and X2, using their combined uncertainties, as if 
they were independent. This approach suffers from two potential problems that we attempt 
to address. First, the values of the x^ may have been obtained using different systematic 
assumptions. For example, different values of the B'^ lifetime may have been assumed in 
separate measurements of the oscillation frequency Arrid. The second potential problem is 
that some contributions of the systematic uncertainty may be correlated between experiments. 
For example, separate measurements of Am^ may both depend on an assumed Monte-Carlo 
branching fraction used to model a common background. 

The problems mentioned above are related since, ideally, any quantity that x^ depends 
on has a corresponding contribution Axi^k to the systematic error which reflects the uncertainty 
Ayi on Hi itself. We assume that this is the case, and use the values of yi and Ayi assumed 
by each measurement explicitly in our averaging (we refer to these values as yi^k and Ayi^k 
below). Furthermore, since we do not lump all the systematics together, we require that each 
measurement used in an average have a consistent definition of the various contributions to the 
systematic uncertainty. Different analyses often use different decompositions of their systematic 
uncertainties, so achieving consistent definitions for any potentially correlated contributions 
requires close coordination between HFAG and the experiments. In some cases, a group of 
systematic uncertainties must be lumped to obtain a coarser description that is consistent 
between measurements. Systematic uncertainties that are uncorrelated with any other sources 



5 




Figure 1: The left-hand plot, (a), compares the 68% confidence-level contours of a hypothetical 
measurement's unconstrained (large ellipse) and constrained (filled ellipse) likelihoods, using 
the Gaussian constraint on yi represented by the horizontal band. The solid error bars represent 
the statistical uncertainties, (j{x) and cr(?/j), of the unconstrained likelihood. The dashed error 
bar shows the statistical error on x from a constrained simultaneous fit to x and yi. The 
right-hand plot, (b), illustrates the method described in the text of performing fits to x only 
with yi fixed at different values. The dashed diagonal line between these fit results has the 
slope p{x,yi)a{yi)/a{x) in the limit of a parabohc unconstrained likelihood. The result of the 
constrained simultaneous fit from (a) is shown as a dashed error bar on x. 



of uncertainty appearing in an average are lumped with the statistical error, so that the only 
systematic uncertainties treated explicitly are those that are correlated with at least one other 
measurement via a consistently-defined external parameter y^. When asymmetric statistical 
or systematic uncertainties are quoted, we symmetrize them since our combination method 
implicitly assumes parabolic likelihoods for each measurement. 

The fact that a measurement of x is sensitive to the value of yi indicates that, in principle, 
the data used to measure x could equally-well be used for a simultaneous measurement of x and 
yi, as illustrated by the large contour in Fig. Wi^) for a hypothetical measurement. However, 
we often have an external constraint Ayi on the value of yi (represented by the horizontal band 
in Fig. ^a)) that is more precise than the constraint cr{yi) from our data alone. Ideally, in 
such cases we would perform a simultaneous fit to x and yi, including the external constraint, 
obtaining the filled {x, y) contour and corresponding dashed one-dimensional estimate of x 
shown in Fig. ^a). Throughout, we assume that the external constraint Ayi on yi is Gaussian. 

In practice, the added technical complexity of a constrained fit with extra free parameters 
is not justified by the small increase in sensitivity, as long as the external constraints Ayi are 
sufficiently precise when compared with the sensitivities cr(yi) to each yi of the data alone. 
Instead, the usual procedure adopted by the experiments is to perform a baseline fit with all yi 
fixed to nominal values yi^Q, obtaining x = Xoi:6x. This baseline fit neglects the uncertainty due 
to Ayi, but this error can be mostly recovered by repeating the fit separately for each external 
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parameter Ui with its value fixed at yi = Q + A?/j to obtain x = ± 5x, as illustrated in 
Fig. ^b). The absolute shift, |xj^o ~ xo\, in the central value of x is what the experiments 
usually quote as their systematic uncertainty Axi on x due to the unknown value of i/i. Our 
procedure requires that we know not only the magnitude of this shift but also its sign. In the 
limit that the unconstrained data is represented by a parabolic likelihood, the signed shift is 
given by 

Axi = p{x, yi)^Y~\ ^yi^ (1) 

where (t{x) and p{x,yi) are the statistical uncertainty on x and the correlation between x and 
yi in the unconstrained data. While our procedure is not equivalent to the constrained fit with 
extra parameters, it yields (in the limit of a parabolic unconstrained likelihood) a central value 
Xq that agrees to 0{Ayi/a{yi)y and an uncertainty 6x © Axi that agrees to 0{Ayi/a{yi)y. 

In order to combine two or more measurements that share systematics due to the same 
external parameters yi, we would ideally perform a constrained simultaneous fit of all data 
samples to obtain values of x and each y^, being careful to only apply the constraint on each yi 
once. This is not practical since we generally do not have sufficient information to reconstruct 
the unconstrained likelihoods corresponding to each measurement. Instead, we perform the 
two-step approximate procedure described below. 

Figs.|2^a,b) illustrate two statistically-independent measurements, xi±(5xi©Axi,i) and X2± 
{SxiQ)Axi^2), of the same hypothetical quantity x (for simplicity, we only show the contribution 
of a single correlated systematic due to an external parameter yi). As our knowledge of the 
external parameters yi evolves, it is natural that the different measurements of x will assume 
different nominal values and ranges for each yi. The first step of our procedure is to adjust the 
values of each measurement to refiect the current best knowledge of the values and ranges 
Ay'i of the external parameters yi, as illustrated in Figs. I2tc,b). We adjust the central values 
Xk and correlated systematic uncertainties Axi^k linearly for each measurement (indexed by A;) 
and each external parameter (indexed by i): 

x'k =Xk + Y^ {y[ - yi,k) (2) 



^yi,k 
^yi,k 



Ax[^, = Ax,,, ■ ^ . (3) 



This procedure is exact in the limit that the unconstrained likelihoods of each measurement is 
parabolic. 

The second step of our procedure is to combine the adjusted measurements, x'^ ± {6xk © 
Ax'). I © Ax'f. 2 © . . .) using the chi-square 



xlomh{x,yi,y2, • • •) = XI ^ 



x'k - I X 



El i\ ^"^i k 



2 



and then minimize this to obtain the best values of x and yi and their uncertainties, as 
illustrated in Fig. El Although this method determines new values for the yi, we do not report 
them since the Axi,k reported by each experiment are generally not intended for this purpose 
(for example, they may represent a conservative upper limit rather than a true refiection of a 
68% confidence level). 
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Figure 2: The upper plots, (a) and (b), show examples of two individual measurements to be 
combined. The large ellipses represent their unconstrained likelihoods, and the filled ellipses 
represent their constrained likelihoods. Horizontal bands indicate the different assumptions 
about the value and uncertainty of Ui used by each measurement. The error bars show the 
results of the approximate method described in the text for obtaining x by performing fits 
with Hi fixed to different values. The lower plots, (c) and (d), illustrate the adjustments to 
accommodate updated and consistent knowledge of yi described in the text. Hollow circles 
mark the central values of the unadjusted fits to x with y fixed, which determine the dashed 
line used to obtain the adjusted values. 
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Figure 3: An illustration of the combination of two hypothetical measurements of x using the 
method described in the text. The ellipses represent the unconstrained likelihoods of each 
measurement and the horizontal band represents the latest knowledge about yi that is used 
to adjust the individual measurements. The filled small ellipse shows the result of the exact 
method using £comb and the hollow small ellipse and dot show the result of the approximate 
method using Xcomb- 
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For comparison, the exact method we would perform if we had the unconstrained hkehhoods 
jCkix,yi,y2, ■ ■ ■) available for each measurement is to minimize the simultaneous constrained 
likelihood 

>Ccomb(a;, 1/1, 1/2, • • •) = n • • •) n ' 

k i 

with an independent Gaussian external constraint on each rji 

(6) 

The results of this exact method are illustrated by the filled ellipses in Figs. El^a,b), and agree 
with our method in the limit that each Ck is parabolic and that each Ay,- ^ cr(yi). In the case 
of a non-parabolic unconstrained likelihood, experiments would have to provide a description 
of £k itself to allow an improved combination. In the case of some a{yi) ~ Ay'^, experiments 
are advised to perform a simultaneous measurement of both x and y so that their data will 
improve the world knowledge about y. 

The algorithm described above is used as a default in the averages reported in the following 
sections. For some cases, somewhat simplified or more complex algorithms are used and noted 
in the corresponding sections. 

Following the prescription described above, the central values and errors are rescaled to a 
common set of input parameters in the averaging procedures, according to the dependency on 
any of these input parameters. We try to use the most up-to-date values for these common 
inputs and the same values among the HFAG subgroups. For the parameters whose averages 
are produced by the HFAG, we use the updated values in the current update cycle. For other 
external parameters, we use the most recent PDG values. 

The parameters and values used in this update cycle are listed in each subgroup section. 

3 Averages of 6-hadron fractions, lifetimes and mixing 
parameters 

Quantities such as 6-hadron production fractions, 6-hadron lifetimes, and neutral 5-meson os- 
cillation frequencies have been measured for many years at high-energy colliders, namely LEP 
and SLC (e+e~ colliders at y/s = mz) as well as Tevatron Run I {pp collider at ^/s = 1.8 TeV). 
More recently, precise measurements of the and lifetimes, as well as of the oscil- 
lation frequency, have also been performed at the asymmetric B factories, KEKB and PEPII 
{e~^e~ colliders at ^/s = mT{AS))- In most cases, these basic quantities, although interesting 
by themselves, can now be seen as necessary ingredients for the more complicated and refined 
analyses being currently performed at the asymmetric B factories and at the Tevatron Run 
II (a/s = 2 TeV), in particular the time-dependent CP measurements. It is therefore impor- 
tant that the best experimental values of these quantities continue to be kept up-to-date and 
improved. 

In several cases, the averages presented in this chapter are indeed needed and used as input 
for the results given in the subsequent chapters. However, within this chapter, some averages 
need the knowledge of other averages in a circular way. This "coupling" , which appears through 
the 6-hadron fractions whenever inclusive or semi-exclusive measurements have to be considered. 



exp 



Vi - Vi 
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has been reduced significantly in the last years with increasingly precise exclusive measurements 
becoming available. To cope with this circularity, a rather involved averaging procedure had 
been developed, in the framework of the former LEP Heavy Flavour Steering Group. This is 
still in use now (details can be found in j2]), although simplifications can be envisaged in the 
future when even more precise exclusive measurements become available. 

3.1 6-hadron production fractions 

We consider here the relative fractions of the different 6-hadron species found in an unbiased 
sample of weakly-decaying b hadrons produced under some specific conditions. The knowledge 
of these fractions is useful to characterize the signal composition in inclusive 6-hadron analyses, 
or to predict the background composition in exclusive analyses. Many analyses in B physics 
need these fractions as input. We distinguish here the following two conditions: T{4S) decays 
and high-energy collisions. 

3.1.1 b-hadron fractions in Y(4S) decays 

Only pairs of the two lightest (charged and neutral) B mesons can be produced in T{4S) decays, 
and it is enough to determine the following branching fractions: 

= TiT{AS)^B+B-)/TUn^S)), (7) 

^00 ^ r(r(45)^i?°5')/rtot(r(4^)). (8) 

In practice, most analyses measure their ratio 

^+-/oo _ /+-//00 _ Y{T{AS) B+B-)/T{T{AS) ^ B"^') , (9) 

which is easier to access experimentally. Since an inclusive (but separate) reconstruction of 
B~^ and B'^ is difficult, specific exclusive decay modes, B~^ — * and S° — ^ x^, are usually 
considered to perform a measurement of i?"*""/*"^, whenever they can be related by isospin 
symmetry (for example B~^ —>■ J /ipK^ and B^ — > J/tpK^). Under the assumption that r(i?"^ — > 
x~^) = r{B^ x^), i.e. that isospin invariance holds in these B decays, the ratio of the number 
of reconstructed B^ x~^ and B^ x^ mesons is proportional to 

/+- B{B+ ^ x+) _ /+- r{B+ ^ x+) t{B+) _ /+- t{B+) 
P^B{B^-^x^) poT{BO xO)t{BO) /oo r(50) ' ^ ' 

where t{B^) and t{B^) are the B^ and B^ lifetimes respectively. Hence the primary quantity 
measured in these analyses is W /'^°r(i?+)/r(-B°), and the extraction of with this 

method therefore requires the knowledge of the t{B~^)/t{B^) lifetime ratio. 

The published measurements of R~^~^^^ are listed in Tabled together with the corresponding 
assumed values of t{B^) /t{B^). All measurements are based on the above-mentioned method, 
except the one from Belle, which is a by-product of the B^ mixing frequency analysis using 
dilepton events (but note that it also assumes isospin invariance, namely r{B~^ — *■ ^"'"X) = 
r[B^ The latter is therefore treated in a slightly different manner in the following 

procedure used to combine these measurements: 
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Table 1: Published measurements of the B^/B^ production ratio in T[4S) decays, together 
with their average (see text). Systematic uncertainties due to the imperfect knowledge of 



T 



B+)/t{B°) are included. 



Experiment Ref. Decay modes Published value of Assumed value 

and year or method = /+"//°° of r(E+)/r(E°) 

CLEO, 2001 4j JJ^JK^ 1.04 ± 0.07 ± 0.04 1.066 ± 0.024 

BABAR, 2002 y (cc)irW 1.10 ± 0.06 ± 0.05 1.062 ± 0.029 

CLEO, 2002 ^ D*eiy 1.058 ± 0.084 ± 0.136 1.074 ±0.028 

Belle, 2003 dilepton events 1.01 ± 0.03 ± 0.09 1.083 ± 0.017 

BABAR, 2004 ^ J/^K 1.006 ± 0.036 ± 0.031 1.083 ± 0.017 

Average 1.010 ±0.038 (tot) 1.081 ± 0.015 



• each published value of W''/^'^ from CLEO and BABAR is first converted back to the 
original measurement of '^{B^) I '^(B^) ? using the value of the lifetime ratio assumed 
in the corresponding analysis; 

• a simple weighted average of these original measurements of R~^~^^^ t{B^)/t{B^) from 
CLEO and BABAR (which do not depend on the assumed value of the lifetime ratio) is 
then computed, assuming no statistical or systematic correlations between them; 

• the weighted average of ^^'^ t{B~^)/t{B^) is converted into a value of using 
the latest average of the hfetime ratios, t{B+)/t{B^) = 1.081 ± 0.015 (see Sec. KT^ : 

• the Belle measurement of is adjusted to the current values of t{B^) = 1.534 ± 
0.013 ps and t{B+)/t{B^) = 1.081 ± 0.015 (see Sec. KT^ . using the quoted systematic 
uncertainties due to these parameters; 

• the combined value of from CLEO and BABAR is averaged with the adjusted value 
of 7?"*""/°° from Belle, assuming a 100% correlation of the systematic uncertainty due to 
the limited knowledge on t{B^)/t{B^)] no other correlation is considered. 

The resulting global average, 

f+- 

= :L_ = 1.010 ±0.038, (11) 

is consistent with an equal production of charged and neutral B mesons. 

On the other hand, the BABAR collaboration has recently performed a direct measurement 
of the /"^ fraction using a novel method, which does not rely on isospin symmetry nor requires 
the knowledge of t{B~^)/t{B^). Its preliminary analysis, based on a comparison between the 
number of events where a single B° D*~i~^u decay could be reconstructed and the number 
of events where two such decays could be reconstructed, yields jH] 

/°° = 0.486 ± 0.010 (stat) ± 0.009 (syst) . (12) 

The two results of Eqs. pi|) and (fT^ are of very different natures and completely indepen- 
dent of each other. Their product is equal to = 0.491 ± 0.023, while another combination 
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of them gives /"*"" + f^^ = 0.977 ± 0.033, compatible with unity. Assuming /+" + Z*"^ = 1, also 
consistent with CLEO's observation that the fraction of T(45') decays to BB pairs is larger 
than 0.96 at 95% CL [TU], the results of Eqs. (fTTj) and (fT^ can be averaged (first converting 
Eq. (fTT|) into a value of = + 1)) to yield the following more precise estimates: 

f+- 

/°° = 0.494 ± 0.008 , /+- = l-/°° = 0.506 ±0.008, = 1.026 ± 0.032 . (13) 

3.1.2 6-hadron fractions at high energy 

At high energy, all species of weakly-decaying b hadrons can be produced. We assume here 
that the fractions of these different species are the same in unbiased samples of high-px b jets 
originating from decays or from pp collisions at the Tevatron, either directly or in strong and 
electromagnetic decays of excited b hadrons. This hypothesis is plausible considering that, in 
both cases, the last step of the jet hadronization is a non-perturbative QCD process occurring 
at a scale of order /Iqcd- On the other hand, there is no strong argument to claim that 
these fractions should be strictly equal, so this assumption should be checked experimentally. 
Although the available data is not quite sufficient at this time to perform a significant check, 
it is expected that the new data from Tevatron Run II will soon improve this situation and 
allow to confirm or infirm this assumption with reasonable confidence. Meanwhile, the attitude 
adopted here is that these fractions are assumed to be equal at all high-energy colliders until 
demonstrated otherwise by experiment.^ 

Contrary to what happens in the charm sector where the fractions of and are 
different, the relative amount of B~^ and B^ is not affected by the electromagnetic decays of 
excited B~^* and B^* states and strong decays of excited 5+** and B^** states. Decays of the 
type B'^** — > B^*'^K also contribute to the B^ and B^ rates, but with the same magnitude if 
mass effects can be neglected. We therefore assume equal production of B^ and B^ . We also 
neglect the production of weakly-decaying states made of several heavy quarks (like B'^ and 
other heavy baryons) which is known to be very small. Hence, for the purpose of determining 
the 6-hadron fractions, we use the constraints 

fu = fd and fu + fd + fs + /baryon = 1 , (14) 

where /„, fd, fs and /baryon are the unbiased fractions of B^ , B^, B^ and 6-baryons, respectively. 

The LEP experiments have measured /, x i3(5° D-i+iyeX) B{b ^ A^) x 13{AI 
A+i-UiX) [Illig and B{b S^) x B{S^ S-i-V^X) from partially reconstructed 

final states including a lepton, /baryon from protons identified in b events [1^], and the production 
rate of charged b hadrons jT7j. The various 6-hadron fractions have also been measured at 
CDF using electron- charm final states ^Hl and double semileptonic decays with (pi and K*i 
final states fHl- AH these published results have been combined following the procedure and 
assumptions described in |2j to yield fu = fd = 0.403 ± 0.011, fs = 0.088 ± 0.021 and /baryon = 
0.107± 0.019 under the constraints of Eq. (fT^ . For this combination, other external inputs are 
used, e.g. the branching ratios of B mesons to final states with a D, D* or D** in semileptonic 
decays, which are needed to evaluate the fraction of semileptonic B'^ decays with a D~ in the 
final state. 

"'^It is not unlikely that the 6-liadron fractions in low-pT jets at a hadronic machine be different; in particular, 
beam-remnant effects may enhance the 6-baryon production. 
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Table 2: Fractions of the different 6-hadron species in an unbiased sample of weakly-decaying 
b hadrons produced at high energy, obtained from both direct and mixing measurements. 



6-hadron 




Fraction 


Correlation 


coefficients 


species 






with fd = /, 


, and fs 




fd fu 


= 0.398 ±0.010 








fs 


= 0.104 ±0.015 


-0.566 




b baryons 


/baryon 


= 0.100 ±0.017 


-0.712 


-0.176 



Time-integrated mixing analyses performed with lepton pairs from bb events produced at 
high-energy colliders measure the quantity 

X = fdXci + f:Xs, (15) 

where and are the fractions of and B^ hadrons in a sample of semileptonic 6-hadron 
decays, and where Xd and Xs are the B'^ and B^ time-integrated mixing probabilities. Assuming 
that all b hadrons have the same semileptonic decay width implies // = /j-Rj, where Ri = Tj/rf, 
is the ratio of the lifetime Tj of species i to the average 6-hadron lifetime = /jTj. Hence 
measurements of the mixing probabilities x, Xd and Xs can be used to improve our knowledge on 
the fu, fd, fs and /baryon fractions. In practice, the above relations yield another determination 
of fs obtained from /baryon and mixing information, 

r -'^ (1 ± (1 /baryon -Rbaryon)Xd (TR^ 

Rs (1 ± r)Xs - Xd 

where r = Ru/Rd = t{B+)/t{B^). 

The published measurements of x performed by the LEP experiments have been combined 
by the LEP Electroweak Working Group to yield x = 0.1257 ± 0.0042 This can be 

compared with a recent measurement from CDF, x = 0.152 ± 0.013 [21], obtained from an 
analysis of the Run I data. The two estimates deviate from each other by 1.9 a, and could be 
an indication that the fractions of b hadrons produced at the Z peak or at the Tevatron are not 
the same. Although this discrepancy is not very significant it should be carefully monitored 
in the future. We choose to combine these two results in a simple weighted average, assuming 
no correlations, and, following the PDG prescription, we multiply the combined uncertainty by 
1.9 to account for the discrepancy. Our world average is then 

X = 0.1282 ±0.0077. (17) 

Introducing the latter result in Eq. (fTK|) . together with our world average Xd = 0.186±0.003 
(see Eq. fjS^ of Sec. I3.3.1|) . the assumption Xs = 1/2 (justified by the large value of Arris, see 
Eq. PHj) in Sec. 13. 3. 2^ . the best knowledge of the lifetimes (see Sec. 13. 2p and the estimate of 
/baryon giveu abovc, yiclds fs = 0.120±0.021, an estimate dominated by the mixing information. 
Taking into account all known correlations (including the one introduced by /baryon), this result 
is then combined with the set of fractions obtained from direct measurements (given above), to 
yield the improved estimates of Table |21 still under the constraints of Eq. (fT^ . As can be seen, 
our knowledge on the mixing parameters substantially reduces the uncertainty on fs, despite 
the rather strong deweighting introduced in the computation of the world average of x- It 
should be noted that the results are correlated, as indicated in Table |21 
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3.2 6-hadron lifetimes 



In the spectator model the decay of 6-flavored hadrons is governed entirely by the flavor 
changing b — ^ Wq transition {q = c,u). For this very reason, lifetimes of all 6-flavored hadrons 
are the same in the spectator approximation regardless of the (spectator) quark content of the 
Hf,. In the early 1990's experiments became sophisticated enough to start seeing the differences 
of the lifetimes among various Hb species. The first theoretical calculations of the spectator 
quark effects on Hb lifetime emerged only few years earlier. 

Currently, most of such calculations are performed in the framework of the Heavy Quark 
Expansion, HQE. In the HQE, under certain assumptions (most important of which is that of 
quark-hadron duality), the decay rate of an Hb to an inclusive final state / is expressed as the 
sum of a series of expectation values of operators of increasing dimension, multiplied by the 
correspondingly higher powers of AQCYy/nib- 

r^,_, = \CKM\'J2clf^(^y{Hb\0^\Hb), (18) 

^ / / vb ' 

n 

where \CKM\'^ is the relevant combination of the CKM matrix elements. Coefficients Cn^ of 
this expansion, known as Operator Product Expansion |221, can be calculated perturbatively. 
Hence, the HQE predicts Th^^j in the form of an expansion in both AQCD/i^b and asirrib). The 
precision of current experiments makes it mandatory to go to the next-to-leading order in QCD, 
i.e. to include correction of the order of as{mb) to the Cn^^s. All non-perturbative physics is 
shifted into the expectation values {Hb\On\Hb) of operators On- These can be calculated using 
lattice QCD or QCD sum rules, or can be related to other observables via the HQE [22] • One 
may reasonably expect that powers of yiqco/'^b provide enough suppression that only the first 
few terms of the sum in Eq. fjl8p matter. 

Theoretical predictions are usually made for the ratios of the lifetimes (with t{B^) chosen 
as the common denominator) rather than for the individual lifetimes, for this allows several 
calculational uncertainties to cancel. Present HQE results including next-to-leading order cor- 
rections (see Refs. [2111231201 for the latest updates) are in some instances already surpassed by 
the experimental measurements, e.g. in the case of t{B~^)/t{B^). Also, HQE calculations are 
not assumption-free. More accurate predictions are a matter of progress in the evaluation of the 
non-perturbative hadronic matrix elements and verifying the assumptions that the calculations 
are based upon. However, the HQE, even in its present shape, draws a number of important 
conclusions, which are in agreement with experimental observations: 

• The heavier the mass of the heavy quark the smaller is the variation in the lifetimes among 
different hadrons containing this quark, which is to say that as rrih ^ oo we retrieve the 
spectator picture in which the lifetimes of all HbS are the same. This is well illustrated by 
the fact that lifetimes in the b sector are all very similar, while in the c sector [rric < rub) 
lifetimes differ by as much factor of 2. 

• The non-perturbative corrections arise only at the order of yigQ^/m^, which translates 
into differences among Hb lifetimes of only a few percent. 

• It is only the difference between meson and baryon lifetimes that appears at the yiQc^/m^ 
level. The splitting of the meson lifetimes occurs at the ^QCD/'^b level, yet it is enhanced 
by a phase space factor IGtt^ with respect to the leading free b decay. 
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To ensure that certain sources of systematic uncertainty cancel, lifetime analyses are some- 
times designed to measure a ratio of lifetimes. However, because of the differences in decay 
topologies, abundance (or lack thereof) of decays of a certain kind, etc., measurements of the in- 
dividual lifetimes are more common. In the following section we review the most common types 
of the lifetime measurements. This discussion is followed by the presentation of the averaging 
of the various lifetime measurements, each with a brief description of its particularities. 

3.2.1 Lifetime measurements, uncertainties and correlations 

In most cases lifetime of an Hh is estimated from a flight distance and a /3'j factor which is used 
to convert the geometrical distance into the proper decay time. Methods of accessing lifetime 
information can roughly be divided in the following five categories: 

1. Inclusive (flavor blind) measurements. These measurements are aimed at extract- 
ing the lifetime from a mixture of 6-hadron decays, without distinguishing the decaying 
species. Often the knowledge of the mixture composition is limited, which makes these 
measurements experiment or accelerator specific. Also, these measurements have to rely 
on Monte Carlo for estimating the (3'-/ factor, because the decaying hadrons are not fully 
reconstructed. On the bright side, these usually are the largest statistics 6-hadron lifetime 
measurements that are accessible to a given experiment, and can, therefore, serve as an 
important performance benchmark. 

2. Measurements in semileptonic decays of a specific Hb- W from h Wc pro- 
duces iui pair {i = e, /x) in about 21% of the cases. Electron or muon from such decays is 
usually a well-detected signature, which provides for clean and efficient trigger, c quark 
from b — s> Wc transition and the other quark(s) making up the decaying Hb combine into 
a charm hadron, which is reconstructed in one or more exclusive decay channels. Know- 
ing what this charmed hadron is allows one to separate, at least statistically, different Hh 
species. The advantage of these measurements is in statistics which is usually superior 
to that of the exclusively reconstructed Hj, decays. Some of the main disadvantages are 
related to the difficulty of estimating lepton-|-charm sample composition and Monte Carlo 
reliance for the (3'-/ factor estimate. 

3. Measurements in exclusively reconstructed decays. These have the advantage of 
complete reconstruction of decaying Hb, which allows one to infer the decaying species as 
well as to perform precise measurement of the /57 factor. Both lead to generally smaller 
systematic uncertainties than in the above two categories. The downsides are smaller 
branching ratios, larger combinatoric backgrounds, especially in Hb — >■ Hc'k{'ktt) and 
multi-body He decays, or in a hadron collider environment with non-trivial underlying 
event. Hb — » J/ipHs are relatively clean and easy to trigger on J/ip but their 
branching fraction is only about 1%. 

4. Measurements at asymmetric B factories. In the T{AS) BB decay, the B 
mesons {B^ or B^) are essentially at rest in the T(4S') rest frame. This makes lifetime 
measurements impossible with experiments, such as CLEO, in which T(45') produced at 
rest. At asymmetric B factories T(4S') is boosted resulting in B and B moving nearly 
parallel to each other. The lifetime is inferred from the distance A2; separating B and 
B decay vertices and T(4S') boost known from colliding beam energies. In order to 
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maximize the precision of the measurement one B meson is reconstructed in the D^*Hh'e 
decay. The other B is typically not fully reconstructed, only position of its decay vertex is 
determined. These measurements benefit from very large statistics, but suffer from poor 
Az resolution. 

5. Direct measurement of lifetime ratios. This method has so far been only applied 
in the measurement of t{B^)/t{B^). The ratio of the lifetimes is extracted from the 
dependence of the observed relative number of B^ and B^ candidates (both reconstructed 
in semileptonic decays) on the proper decay time. 

In some of the latest analyses, measurements of two {e.g. t{B^) and t{B^)/t{B^)) or three 
{e.g. t{B'^), t{B^)/t{B^), and Am^) quantities are combined. This introduces correlations 
among measurements. Another source of correlations among the measurements are the sys- 
tematic effects, which could be common to an experiment or to an analysis technique across 
the experiments. When calculating the averages, such correlations are taken into account per 
general procedure, described in Ref. |77j . 

3.2.2 Inclusive b-hadron lifetimes 

The inclusive b hadron lifetime is defined as Tb = Yli fi'^i where Ti are the individual 6-hadron 
lifetimes and fi the fractions of the various species present in an unbiased sample of weakly- 
decaying b hadrons produced at a high-energy collider.^ This quantity is certainly less funda- 
mental than the lifetimes of the individual b hadron species, the latter being much more useful 
in the comparison of measurements with theoretical predictions. Nonetheless, we present the 
measurements of the inclusive lifetime for completeness. 

In practice, an unbiased measurement of the inclusive lifetime is difficult to achieve, because 
it would imply an efficiency which is guaranteed to be the same across species. So most of the 
measurements are biased. In an attempt to group analyses which are expected to select the 
same mixture of b hadrons, the available results (given in Table Ej) are divided into the following 
three sets: 

1. measurements at LEP and SLD that accept any 6-hadron decay, based on topological 
reconstruction (secondary vertex or impact parameters) using charged tracks; 

2. measurements at LEP based on the identification of a lepton from a b decay; and 

3. measurements at the Tevatron based on inclusive b ^ J/tp reconstruction, where the J/ip 
is fully reconstructed. 

The measurements of the first set are generally considered as estimates of Tb, although the 
efficiency to reconstruct a secondary vertex most probably depends, in an analysis-dependent 
way, on the number charged tracks from a b hadron decay, which in turn should depend on the 
type of b hadron. Even though these efficiency variations can in principle be accounted for using 
Monte Carlo simulations (which inevitably contain assumptions on branching fractions), the 
b hadron mixture in that case can remain somewhat ill-defined and could be slightly different 
between analyses in that set. 

^In principle such a quantity could be slightly different in Z decays and a the Tevatron, in case the fractions 
of 6-hadron species are not exactly the same; see the discussion in Sec. 13.1.21 
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Table 3: Measurements of average 6-hadron lifetimes. 



Experiment 


Method 


Data set 


n (ps) 


Ref. 


ALEPH 


Dipole 


91 


1.511 ±0.022 ±0.078 


I2H! 


DELPHI 


All track i.p. (2D) 


91-92 


1.542 ±0.021 ±0.045 




DELPHI 


Sec. vtx 


91-93 


1.582 ±0.011 ±0.027 




DELPHI 


Sec. vtx 


94-95 


1.570 ±0.005 ±0.008 


m 


L3 


Sec. vtx + i.p. 


91-94 


1.556 ±0.010 ±0.017 


OPAL 


Sec. vtx 


91-94 


1.611 ± 0.010 ± 0.027 




SLD 


Sec. vtx 


93 


1.564 ± 0.030 ± 0.036 


M 


Average set 1 [b vertex) 




i.O (0 It U.UUo 




ALEPH 


Lepton i.p. (3D) 


yi— yo 


i.ooo it U.Uio it u.uzz 


I'J CI 


L3 


Lepton i.p. (2D) 


91-94 


1.544 ±0.016 ±0.021 


1221'' 


OPAL 


Lepton i.p. (2D) 


90-91 


1.523 ±0.034 ±0.038 


isni 


Average set 2 (6 — > ^) 




1.537 ±0.020 




CDF 


J/ip vtx 


92-95 


1.533 ±0.015+^:^;^^ 


m 


Average of all above 




1.574 ±0.008 





The combined DELPHI result quoted in [201 is 1.575 ± 0.010 ± 0.026 ps. 
^ The combined L3 resuh quoted in ^ is 1.549 ± 0.009 ± 0.015 ps. 



On the contrary, the mixtures corresponding to the other two sets of measurements are 
better defined, in the limit where the reconstruction and selection efficiency of a lepton or a 
J /ip from a b hadron does not depend on the type of that b hadron. These mixtures are given 
by the production fractions and the inclusive branching fractions for each b hadron species 
to give a lepton or a J/ip. In particular, under the assumption that all b hadrons have the 
same semileptonic decay width, the analyses of the second set should measure T{b ^ i) = 
(Yli fi^i) / iYli fi'^i) which is necessarily larger than Tf, if lifetime differences exist. Given the 
present knowledge on Tj and r(6 —>•£) — is expected to be of the order of 0.01 ps. 

For the averaging, correlated systematic are taken into account, which are due to b and c 
fragmentation, b and c decay models, B{B — > i), B{B ^ c ^ i), Tc, B{c — > i), and B charged 
track decay multiplicity. The averages for the sets defined above (also given in Table |21) are 

r(6 vertex) = 1.576 ± 0.008 ps , (19) 
T{b^£) = 1.537 ± 0.020 ps, (20) 
r(6^J/V0 = 1.533l°;°ips, (21) 

whereas an average of all measurements, ignoring mixture differences, yields 1.574 ± 0.008 ps. 
3.2.3 and B+ lifetimes 

In LEP and CDF experiments, the most precise measurements of the B^ and B~^ lifetimes 
have originated from two classes of partially reconstructed decays. In the first class the decay 
B ^ D I'^v^X is used in which the charge of the charmed meson distinguishes between neutral 
and charged B mesons; the ratio B^ / B^ lifetime ratio is usually also extracted directly, as the 
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Table 4: Measurements of the lifetime. 



Experiment 


Method 


Data set 


r(fi°) (ps) 


Ref. 


ALEPH 


D^*H 


91- 


-95 


1.518 ±0.053 ±0.034 


m 


ALEPH 


Exclusive 


91- 


-94 


1.25+°:^^ ±0.05 


m 


ALEPH 


Partial rec. vt+tt^ 


91- 


-94 


1 /IQ+O.17+0.08 
-■-•^^-0. 15-0. 06 


m 


CDF 


Di*)i 


92- 


-95 


1.474 ±0.039l[];[!^? 


m 


CDF 


Excl. J/ijK 


92- 


-95 


1.497 ±0.073 ±0.032 


m 


CDF 


Excl. J/ipK 


02- 


-03 


1.49 ±0.06 ±0.062 




DELPHI 




91- 


-93 


1.6ll{];J^ ± 0.08 


m 


DELPHI 


Charge sec. vtx 


91- 


-93 


1.63 ±0.14 ±0.13 


m 


DELPHI 


Inclusive D*i 


91- 


-93 


1.532 ±0.041 ±0.040 


m 


DELPHI 


Charge sec. vtx 


94- 


-95 


1.531 ±0.021 ±0.031 


m 


L3 


Charge sec. vtx 


94- 


-95 


1.52 ±0.06 ±0.04 


m 


OPAL 




91- 


93 


1.53 ±0.12 ±0.08 


m 


OPAL 


Charge sec. vtx 


93- 


-95 


1.523 ±0.057 ±0.053 


m 


OPAL 


Inclusive D*i 


91- 


-00 


1.541 ±0.028 ±0.023 


m 


SLD 


Charge sec. vtx i 


93- 


-95 


1.56l°:i^±0.10 




SLD 


Charge sec. vtx 


93- 


-95 


1.66 ±0.08 ±0.08 


ISO]" 


BABAR 


Exclusive 


99- 


-00 


1.546 ±0.032 ±0.022 


M 


BABAR 


Inclusive D*i 


99- 


-01 


1.529 ±0.012 ±0.029 




BABAR 


Exclusive D*i 


99- 


-02 


1.523l[]:[!^^± 0.022 


|SS! 


BABAR 


Inch D*TT, D*p 


99- 


-01 


1.533 ±0.034 ±0.038 


P! 


Belle 


Exclusive 


00- 


-01 


1.554 ±0.030 ±0.019 




Average 








1.534 ±0.013 





Recent measurements not yet included in the average 
CDF Excl. J/ipK 02-04 1.539 ± 0.051 ± 0.008 §^ 
BABAR Inclusive 99-04 1.501 ± 0.008 ± 0.030 jST] 

Belle Exclusive 00-03 1.534 ± 0.008 ± 0.010 

" The combined SLD result quoted in ^ is 1.64 ± 0.08 ± 0.08 ps. 
^ To replace [HHl- To replace |nS]. ^ Preliminary. 



B^ and B~^ lifetime measurements are correlated. In the second class the charge attached to 
the 6-decay vertex is used to achieve this separation. 

With the benefit from very large statistics, the asymmetric i?-factory experiments, BABAR 
and Belle, now provide more precise lifetime measurements using exclusively reconstructed de- 
cays, as well as partially reconstructed decays. With increased data sample, CDF also provides 
improved measurements using large samples of exclusive B^ J/ipK^*^^ and B^ JipK^ 
decays. 

The averaging is summarized in Tables |31 13 and IHl The following sources of correlated 
systematic uncertainties have been considered: D** branching ratio uncertainties j2], mo- 
mentum estimation of B mesons from decays (6-quark fragmentation parameter (Xe) = 
0.702 ± 0.008 and 6-baryon lifetimes (see Sees. 13.2.41 and I3.2.6|) . and 6-hadron fractions 
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Table 5: Measurements of the 2?+ lifetime. 



Experiment 


Method 


Data set 


r{B+) (ps) 


Ref. 


ALEPH 




91- 


-95 


1.648 ± 0.049 ± 0.035 


m 
m 


ALEPH 


Exclusive 


91- 


-94 


1 c;o+0.21+0.04 
-■-•^O-O.IS-O.OS 


CDF 




92- 


-95 


1.637 ±0.058l[5:[;^^ 


m 


CDF 


Excl. J/tpK 


92- 


-95 


1.636 ± 0.058 ± 0.025 


m 


CDF 


Excl. (J/ipK) 


02- 


-03 


1.64 ±0.05 ±0.02 




DELPHI 


£,(*)£ 


91- 


93 


1.61 ±0.16 ±0.12 




DELPHI 


Charge sec. vtx 


91- 


-93 


1.72 ±0.08 ±0.06 


im" 


DELPHI 


Charge sec. vtx 


94- 


-95 


1.624 ±0.014 ±0.018 


m 


L3 


Charge sec. vtx 


94- 


-95 


1.66 ±0.06 ±0.03 


m 


OPAL 




91- 


-93 


1.52 ±0.14 ±0.09 


m 


OPAL 


Charge sec. vtx 


93 


-95 


1.643 ±0.037 ±0.025 


gHj 


SLD 


Charge sec. vtx i 


93- 


-95 


1.611^;}^ ±0.07 




SLD 


Charge sec. vtx 


93- 


-95 


1.67 ±0.07 ±0.06 


BABAR 


Exclusive 


99- 


-00 


1.673 ±0.032 ±0.023 


EH 


Belle 


Exclusive 


00- 


-01 


1.695 ±0.026 ±0.015 




Average 








1.652 ±0.014 





Recent measurements not yet included in the average 
CDF Excl. J/^jK 02-04 1.662 ± 0.033 ± 0.008 
Belle Exclusive 00-03 1.635 ± 0.011 ± 0.011 

The combined DELPHI result quoted in is 1.70 ± 0.09 ps. 
The combined SLD result quoted in ^ is L66 ± 0.06 ± 0.05 ps. 
To replace [HHI- To replace ^ Preliminary. 



at high energy (see Table |21). The world averages are: 

r(fi°) = 1.534 ± 0.013 ps, (22) 
r(fi+) = 1.653 ± 0.014 ps, (23) 
r(5+)/r(5°) = 1.081 ±0.015. (24) 

As indicated in the tables, these averages do not include recent CDF jSH] , D0 jSH] , .BABAR [SZ] , 
and Belle j^Hl results; it is planned to incoporate them in the next update of this report. 

3.2.4 lifetime 

Some of the more precise measurements of the B^ lifetime originate from partially reconstructed 
decays in which a Dj meson has been completely reconstructed. 

The following correlated systematic errors were considered: average B lifetime used in 
backgrounds, B^ decay multiplicity, and branching ratios used to determine backgrounds {e.g. 
B{B ^ DsD)). 

A knowledge of the multiplicity of B^ decays is important for measurements that partially 
reconstruct the final state such as i? — > D^X (where X is not a lepton). The boost deduced 
from Monte Carlo simulation depends on the multiplicity used. Since this is not well known. 
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Table 6: Measurements of the ratio t[B^)/t{B^). 



Experiment 


Method 


Data set 


Ratio t{B+)/t{B°) 


Ref. 


ALEPH 




91- 


-95 


1.085 ±0.059 ±0.018 


m 


ALEPH 


Exclusive 


91- 


-94 




1 97+0.23+0.03 
-■-•^ ' -0.19-0.02 




CDF 




92- 


-95 




1.110 ±0.056+^30 




CDF 


Excl. J/^K 


92- 


-95 


1 


.093 ±0.066 ±0.028 


m 


DELPHI 




91- 


-93 




1.00t|];J^±0.10 




DELPHI 


Charge sec. vtx 


91- 


-93 




1.06l°:J?±0.10 


mi 


DELPHI 


Charge sec. vtx 


94- 


-95 


1, 


.060 ± 0.021 ± 0.024 


EH 


L3 


Charge sec. vtx 


94- 


-95 




1.09 ±0.07 ±0.03 


m 


OPAL 




91- 


-93 




0.99±0.14t°;°| 


m 


OPAL 


Charge sec. vtx 


93- 


-95 


1 


.079 ±0.064 ±0.041 




SLD 


Charge sec. vtx i 


93- 


-95 




1.03+g-J^ ± 0.09 




SLD 


Charge sec. vtx 


93- 


-95 




imt^f, ± 0.05 




BABAR 


Exclusive 


99- 


-00 


1, 


.082 ±0.026 ±0.012 


EH 


Belle 


Exclusive 


00- 


-01 


1 


.091 ±0.023 ±0.014 


m 


Average 










1.081 ±0.015 





Recent measurements not yet included in the average 



D0 D*+i2 ratio 02-04 1.080 ± 0.016 ± 0.014 0~ 
CDF Excl. J/^pK 02-04 1.080 ± 0.042 

Belle Exclusive 00-03 1.066 ± 0.008 ± 0.008 

" The combined SLD result quoted in ^ is 1.01 ± 0.07 ± 0.06. 
To replace |55j . ^ Preliminary. 



the multiplicity in the simulation is varied and this range of values observed is taken to be a 
systematic. 

Similarly not all the branching ratios for the potential background processes are measured. 
Where they are available, the PDG values are used for the error estimate. Where no mea- 
surements are available estimates can usually be made by using measured branching ratios of 
related processes and using some reasonable extrapolation. 

The inputs used to form the average B^ lifetime are given in Table [7| 

It is important to note that similar to the kaon system, neutral B mesons contain short- 
and long-lived components, since the two mass eigenstates B^ and Bh differ not only in their 
masses, but also in their widths with AP = Tl — Th- In the Standard Model for the B^ system, 
this difference can be large, i.e., ATs/Tg = 0.12 ± 0.06 jHSl- Specific measurements of this 
parameter are explained in more detail in Sec. 13.3.21 

Flavor-specific decays, such as semileptonic Bg Dgiu will have equal contributions of 
th = and = I/Pl- For this reason, lifetime measurements via these decays are broken 

out separately as given in Table [71 and their world average is: 

t{BI)d,i = 1.442 ± 0.066 ps . (25) 

Any final state can be decomposed into its CP-even and CP-odd component, and in the 
Standard Model, AP^cp = AP^. Fully exclusive decays of P° into P° — J/tlxj) are expected to 
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Table 7: Measurements of the lifetime. 



Experiment 


Method 


Data set 


r{B'g) (ps) 


Ref. 


ALEPH 


Dsl 


91-95 


1.54l^:i| ± 0.04 


m 


CDF 


Dsi 


92-96 


1.36 ± 0.09tl f5 


m 


DELPHI 


DJ 


91-95 


lA2toit ± 0.03 


m 


OPAL 


DJ 


90-95 


1.50lgi^ ± 0.04 


m 


Average of Dgl measurements 


1.442 ±0.066 




ALEPH 


Dsh 


91-95 


1.47 ±0.14 ±0.08 


m 


DELPHI 


Dsh 


91-95 


1.53l|]:J^ ±0.07 


m 


DELPHI 


Dg inch 


91-94 






OPAL 


Dg inch 


90-95 


1 79+0.20+0.18 
^-0.19-0.17 


inn 


Average of all above Dg measurements 


1.469 ±0.059 




CDF 




92-95 


1.34+^:^^ ± 0.05 


EH 


CDF 




02-04 


1.369 ±o.iool|]:[J?^ 




D0 




02-04 


1.4441°;°^^ ± 0.02 


m 



Average of J/ip4> measurements 1.404 ± 0.066 



Preliminary. 



be dominated by the CP-even state and its lifetime. First measurements of the CP mix for this 
decay mode are outlined in Sec. 13.3.21 CDF and D0 measurements from this particular mode 
— »• J/ ipcf) are combined into an average given in Table [7| There are no correlations between 
the measurements for this fully exclusive channel, and the world average for this specific decay 
is: 

r(E°) j/^,^ = 1.404 ± 0.066 ps . (26) 

Finally, the remaining measurements are variations of lifetimes determined more inclusively 
via Dg plus hadrons, and hence into a more unknown mixture of flavors and/or CP-states. A 
lifetime weighted this way can still be a useful input for analyses examining such an inclusive 
sample. These are separated in Tableland combined with the semileptonic lifetime to obtain: 

r(5°) = 1.469 ± 0.059 ps . (27) 

3.2.5 B+ lifetime 

c 

There are currently two measurements of the lifetime of the B'^ meson from CDF (70] and 
D0 [71] using the semileptonic decay mode B^ — >■ J/ ipi and fitting simultaneously to the mass 
and lifetime using the vertex formed from the leptons from the decay of the J/ ip and the third 
lepton. Correction factors to estimate the boost due to the missing neutrino are used. Mass 
values of 6.40 ± 0.39 ± 0.13 GeV/c^ and 5.951°;}^ ± 0.34 GeV/c^, respectively, are found by 
fitting to the tri-lepton invariant mass spectrum. These mass measurements are consistent 
to within errors, and no adjustments to the lifetimes are made. Correlated systematic errors 
include the impact of the uncertainty of the B^ px spectrum on the correction factors, the 
level of feed-down from %l)[2S), MC modeling of the decay model varying from phase space to 
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Table 8: Measurements of the lifetime. 



Experiment 


Method 


Data set 


riB+) (ps) 


Ref. 


CDF 
D0 


J/ipfi 


92-95 
02-04 


0.46+^:i«± 0.03 
0.448t!|:j^^± 0.121 


[70] 


Average 






0.45 ± 0.12 





Preliminary. 



the ISGW model, and mass variations. Inputs are given in Table |H1 and the world average is 
determined to be: 

r(fi+) = 0.45 ±0.12 ps. (28) 
3.2.6 and b-baryon lifetimes 

The most precise measurements of the 6-baryon lifetime originate from two classes of partially 
reconstructed decays. In the first class, decays with an exclusively reconstructed A'^ baryon 
and a lepton of opposite charge are used. These products are more likely to occur in the decay 
of Al baryons. In the second class, more inclusive final states with a baryon (p, p, A, or A) and 
a lepton have been used, and these final states can generally arise from any b baryon. 

The following sources of correlated systematic uncertainties have been considered: exper- 
imental time resolution within a given experiment, 6-quark fragmentation distribution into 
weakly decaying b baryons, A^ polarization, decay model, and evaluation of the 6-baryon purity 
in the selected event samples. In computing the averages the central values of the masses are 
scaled to M(yl°) = 5624 ± 9 MeV/c^ [72] and M(&-baryon) = 5670 ± 100 MeV/c^. 

The meaning of decay model and the correlations are not always clear. Uncertainties related 
to the decay model are dominated by assumptions on the fraction of n-body decays. To be 
conservative it is assumed that it is correlated whenever given as an error. DELPHI varies the 
fraction of 4-body decays from 0.0 to 0.3. In computing the average, the DELPHI result is 
corrected for 0.2 ± 0.2. 

Furthermore, in computing the average, the semileptonic decay results are corrected for a 
polarization of — 0.451q;}7 [2] and a A^ fragmentation parameter (Xe) = 0.70 ± 0.03 [75] . 

Inputs to the averages are given in Table |01 The world average lifetime of b baryons is then: 

(r(6-baryon)) = 1.210 ± 0.048 ps . (29) 

Keeping only Afi"^ and Ai~i~^ final states, as representative of the A^ baryon, the following 
lifetime is obtained: 

r(/l^) = 1.232 ±0.072 ps. (30) 

Averaging the measurements based on the S^i^ final states fS] gives a lifetime value 
for a sample of events containing and baryons: 

(r(-,)) = 1.391°:^^ ps. (31) 
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Table 9: Measurements of the 6-baryon lifetimes. 



Experiment 


Method 


Data set 


Lifetime (ps) 


Ref. 


ALEPH 




91- 


-95 


1 1 8+0.13 


±0.03 


m 


ALEPH 




91- 


-95 




±0.04 


m 


CDF 




91- 


-95 


1.32 ±0.15 ±0.06 


m 


CDF 


j/i)A 


02- 


-03 


1.25 ±0.26 ±0.10 




D0 


J/lpA 


02- 


-04 


1 99+0.22 
-■-•^^-0.18 


±0.04 


m 


DELPHI 


Ati 


91- 


-94 


1 1 1 +0.19 
-'-•-'--'--0.18 


±0.05 


[771'' 


OPAL 


At I, Al-i+ 


90- 


-95 


-'-•^^-0.22 


±0.06 


m 


Average of ab 


ove 7 [Al lifet: 


ime) 




1.232 ± 


0.072 




ALEPH 


Al 


91- 


-95 


-'-•^'-'-0.08 


±0.06 


m 


DELPHI 


AiTl vtx 


91- 


-94 


1.16 ±0.20 ±0.08 




DELPHI 


Afi i.p. 


91- 


-94 


1 in+oi9 

^■^^-0.17 


±0.09 




DELPHI 


pi 


91- 


-94 


1.19 ±0.14 ±0.07 


[77j° 


OPAL 


A£ i.p. 


90- 


-94 


1 91+0.15 
-'-•^-'--0.13 


±0.10 




OPAL 


Ai vtx 


90- 


-94 


1.15 ±0.12 ±0.06 




Average of ab 


ove 13 (6-baryon lifetime) 


1.210 ±0.048 




ALEPH 
DELPHI 


Si 


90- 
91- 


-95 
-93 


1 qc+O.37+0.15 
±.00_Q 28-0.17 

1.5t°i±0.3 


m 


Average of ab 


ove 2 {Ef, lifet: 


ime) 




1 30+0.34 

i.OJ_Q 28 





The combined DELPHI result quoted in |771 is 1.14+0.08 ± 0.04 ps. 
The combined OPAL result quoted in ^ is 1.16 ± 0.11 ± 0.06 ps. 



Preliminary. 



3.2.7 Summary and comparison to theoretical predictions 

Averages of lifetimes of specific b hadron species are collected below in Table E3 

As described in Sec. 13.21 Heavy Quark Effective Theory can be employed to explain the 
hierarchy of T(i?+) ^ T(yl°) < T(i?°) ^ t{B^) < t{B^), and used to predict the ratios between 
lifetimes. A recent prediction of the ratio between the and 5° lifetimes, is 1.06 ± 0.02 [21]. 
The ratio r(yl^)/r(-B°) has particularly been the source of theoretical scrutiny since earlier 
calculations jHO] predicted a value greater than 0.90, almost two sigma higher than the world 
average at the time. Recent calculations of this ratio that include higher order effects predict a 
ratio between the A^ and lifetimes of 0.86 ± 0.05 |21] and reduces this difference. Ref. 
presents probability density functions of its predictions with variation of theoretical inputs, and 
the indicated errors (and ranges in Table ITD below) are the RMS of the distributions. Measured 
lifetime ratios compared to predicted ranges are given in Table ^2 

3.3 Neutral S-meson mixing 

There are two neutral B — B systems, B'^ — B^ and B^^ — B^^, which both exhibit the phenomenon 
of particle-antiparticle mixing. For each of these systems, there are two mass eigenstates which 
are linear combinations of the two flavour states, B ot B. We consider the case where a neutral 
B meson is produced and detected in a flavour state, through its decay to a flavour-specific 



24 



Table 10: Summary of lifetimes of different b hadron species. 



b hadron species Measured lifetime 



B+ 


1.653 ±0.014 ps 




1.534 ±0.013 ps 


5° (— s> flavor specific) 


1.442 ±0.066 ps 


5° J/^0) 


1.404 ±0.066 ps 




0.45 ±0.12 ps 


At 


1.232 ±0.072 ps 


Sh mixture 


1.39l^Q 28 


6-baryon mixture 


1.210 ±0.048 ps 


6-hadron mixture 


1.574 ±0.008 ps 



Table 11: Ratios of 6-hadron lifetimes relative to the lifetime and theoretical ranges predicted 
by theory 



Lifetime ratio 


Measured value 


Predicted range 


r{B+)/r{B') 


1.081 ±0.015 


1.04 - 1.08 


T{B'i)/T{By 


0.939 ± 0.044 


0.99 - 1.01 


T{Al)/r{B') 


0.803 ±0.047 


0.81 - 0.91 


r(6-baryon)/r(fi°) 


0.789 ±0.032 


0.81 - 0.91 



Using the flavor specific lifetime for definiteness. 



final state. There are four different time-dependent probabilities; if CPT is conserved (which 
will be assumed throughout), they can be written as 



' V{B - 


-B) = 


^ [cosh(f t) + 


cos(Amt)] 


2 


V{B - 

< 


-B) = 


^ [cosh(f t) - 


cos(Amt)] 


£ 

P 


ViB - 


-B) = 


[cosh(f t) - 


cos(Amt)] 


2 

P 
<? 




-B) = 


^ [cosh(f t) ± 


cos(Amt)] 





(32) 



where t is the proper time of the system {i.e. the time interval between the production and 
the decay in the rest frame of the B meson) and F = 1/t{B) is the average decay width. At 
the B factories, only the proper-time difference At between the decays of the two neutral B 
mesons from the T(45') can be determined, but, because the two B mesons evolve coherently 
(keeping opposite flavours as long as none of them has decayed), the above formulae remain 
valid if t is replaced with At and the production flavour is replaced by the flavour at the time 
of the decay of the accompanying B meson in a flavour specific state. As can be seen in the 
above expressions, the mixing probabilities depend on the following three observables: the mass 
difference Am and the decay width difference AF between the two mass eigenstates, and the 
parameter which signals CP violation in the mixing if \q/p\'^ 7^ 1. 

In the following sections we review in turn the experimental knowledge on these three 
parameters, separately for the 5° meson (Am^, AF^, \q/p\d) and the B^ meson (Am^, AF^, 
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\(1/P\s)- 



3.3.1 mixing parameters 



CP violation parameter \q/p\d 



Evidence for CP violation in 5° mixing has been searched for, both with flavor-specific 
and inclusive decays, in samples where the initial flavor state is tagged. In the case of 
semileptonic (or other flavor-specific) decays, where the final state tag is also available, the 
following asymmetry 



A 



SL 



N{B (t) ^ f+z/,X) - N{B\t) ^ t-V,X) _ \p/q\[ 



\p/(l\l- 



\<l/p\l 



(33) 



has been measured, either in time- integrated analyses at CLEO [S2 IH21 EH] and CDF [HI] , or 
in time-dependent analyses at OPAL jHS], ALEPH BABAR [HZl Ei and Belle |Hi. In the 
inclusive case, also investigated and published at ALEPH (HH] and OPAL [^0], no final state 
tag is used, and the asymmetry [21] 



N{B^{t) all) - N(W{t) all) 



N{B^{t) 



all) + N{B\t) 



all) 



A 



SL 



2r, 



sm(Amd t) — sin^ 



Am^ t 



(34) 



must be measured as a function of the proper time to extract information on CP violation. In 
all cases asymmetries compatible with zero have been found, with a precision limited by the 
available statistics. A simple average of all published results for the B^ meson 
I57| 1551 EO] and of the preliminary Belle result [55] yields 



= -0.0026 ± 0.0067 



or, equivalently through Eq. (jHSl), 



\<1/P\c 



1.0013 ±0.0034. 



(35) 



(36) 



This result'^, summarized in Table IT^ is compatible with no CP violation in the mixing, an 
assumption we make for the rest of this section. 

Mass and decay width differences Am^ and AF^ 

Many time-dependent B^-^ oscillation analyses have been published by the ALEPH, 
BABAR, Belle, CDF, DELPHI, L3 and OPAL collaborations. The corresponding measure- 
ments of Arrid are summarized in Tables and ^3 where only the most recent results are 
listed {i.e. measurements superseded by more recent ones have been omitted). Although a 
variety of different techniques have been used, the individual Ama results obtained at high- 
energy colliders have remarkably similar precision. Their average is compatible with the recent 



^Early analyses and (perhaps hence) the PDG use the complex parameter es = {p — q)/{p+q)', if CP violation 
in the mixing in small, .4sl — 4Re(eB)/(l + |e_Bp) and our current world average is Re(eB)/(l + |esp) = 
-0.0007 ±0.0017. 
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Table 12: Measurements of CP violation in mixing and their average in terms of both 
and \q/p\d- The individual results are listed as quoted in the original publications, or converted^ 
to an value. When two errors are quoted, the first one is statistical and the second one 
systematic. 



Exp. & Ref. 


Method 




Measured ^sl 


Measured \q/p\d 


CLEO 




partial hadronic rec. 


+0, 


.017 ±0.070 ±0.014 




CLEO 


M 


dileptons 


+0, 


.013 ±0.050 ±0.005 




CLEO 




average of above two 


+0, 


.014 ±0.041 ±0.006 




OPAL 


IBS] 


leptons 


+0, 


.008 ±0.028 ±0.012 




OPAL 


jno] 


inclusive (Eq. (jSlI)) 


+0, 


.005 ±0.055 ±0.013 




ALEPH 


M 


leptons 


-0, 


.037 ±0.032 ±0.007 




ALEPH 


M 


inclusive (Eq. (jMI)) 


+0, 


.016 ±0.034 ±0.009 




ALEPH 


M 


average of above two 


-0, 


.013 ± 0.026 (tot) 




BABAR 


|HH] 


dileptons 


+0, 


.005 ±0.012 ±0.014 


0.998 ±0.006 ±0.007 


BABAR 




full hadronic rec. 






1.029 ±0.013 ±0.011 


Belle 


m 


dileptons (prel.) 


-0, 


.0013±0.0060±0.0056 


1.0006 ±0.0030±0.0028 



Average of all above -0.0026 ± 0.0067 (tot) 1.0013 ± 0.0034 (tot) 



and more precise measurements from the asymmetric B factories. The systematic uncertain- 
ties are not negligible; they are often dominated by sample composition, mistag probability, 
or 6-hadron lifetime contributions. Before being combined, the measurements are adjusted on 
the basis of a common set of input values, including the averages of the 6-hadron fractions 
and lifetimes given in this report (see Sees. 13.11 and 13. 2|) . Some measurements are statisti- 
cally correlated. Systematic correlations arise both from common physics sources (fractions, 
lifetimes, branching ratios of b hadrons), and from purely experimental or algorithmic effects 
(efficiency, resolution, flavour tagging, background description). Combining all published mea- 
surements listed in Table ^1 and accounting for all identifled correlations as described in [2] 
yields Am^ = 0.502 ± 0.004 ± 0.005 ps'^ 

On the other hand, ARGUS and CLEO have published measurements of the time-integrated 
mixing probability Xd [II21 IHH IH2] , which average to Xd = 0.182±0.015. Following Ref. [H2!, the 
width difference AF^ could in principle be extracted from the measured value of F^ = l/r(i?°) 
and the above averages for Am^ and Xd (provided that AF^ has a negligible impact on the 
Amd analyses that have assumed AF^ = 0), using the relation 

^d + Vd ... ^^d , AFrf 

^^ = 2MTi) -^h., = ^ ^"^^^=2f;- ^'^^ 

However, direct time-dependent studies yield stronger constraints: DELPHI published the re- 
sult \ATd\/Td < 18% at 95% CL (OH, while BABAR recently obtained -8.4% < sign(ReAcp) AFd/Fd < 
6.8% at 90% CL [HZ]. 

Assuming AF^ = and using l/F^ = t{B^) = 1.534 ±0.013 ps, the Am^ and Xd results are 
combined through Eq. (jHTj) to yield the world average 

Anid = 0.502 ± 0.006 ps"^ , (38) 
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Table 13: Time-dependent measurements included in the Am^ average. All these measurements 
have been adjusted to a common set of physics parameters before being combined. The CDF2 
and DO measurements are preliminary. The new BABAR and Belle [HH] measurements are 
not included here, but listed in Table IT^ 



Experiment 




Met no a 


Published value 


Adjusted value 


and Ref. 


rec. 


tag 


of A?TT,(^ in ps ^ 




of AttIj^ in ps~^ 




OO 


t 


Vjct 
1 


U.4U4±U.U45±U, 


,027 




AL/llvr rl 


[\)Z\ 


t 


n 4c;o-i-n ir^Q-i-n 
u.^ozinu.Udyinu. 


, U4^ 






1 CIO 1 

[92\ 


auove 


Lwo comuineu. 


n /loo-i-n (T^o-i-n 




n 441 -i-n (T^o-i-n noi 




no 


n* 


fi D 


U.4oz±U.U44±U. 


,024 


n /ico-i-n n/i/i-i-n no/i 

U.4oZlCU.U44lCU.UZ4 


DELPHI 


|,93j 


n 
I 


Vjct 


n /lOQ-un n/io i n 
U.4yo±U.U4Z±U. 


no7 
,UZ r 


0.504±0.04z±0.025 


DELPHI 




or* 


/n 

Vjct 

n 
I 


0.499±0.053±0. 


,015 


U.OUiitU.UOoitU.UiO 


DELPHI 




n 
I 


n A^D-i-n nAD-i-n 

U.^OUIlU.U^UIlU. 


,UOi 


f\ A c^'-r 1 /~i /~i /I /~i 4-0.049 

0.487±0.040 Zosm 


DELPHI 




D* 


D 

Vjct 


U.OZo±U.U z±U. 


n /I Q 
,U4o 


n p;i7-i-n n7o-i-n 04^? 


DELPHI 


^ 


vtx 


comb 


U.531±(j.(j25±(J. 


,007 


O.oo0±0.0zo±0.00o 


L3 




n 
I 




U.40oltU.U40ltU. 


,Uo/ 


0.4d9±0.04d±0.029 


L3 




g 


Vjet 


U.4z ±U.U44±U. 


C\AA 
,U44 


n /|Q«_Ln n/|/|_Ln n/IO 
U.4oOltU.U44ltU.U4Z 


L3 


ins 




I 


^(IPj 


U.4dz±U.UDo±U. 


,U0o 


n /|on_i_n ncQ_i_n n/i7 
U.4oU±U.UDo±U.U4 


OPAL 




D 


1] 
I 


U.4oU±U.U4o _Q 


,028 
,030 


n /ii^o_i_n n/iQ +0.018 
U.4DZ±U.U4o _o 017 


OPAL 


inzi 





Vjot 


n /i/i/i_i_n non +0 
U.444±U.Uzy _Q 


,020 
,017 


n /i^p;_i_n non +0.015 
U.400±U.Uzy _o.oi4 


OPAL 


inni 


L> I 


Vjet 


u.ooy±u.uou±u. 


no/1 
,U/4 


n F;/ip;_i_n n^n_i_n noo 
U.54o±U.UoU±U.Uzo 


OPAL 


inni 


L> 


/; 


U.ODttU.Uoy _Q 


,029 
,023 


n F^Ti _i_n non +0.028 
U.Ori±U.Uoy _o.o22 


OPAL 




TT I 


Vjet 


n /I n? 1 n no/i i n 
U.4y ( ±U.Uz4±U. 


,U/0 


n /in^_i_n no/i_i_n noc^ 
U.4yo±U.Uz4±U.UZO 


CDFl "TUU] 




QQT 


n /171 +0.078 +0, 
^ -0.068 -0, 


,033 
,034 


0/171 +0.078 +0.033 
^■^1^ -0.068 -0.034 


CDFl i 


lOlj 






U.0UoitU.U04±U. 


,U ( i 


n c;i/ij_n nfi/i-un n7n 
U.Oi4±U.Uo4±U.U ( U 


CDFl 


TM 


i 


^1 Qjet 


n c;nn-i-n nfi9-i-n 


, u^o 


n c^Qfi-un nc;o-i-n nQ7 


CDFl [Ml 


D*i 


0.516±0.099 


,029 
,035 


n coQ-Lfi nnn +0.029 
u.o,^oiLU.uyy _q q^c^ 


CDF2 [ 


TTH] 




SST,SMT,JQT 


0.536±0.037±0^ 


,017 


0.536±0.037±0;017 


CDF2 dnsi 




SST 


0.526±0.056±0, 


,005 


0.526±0.056±0.005 


D0 \im\ 




comb 


0.456±0.034±0, 


,025 


0.456±0.034±0.025 


BABAR [ 


107j 


B° 


I, K, NN 


0.516±0.016±0, 


,010 


0.518±0.016±0.008 


BABAR [ 


IHHl 


i 




0.493±0.012±0. 


,009 


0.489±0.012±0.007 


BABAR 


m 


D*iu 


/sT, NN 


0.492±0.017±0. 


,014 


0.490±0.017±0.013 


Belle 




comb 


0.528±0.017±0. 


,011 


0.529±0.017±0.011 


Belle mni 


D*£u 


comb 


0.494±0.012±0. 


,015 


0.496±0.012±0.014 


Belle mU 


L>*7r(part) i 


0.509±0.017±0. 


,020 


0.511±0.017±0.019 


Belle dl 


i 




0.503±0.008±0, 


,010 


0.504±0.008±0.009 


World averag 


;e (all above measurements included): 




0.502±0.004±0.005 


- ALEPH 


, DELPHI 


:, L3, OPAL and CDFl only: 




0.496±0.010±0.009 



- Above measurements of BABAR and Belle only: 0.503±0.005±0.005 
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Table 14: Simultaneous measurements of Am^ and t{B^), with the statistical correlations pstat 
between them. The Belle analysis also measures t{B^) at the same time, but it is converted 
here into a two-dimensional measurement of Am^ and t{B^), for an assumed value of t{B^). 
The first error on Am^ and t{B^) is statistical and the second one systematic; the latter includes 
a contribution obtained from the variation of t{B~^) or t{B^)/t{B^) in the indicated range. 
Units are ps~^ for Am,^ and ps for the lifetimes. The second BABAR result is still preliminary. 

Exp. & Ref. Measured Arrid Measured t{B^) pstat Assumed ri^B'^) 

BAR4R 0.492 ±0.018 ±0.013 1.523 ± 0.024 ± 0.022 -0.22 (1.083 ± 0.017)r(fi") 

BABAR ^ 0.523 ±0.004 ±0.007 1.501 ± 0.008 ± 0.030 -0.012 1.671 ± 0.018 

Belle ISHl 0.511 ±0.005 ±0.006 1.534 ± 0.008 ± 0.010 -0.27 1.635 ±0.011 



or, equivalently, 

Xrf = 0.770 ±0.011 and = 0.186 ± 0.003 . (39) 

Figure E] compares the Am^ values obtained by the different experiments. 

The B^ mixing averages given in Eqs. and and the 6-hadron fractions of Table El 
have been obtained in a fully consistent way, taking into account the fact that the fractions are 
computed using the Xd value of Eq. and that many individual measurements of Am^ at 
high energy depend on the assumed values for the 6-hadron fractions. Furthermore, this set of 
averages is consistent with the lifetime averages of Sec. 13.21 

It should be noted that the above average of Am^, Table ^1 and Figure |3] do not include 
two precise measurements that have been released for the Summer 2004 conferences by BABAR 
and Belle. The new BABAR analysis [SZI, based on partially reconstructed B° — s> D*£z/ decays, 
extracts simultaneously Am^ and t{B^) (similarly to the BABAR published result based on 
fully reconstructed B" D*iu decays jnSj while the new Belle analysis [SB] based on fully 
reconstructed hadronic B° decays and B° D*lv decays, extracts simultaneously Ama, t(B°) 
and r(B+). 

The results of the three B-factory analyses extracting Ama and t{B^) at the same time 
are listed in Table El Performing a two-dimensional combination of these results needs a 
careful assessment of the statistical and systematic correlations between Arrid and t{B^) in 
each analysis and between the measurements in the different analyses. This will be completed 
for a next version of this report. 

3.3.2 mixing parameters 

CP violation parameter \q/p\s 

No measurement or experimental limit exists on \q/p\si except in the form of a relatively 
weak constraint from CDF on a combination of \q/p\d and \q/p\s, f'dXdi)- — \q/'p\'d) + f'sXsi)- — 
\q/p\1) = 0.006 ± 0.017 [Hll, using inclusive semileptonic decays of b hadrons. The result is 
compatible with no CP violation in the mixing, an assumption made in all results described 
below. 

Mass difference Am^ 
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ALEPH 
(3 analyses) 

DELPHI * 
(5 analyses) 

L3 

(3 analyses) 

OPAL 
(5 analyses) 

CDFl * 
(4 analyses) 

CDF2* 
(2 prel. analyses) 

DO 

(1 prel. analysis) 

BABAR * 
(3 analyses) 

BELLE * 
(4 analyses) 

average of above 
after adjustments 

ARGUS+CLEO 
(Xj measurements) 

world average 



working group average 
without adjustments 

^0 




0.45 0.5 0.55 

Am . (ps"^) 



0.446±0.026±0.019 ps 
0.519±0.018±0.011 ps 
0.444±0.028±0.028 ps 
0.479±0.018±0.015 ps 
0.495±0.033±0.027 ps 
0.533±0.031±0.012 ps 
0.456±0.034±0.025 ps 
0.500±0.008±0.006 ps 
0.506±0.006±0.008 ps 

0.502±0.006 ps"^ 
0.493±0.032 ps'^ 

0.502±0.006 ps"^ 



Figure 4: The B^-B oscillation frequency Am^ as measured by the different experiments. The 
averages quoted for ALEPH, L3 and OPAL are taken from the original publications, while the 
ones for DELPHI, CDF, BABAR, and Belle have been computed from the individual results 
listed in Table IT^ without performing any adjustments. The time-integrated measurements of 
Xd from the symmetric B factory experiments ARGUS and CLEO have been converted to a 
Arrid value using t{B^) = L534 ± 0.013 ps. The two global averages have been obtained after 
adjustments of all the individual Ama results of Table IT^ (see text). 
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The time-integrated measurements of x (see Sec. 13.1.^ . when compared to our knowledge 
of Xd and the 6-hadron fractions, indicate that mixing is large, with a value of Xs close to 
its maximal possible value of 1/2. However, the time dependence of this mixing (called 
oscillations) has not been observed yet, mainly because the period of these oscillations turns 
out to be so small that it can't be resolved with the proper-time resolutions achieved so far. 

The statistical significance S of a. B^ oscillation signal can be approximated as |113j 



where and /sig are the number of B^ candidates and the fraction of B^ signal in the selected 
sample, w is the total mistag probability, and at is the resolution on proper time. As can be 
seen, the quantity S decreases very quickly as Am^ increases: this dependence is controlled 
by at, which is therefore the most critical parameter for Arris analyses. The method widely 
used for B^ oscillation searches consists of measuring a B^ oscillation amplitude A at several 
different test values of Arrig, using a maximum likelihood fit based on the functions of Eq. ()32|) 
where the cosine terms have been multiplied by A. One expects ^ = 1 at the true value of 
ArUs and to ^ = at a test value of Arris (far) below the true value. To a good approximation, 
the statistical uncertainty on A is Gaussian and equal to 1/iS |113j . 

Figures El and IHl show the amplitude spectra published by ALEPH |114j , CDF |115j , DEL- 
PHI inainniinnilllll, opal (nilllin] and SLD [nniinil-^ in each analysis, a particular 
value of Arris can be excluded at 95% CL if ^+ 1.645 a^ < 1, where cr^ is the total uncertainty 
on A. Because of the proper time resolution, the quantity cr^(Ams) is an increasing function of 
Arris (see Eq. which merely models l/a_A{Ams) since all results are limited by the available 
statistics). Therefore, if the true value of Arris were infinitely large, one expects to be able to 
exclude all values of Arris up to Am^^'^^, where Aml'^'^, called here the sensitivity of the analysis, 
is defined by 1.645 cr^(Am.^'^'^^) = 1. The most sensitive analyses appear to be the ones based 
on inclusive lepton samples at LEP, where reasonable statistics is available. Because of their 
better proper time resolution, the small data samples analyzed inclusively at SLD, as well as 
the few fully reconstructed B^ decays at LEP, turn out to be also very useful to explore the 
high ArUs region. 

These oscillation searches can easily be combined by averaging the measured amplitudes A 
at each test value of Arris- The combined amplitude spectra for the individual experiments are 
displayed in Fig. [71 and the world average spectrum is displayed in Fig.|Hl The individual results 
have been adjusted to common physics inputs, and all known correlations have been accounted 
for; in the case of the inclusive analyses, the sensitivities {i.e. the statistical uncertainties on 
A), which depend directly through Eq. (jiUj) on the assumed fraction /gig ~ fs of 5° mesons in 
an unbiased sample of weakly-decaying b hadrons, have also been rescaled to a common average 
of fs = 0.104 ± 0.015. The combined sensitivity for 95% CL exclusion of Arris values is found 
to be 18.2 ps~^. All values of Arris below 14.5 ps~^ are excluded at 95% CL, which we express 
as 



^ An unpublished analysis from SLD |122| , based on an inclusive reconstruction from a lepton and a topolog- 
ically reconstructed D meson, is not included in the plots or combined results quoted in this section. However, 
nothing is known to be wrong about this analysis, and including it would increase the combined Arris limit of 
Eq. l|lT|l by less than 0.1 ps^^ and the combined sensitivity by 0.9 ps^^. 




(40) 



Am, > 14.5 ps"^ at 95% CL . 



(41) 
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ALEPH lepton (13.1 ps 



a> 4 
%2 



DELPHI lepton ( 9.0 ps ') 




■2 




ALEPH D -lepton ( 7.5 ps 




DELPHI D -lepton + ^-lepton ( 8.6 ps ) 




OPAL leptons ( 7.2 ps ) 



OPAL D -lepton ( 4.2 ps ) 




Am fps ) 



Am fps ) 



Figure 5: i?°-oscillation amplitude spectra, displayed separately for each oscillation analysis. 
The points and error bars represent the measurements of the amplitude A and their total 
uncertainties cr^, adjusted to a set of physics parameters common to all analyses (including 
fs = 0.104±0.015). Values of Arris where the solid curve (^+ 1.645 cr^) is below 1 are excluded 
at 95% CL. The dashed curve shows 1.645 o"^; the number in parenthesis indicates where this 
curve is equal to 1, and is a measure of the sensitivity of the analysis, a) ALEPH inclusive 
lepton |114j . b) DELPHI inclusive lepton |117j . c) OPAL inclusive lepton and dilepton |118j . 
d) ALEPH Ds-e dHI, e) DELPHI Ds-E UHl and 0-£ [1111, f) OPAL D^-i [112] (continued on 
Fig.©. 
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ALEPH ( 0.4 ps ') 



CDFl ^-lepton ( 5.1 ps ) 




DELPHI B + D -hadron ( 3.2 ps ) DELPHI incl. vtx ( 6.9 ps ) 



4> t 








SLD D -tracks ( 3.3 ps ^) 



SLD incl. vtx ( 9.0 ps ^) 




15 20 

Am^ (ps"^) 



15 20 

Am^ (ps"^) 



Figure 6: (continuation of Fig. El) -B°-oscillation amplitude spectra, displayed separately for 
each oscillation analysis, in the same manner as in Fig. El a) ALEPH fully reconstructed 
5° [HI, b) DELPHI fully reconstructed E° and D.-hadron c) SLD +tracks [mj, d) 
CDF [nni, e) DELPHI inclusive vertex f) SLD inclusive vertex dipole [T^ . 
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ALEPH (14.4 ps ) 




SLD (11.5 ps ') 



d) 




DELPHI (11.9 ps ') 



CDFl ( 5.1 ps ') 




OPAL (7.9ps ^) 




all (18.2 ps"^) 



f) 




15 20 

Am^ (ps"^) 



10 15 20 
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Figure 7: Combined i?°-oscillation amplitude spectra, displayed separately for each experiment, 
in the same manner as in Fig. El a) ALEPH b) DELPHI [0311131^11111!, c) OPALfTTSl 

nT^ . d) SLD [Mlinni, e) CDF [nn!, f) aU experiments together. 
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' ' ' I ' ' ' ' I ' ' ' ' I 
- World average 




data±lo A 95% CL limit 14.5 ps"^ 
1.645 a -&■ sensitivity 18.2 ps" 

data ± 1.645 a 
data ± 1.645 a (stat only) 



10 12.5 15 17.5 20 22.5 25 



Am^ (ps"^) 



Figure 8: Combined measurements of the oscillation amplitude as a function of Arris, 
including all results published by Summer 2004 [113 El Eni ESI IM 1113 EZl EHl EHl 
I121j . The measurements are dominated by statistical uncertainties. Neighboring points are 
statistically correlated. 



35 



The values between 14.5 ps~^ and 21.7 ps~^ cannot be excluded, because the data is compatible 
with a signal in this region. However, no deviation from ^ = is seen in Fig. |H1 that would 
indicate the observation of a signal. 

It should be noted that most Am^ analyses assume no decay-width difference in the 
system. Due to the presence of the cosh terms in Eq. (|32p . a non-zero value of AF^ would 
reduce the oscillation amplitude with a small time-dependent factor that would be very difficult 
to distinguish from time resolution effects. 

Convoluting the average S° lifetime, 1.469 ± 0.059 ps, with the hmit of Eq. dH)) yields 

Xs = Am, r(5°) > 20.8 at 95% CL . (42) 

Under the assumption AF, = 0, i.e. Us = AFs/(2Fs) = (and no CP violation in the mixing), 
this is equivalent to 

2 2 

X, = > 0.49885 at 95% CL . (43) 

2{xi + 1) 

Decay width difference AF^ 

Information on AF^ can be obtained by studying the proper time distribution of untagged 
data samples enriched in mesons |123j . In the case of an inclusive B^ selection |^ or a 
semileptonic B^ decay selection |116U61j . both the short- and long-lived components are present, 
and the proper time distribution is a superposition of two exponentials with decay constants 
Ts ± AFs/2. In principle, this provides sensitivity to both Tg and (AFg/F^)^. Ignoring AF^ 
and fitting for a single exponential leads to an estimate of F^ with a relative bias proportional 
to (AF^/Fs)^. An alternative approach, which is directly sensitive to first order in AFg/F^, is 
to determine the lifetime of B^ candidates decaying to CP eigenstates; measurements exist for 
5° ^ J/^^ [271 Eg Ei and 5° ^ Di*'^^Di*^~ [1211, which are mostly CP-even states [T^ . 
However, a time-dependent angular analysis of B^ —>■ J / ipcf) allows the simultaneous extraction 
of AFs/F<i and the CP-even and CP-odd amplitudes |127j . An estimate of AFs/F^ has also 
been obtained directly from a measurement of the 5^ Di*'''^Di*^~ branching ratio |124j . 
under the assumption that these decays account for all the CP-even final states (however, no 
systematic uncertainty due to this assumption is given, so the average quoted below will not 
include this estimate). 

Present published data is not precise enough to efficiently constrain both F^ and AF^/F^; 
since the B^ and B^ lifetimes are predicted to be equal within a percent |1261 l24j. an expectation 
compatible with the current experimental data (see Table ITT|) . the constraint F^ = F^ can also 
be used to improve the extraction of AFg/F,,. Applying the combination procedure of Ref. [2j 
on the published results ED EH II211 EOl E21 yields 

|AF,|/F, < 0.54 at 95% CL (44) 

without external constraint, or 

|AF,|/F, < 0.29 at 95% CL (45) 

when constraining l/F^ to the measured B^ lifetime. This can be compared with the recent 
preliminary measurement of CDF |127j obtained from the time-dependent angular analysis of 
B^ — > J/ip(f) decays: 

AF,/F, = 0.65l[j:| ± 0.01 . (46) 
These results are not yet precise enough to test the Standard Model predictions. 
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4 Semileptonic B decays 



The original charge of the "Semileptonic B decays" working group consisted in 

• the determination of the best values of the inclusive and exclusive semileptonic B de- 
cay branching ratios from the combined data, both for Cabibbo-favored and Cabibbo- 
suppressed decays, fully taking into account correlations between the experiments; 

• determining the best possible values of F(l)|V"cfe| from exclusive semileptonic B decays 
and of |V"cfe| from the inclusive semileptonic B decay rate; 

• a detailed understanding of the correlated theoretical errors and the encouragement of a 
uniform and consistent error estimation among the active experiments. 

Currently, this program is implemented to varying degrees. The work of the "LEP Heavy 
Flavor Steering Group" has resulted in averages for Cabibbo-favored exclusive and inclusive 
semileptonic decays (see, e.g., Ref. |128p . Recently, the determination of \Vcb\ from inclusive 
decays has entered a new level precision with the inclusion of data on the moments of invari- 
ant mass and lepton energy distribution in semileptonic B decays. This approach is not yet 
implemented in HFAG. 

In this edition of the HFAG updates, we present an average of inclusive determinations 
of \Vub\. Work is in progress to achieve the same for exclusive decays and eventually for a 
combination of inclusive and exclusive results. See Ref. |129j for an alternative best estimate 
of iKfel- 

In the following a detailed description of all parameters and published analyses (including 
preliminary results) relevant for the determination of the combined results is provided. The 
description is based on the information available on the web-page at 

http : //www . slac . Stanford . edu/xorg/hf ag/ semi/ summer04/ summer04 . shtml 

In the combination of the published results, the central values and errors are rescaled to a 
common set of input parameters, summarized in Table IT^ and provided in the file common . par am 
(accessible from the web-page). All measurements with a dependency on any of these parame- 
ters are rescaled to the central values given in Table El and their error is recalculated based on 
the error provided in the column "Excursion" . The detailed dependency for each measurement 
is contained in files (provided by the experiments) accessible from the web-page. 

4.1 Exclusive Cabibbo-favored decays 

Aspects of the phenomenology of exclusive Cabibbo-favored B decays and their use in the 
determination of \Vcb\ in the context of Heavy Quark Effective Theory (HQET) are described 
in many places, e.g., in Ref. |128j and will not be repeated here. 

Averages are provided for both the branching ratios B{B^ — > D*~^i~V) plus B{B^ — * D^i~V) 
and the CKM matrix element \Vcb\ multiplied by the form factor at zero recoil of the decay 
B° D*^i~V and — ^ D^£~V, respectively. 
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Table 15: Common input parameters for the combination of semileptonic B decays. Most of 
the parameters are taken from Ref. [3]. This table is encoded in the file common, par am. The 
units are picoseconds for lifetimes and percentage for branching fractions. 



Parameter 


Assumed Value 


Excursion 


Description 


rb 


21.646 


±0.065 


Rb 


bdst 


1.27 


±0.021 


B{B D*tV) 


bdsd 


1.62 


±0.040 


B{B ^ D*D) 


bdst2 


0.65 


±0.013 


B{b D*t) (OPAL incl) 


bdsd2 
bdsdS 


4.2 
0.87 


±1.5 

+0.23 
-0.19 


Bib D*D) (OPAL incl) 
Bib D*D) (DELPHI incl) 


xe 


0.702 


±0.008 


B fragmentation: {Eb) / E^eam 


bdsi 


17.3 


±2.0 


B{b D*+ incl) 


cdsi 


22.6 


±1.4 


B{c D*+ incl) 


tbO 


1.534 


±0.013 




tbplus 


1.653 


±0.014 


t{B^) 


tbps 


1.442 


±0.066 


rm 


fbd 


39.8 


±1.0 


B^ fraction at i/i = m^o 


fbs 


10.5 


±1.5 


B^ fraction at = m^o 


fbar 


9.9 


±1.7 


Baryon fraction at ^/s = m^o 


dst 


67.7 


±0.5 


B{D*+ ^ 0^11+) 


dkpp 


9.20 


±0.6 


b\d^ A'-7r+7r+) 


dkp 


3.80 


±0.09 


b\d^ K~n+) 


dkpzp 


13.0 


±0.8 


i3(D0 ^ i^-7r+7r°) 


dkppp 


7.46 


±0.31 


bId° K-n+n+n~) 


dkzpp 


2.99 


±0.18 


b\d^ K%+Tr~) 


dkln 


7.0 


±0.4 


b\d^ K~(i+u) 


dkk 


4.3 


±0.2 


b\d^ K~K+) 


dkx 


1.100 


±0.025 


K^-K^X rates 


dkox 


0.42 


±0.05 


B{D^ K^X) 


dnlx 


6.87 


±0.28 


b\d^ xeu) 


dkpcl 


61.2 


±2.9 




dssR 


0.64 


±0.11 


Bib D**£V) X B{D** ^ D*+X) 


fbO 


49.4 


±0.8 


foo ^ siTiAS) B'^B'^) 


chid 


0.186 


±0.004 


Xdi time-integrated probability for B^ mixing 


chi 


0.092 


±0.002 


X = X (rVlOO) 



4.1.1 B ^D*+rv 

The measurements included in the average, shown in Table ITBl are scaled to a consistent set of 
input parameters and their errors. Therefore some of the (older) measurements are subject to 
considerable adjustments. 

• In order to reduce the dependence on theoretical error estimates, the central values and 
errors for the form factors i?i and i?2 are taken from the measurement by CLEO jL3(Jj . 
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However, all experiments (except for the CLEO |135j . the recent DELPHI and the 
BABAR |137j measurements) quote in their abstracts based on form factors (and 

their respective errors) from theory. Belle provides a second result evaluated with the 
CLEO form factors. All other experiments have recalculated F(l)|V^fe| to rely on the 
CLEO form factors. In the future, a substantial improvement in the error associated with 
form factors is expected by including form factor measurements at the B factories. 

• Updates in the branching fractions of D mesons and the production of D** mesons in B 
decays have generally lead to increased rates for B{B — > D*^l~V). 

• The average B^ lifetime has changed considerably since the ALEPH measurement, es- 
pecially with the much higher precision available at the B factories. This effect is less 
visible in the other measurements as they are more recent. 

• The production of B^ mesons at ^/s = rrizo has a direct impact on all LEP measurements. 
Adjusting results on the T(45') for the branching fraction of B{T{4S) B^B^) = 
yields an increase compared to the assumption of /''^ = 0.5. 

• Many input parameters are now known with a much increased precision — this decreases 
some of the systematic errors of the rescaled results with respect to the original publica- 
tion. 

The largest correlated errors are the fraction of mesons (fbd and fbO, respectively), the 
form factors -Ri(l) and -R2(l) at zero recoil, B^ meson lifetime, branching fractions of D mesons, 
and the details of D** modeling. 

At LEP, the measurements of 5^ ^ D*~^i~T' decays have been done both with "inclusive" 
analyses based on a partial reconstruction of the B D*^i^V decay and a full reconstruc- 
tion of the exclusive decay. The average branching ratio B{B — > D*'^i~V) is determined in 
a one- dimensional fit from the measurements provided in Table E| The statistical correla- 
tion between two analyses from the same experiment (DELPHI and OPAL, respectively) is 
taken into account. Figure IHl^a) illustrates the measurements and the resulting average. The 
X^/dof = 14.7/7 is slightly above the level where the PDG starts to introduce scale factors. 
The measurement by CLEO provides the largest contribution to the x^- 

The average for F(l)|V^b| is determined by the two-dimensional combination of the results 
provided in Table IT7I This allows to fully incorporate the correlation between F(l)|V"cfe| and p^. 
Figure ITUT b) illustrates the average F(l)|Vcfe| and the measurements included in the average. 
Figure imf a) provides a one-dimensional projection for illustrative purposes. The x^/dof = 
27.5/14 is below the level where the PDG starts to introduce scale factors. The measurement 
by CLEO provides the largest contribution to the x^- 

For a determination of \Vcb\, the form factor at zero recoil F{1) needs to be computed. A 
possible choice is F{1) = 0.91 ± 0.04thco |138j . resulting in 

\Vcb\ = (41.4 ± l.Oexp ± l.Sthco) X lO-l 

The value for F{1) and its error is based on a comparison of estimates using OPE sum rules 
and with an HQET based lattice gauge calculation (see Ref. jl38j for more details). 
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Table 16: Average branching ratio B{B D*^i v) and individual results. See the description 
in the text for explanations why the published results are lower than the rescaled results. 



Experiment 



D* 



V 



(rescaled) 



(published) 



ALEPH (excl) HHH 
OPAL (excl) 
OPAL (incl) 

DELPHI (incl) jH 

Belle (excl) ^EM 

cLEo (excl) unni 

DELPHI (excl) [THT)] 

BABAR (excl) [T!?7j 



syst 



5.78 ± 0.26stat ± 0.36, 
5.51±0.19stat±0.40syst 
6.20 ± 0.27stat ± 0.58syst 
5.04±0.13stat±0.36syst 
4.72 ± 0.23stat ± 0.42syst 
6.26 ± 0.19stat ± 0.39syst 
5.80 ± 0.22stat ± 0.47syst 
4.85 ± 0.07stat ± OMsyst 



5.53 ± 0.26stat ± 0.52syst 
5.11±0.19stat±0.49,yst 
5.92 ± 0.27stat ± 0.68,yst 
4.70 ± 0.13stat -0.31 syst 
4.60 ± 0.23stat ± 0.40syst 
6.09±0.19stat±0.40syst 
5.90 ± 0.22stat ± 0.50syst 
4.90 ± 0.07stat ± 0.36syst 



Average 



5.34 ± 0.20 



xVdof = 14.7/7 



Table 17: Average of -F(l)|Vcfe| determined in the decay B D*^i v and individual results. 
The fit for the average has x^/dof = 27.5/14. The total correlation between the average 
F{\)\Va\ and is 0.57. 



Experiment 


F(l)|Kb|[10-3] (rescaled) 
F(l)|Kfe|ilO^^] (published) 


(rescaled) 
(published) 


ALEPH (excl) ^ 


33.7 lb 2.1stat i 1.6syst 

31.9 ± 1.8stat ± 1.9syst 


0.75 ± 0.25stat ± 0.37,yst 
0.37±0.26stat±0.14,y,t 


OPAL (incl) 


38.6±1.2stat±2.4,y,t 
37.5 lb 1.2stat i 2.5syst 


1.25±0.14,tat±0.39,yst 
1.12±0.14stat±0.29syst 


OPAL (excl) [ng 


39.3 ± 1.6,tat ± l.Ssyst 
36.8 ± 1.6stat ± 2.0syst 


1.49 ± 0.21stat ± 0.26,yst 

1.31 ± 0.21stat ± 0.16syst 


DELPHI (incl) [HS] 


37.0 ± 1.4,tat ± 2.5syst 

35.5 lb 1.4stat _2.4syst 


1.50±0.14,tat±0.37,yst 

1.34 lb 0.14stat -0 22syst 


Belle (excl) [IMj 


36.5 ± 1.9stat ± 1.9syst 
35.8 ± 1.9stat ± 1.9syst 


1.45±0.16,tat±0.20,yst 
1.45±0.16,tat±0.20syst 


CLEO (excl) [133] 


43.7±1.3stat±1.8syst 

43.1±1.3stat±1.8syst 


1.61±0.09stat±0.21,yst 
1.61±0.09stat±0.21,yst 


DELPHI (excl) [HH] 


38.9 ± 1.8stat ± 2.1gyst 
39.2 ± 1.8stat =b 2.3syst 


1.32±0.15stat±0.33syst 
1.32 ± 0.15stat ± 0.33syst 


BABm (excl) [T37j 


35.4±0.3stat±1.6syst 
35.5 lb 0.3stat i 1.6syst 


1.29±0.03stat±0.27syst 
1.29±0.03stat±0.27syst 


Average 


37.7 ± 0.9 


1.55 ± 0.14 



4.1.2 B D+i-V 

The average branching ratio B{B^ — > D^i'V) is determined by the combination of the results 
provided in Table [TSl The error sources here are the same as discussed in Sec. I4.1.H but 
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ALEPH 

5.78 + 0.26 + 0.36 

OPAL (excl) 
5.51+019 + 0.40 

OPAL (partial reco) 
6.20 + 027 + 0.58 

DELPHI (partial reco) 
5.04 + O13+0.36 

BELLE 

4.72 + 0.23 + 0.42 
CLEO 

6.26 + 019 + 0.39 

DELPHI (cxcl) 
5.80 + 0.22 + 0.47 

BABAR 

4.85 ± 0.07 ± 0.34 

Average 
5.34 ± 0.20 



HFAG 



I 2004 I 
x'/Llor= l4.7/7(CL = 4.n^') 



a) 



B(B ^D*"rv)[%] 



ALEPH 

2.39 + 0.20 + 0.50 



CLEO 

2.10 + 0.13 + 0.17 



BELLE 

2.18+0.12 + O39 



Average 
2.13+0.20 



HFAG 



I 2004 I 

X^/dof = 0.3/2(CL = 86.6'?() 



1.5 



b) 



2.5 



B(B ^D"rv)[%] 



Figure 9: Average branching ratio of exclusive semileptonic B decays, (a) B D*~^i v and 
(b) ^ D^tV. 



ALEPH 

33.7+ 2.1 ± 1.6 * 

OPAL (partial reco) 

3S.6+ 1.2+ 2.4 

OPAL (excl) 

39.3+ 1.6+ 1.8 

DELPHI (partial reco) 

37.0+ 14+ 2.5 

BELLE 

36.5 ± 1.9 ± 1.9 

CLEO 

43.7+ 1.3+ 1.8 
DELPHI (excl) 
38.9+ 1.8+ 2.1 
BABAR 
35.4+ 0.3 + 1.6 
Average 
37.7+ 0.9 
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Figure 10: (a) Illustration of the average F(l)|V^b| and rescaled measurements of exclusive 
B D*'^l~T' decays determined in a two-dimensional fit. (b) Illustration of F(l)|Vcf,| vs. p^. 
The error ellipses correspond to A^^ = 1- 



generally at a higher level due to larger background levels, less stringent kinematic constraints, 
and larger kinematic suppression at the endpoint. Figure Efb) illustrates the measurements 
and the resulting average. 
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Table 18: Average of the branching ratio B{B D'^i u) and individual results. 



Experiment 


B{b'' D+e-V)[%] (rescaled) 


B(b" D+i-V)[%] (pubhshed) 


ALEPH [IM, 

CLEO dnn] 

Belle mo] 


2.39 ± 0.20stat ± O.SOsyst 

2.10±0.13stat±0.17syst 

2.18±0.12stat±0.39syst 


2.35 ± 0.20stat ± 0.44syst 
2.20±0.16stat±0.19syst 
2.13±0.12stat±0.39syst 


Average 


2.13 ± 0.20 


xVdof = 0.3/2 



The average for G(l)|V^b| is determined by the two-dimensional combination of the results 
provided in Table IT^ Figure HTJ^b) illustrates the average F(l)|Vcfe| and the measurements 
included in the average. Figure ITlT a) provides a one-dimensional projection for illustrative 
purposes. 

Table 19: Average of ^(l)!^;,! determined in the decay 'B^ D^l~V and individual results. 
The fit for the average has x^/dof = 0.3/4. 



Experiment 


G(l)|Kfe|[10-^] (rescaled) 
G'(l)|Kb|[10~^] (published) 


(rescaled) 
(published) 


ALEPH [131J 


39.9 ± lO.Ostat ± 6.4syst 
31.1 ± 9.9stat ± 8.6syst 


1.01 ± 0.98stat ± 0.37syst 
0.20 ± 0.98stat ± O.SOsyst 


CLEO [nni 


44.9 ± 5.8,tat ± 3.5syst 
44.8 ± e.lstat ± 3.7syst 


1.27±0.25stat±0.14syst 
1.30±0.27stat±0.14syst 


Belle [140] 


41.6±4.4,tat±5.2syst 

41.1±4.4,tat±5.1syst 


1.12±0.22stat±0.14syst 
1.12±0.22stat±0.14syst 


Average 


42.1 ± 3.7 


1.15 ± 0.16 



For a determination of \Vcb\i the form factor at zero recoil G{\) needs to be computed. A 
possible choice is G{1) = 1.04 ± 0.06theo |138j . resulting in 

iKfel = (40.4 ± 3.6exp ± 2.3thco) X 10"^ 
4.2 Inclusive Cabibbo-favored decays 

Aspects of the theory and phenomenology of inclusive Cabibbo-favored B decays and their use 
in the determination of |V^fe| in the context of the Heavy Quark Expansion (HQE), an Operator 
Product Expansion based on HQET, are described in many places (see, e.g., Ref. |141j and 
references therein). 

Averages are provided for the total semileptonic branching ratio B{B XtP) and the 
partial semileptonic branching fraction B{B — > XtP; > 0.6 GeV). 

The measurements of the total semileptonic branching ratio B{b XlV) at LEP (see, e.g., 
Ref IS] or Ref. |142j [) represent a different analysis class with a more explicit model dependence 



42 




Figure 11: (a) Illustration of the average G(l)|V^fe| and rescaled measurements of exclusive 
B — s> D^i~V decays determined in a two-dimensional fit. (b) Illustration of F(l)|Vc6| vs. p^. 
The error ellipses correspond to Ax^ = 1- 



than the (lepton-)tagged analyses used at the T{4S). Therefore the LEP measurements are not 
used in the averages computed here. 

4.2.1 Total semileptonic branching fraction 

The average for the total branching ratio B{B XtP) is determined by the combination of 
the results provided in Table QUI In this average, the extrapolation of the measured rate to the 
total decay rate is performed by each experiment, usually with a fit of several components to 
the experimental spectrum. 



Table 20: Average of the total semileptonic branching fractions B{B Xiu) determined in 
tagged measurements on the T{4S). 



Experiment 


Btot{B Xei7)[%] (rescaled) 


Btot{B Xev)[%] (published) 


ARGUS (£-tag) [US] 
BABAR (Breco-tag) [Hi] 
Belle (^tag) [TiH] 
BABAR (e-tag) [HH] 
Belle (B,eco-tag) [HI] 
CLEO (£-tag) [HE] 


9.76 ±0.50 ±0.39 
10.40 ±0.50 ±0.46 
10.97 ±0.12 ±0.49 
10.92 ±0.18 ±0.29 
11.19 ±0.20 ±0.31 
10.91 ±0.08 ±0.30 


9.70 ±0.50 ±0.60 
10.40 ±0.50 ±0.46 
10.90 ±0.12 ±0.49 

10.87 ±0.18 ±0.30 
11.19 ±0.20 ±0.31 

10.88 ±0.08 ±0.33 


Average 


10.90 ±0.23 


xVdof = 5.0/5 



It is possible to determine \Vcb\ from the total semileptonic branching fraction and the 
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lifetime of B mesons. However, this approach has limitations so that it is no longer applied. 
Recently, several groups have published global fits to the partial branching fraction plus the 
moments of the hadronic mass distribution and the lepton energy has been done by different 
groups [THTIITH^ . 

4.2.2 Partial semileptonic branching fraction 

The average for the visible branching ratio B{B XtP^pi > 0.6GeV/c) is determined by the 
combination of the results provided in Table For the determination of \ Vcb\, the extrapolation 
to the full spectrum is not necessary, as the HQE allows the direct determination of |V"cfe| from 
a partial rate. 



Table 21: Average of B{B X£u]pi > 0.6GeV/c) determined in (model-independent) lepton- 
tagged measurements on the T(4S'). 



Experiment 


B,.i,{B Xii7)[%] (rescaled) 


B^is{B XtU)[%] (pubhshed) 


Belle (£-tag) [145 
BABAR (e-tag) ^ 

CLEO (£-tag) [minn] 


10.31 ±0.11 ±0.47 
10.37 ±0.06 ±0.23 
10.23 ±0.08 ±0.22 


10.24 ±0.11 ±0.46 
10.36 ±0.06 ±0.23 
10.21 ±0.08 ±0.22 


Average 


10.29 ± 0.18 


X^/dof = 0.22/2 



4.2.3 Determination of \Vcb\ 

The determination of \Vcb\ directly from the visible branching fraction in combination with 
moments of the hadronic mass distribution and the lepton energy has been done by different 
groups |150t 11511 1152j . It is foreseen to use a similar approach here as well in the future. 

4.3 Exclusive Cabibbo-suppressed decays 

Here we list results on exclusive determinations of \ Vub\- An average of (these) exclusive h — s> utP 
results is envisioned for the future. The measurements are separated into two classes: a first 
one averaging over the entire range (shown in Table 0^ . and a second class, where the decay 
rate is measured differentially in (few) bins of (shown in Table . 

Table 22: Summary of exclusive determinations of B{B XtP) and \Vub\ using the entire 
q^ range. The errors quoted on \Vub\ correspond to statistical, experimental systematic and 
theoretical systematic, respectively. 



Experiment 


Mode 


B[10-4] 


\Vub\ [10-^] (rescaled) 


CLEO 
BABAR 
Belle USSl 


B° ^ p-^+i/ 
B° ^ p~l+v 


2.69 ±0.41 +^;^[;±0.50 
3.29 ± 0.42 ± 0.47 ± 0.60 
1.3 ±0.4 ±0.2 ±0.3 


3.24 ±0.25 +^:2i±0.58 
3.59 ± 0.23 ± 0.26 ± 0.66 
3.1 ±0.2 ±0.2 ±0.6 
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Table 23: Summary of exclusive determinations of B{B XtP) and \Vub\ binned in q^. The 
errors quoted on \ Vuh\ correspond to statistical, experimental systematic, theoretical systematic, 
and signal form-factor shape, respectively. 



Experiment 


Mode 




iKftI [10"^] (rescaled) 


CLEO jl56, 

CLEO dsni 

Belle dnZI 
Belle [T57j 
BABAR USHl 


B° p-~£+iy 

B° TT~i+U 

B^ p~i+u 
B^ 7r~e+u 


1.33 ± 0.18 ± 0.11 ± 0.01 ± 0.07 
2.17 ±0.34 ± 0.41 ± 0.01 
1.76 ±0.28 ±0.20 ±0.03 
2.54 ±0.78 ±0.85 ±0.30 
1.46 ±0.27 ±0.28 


2.88 ± 0.55 ± 0.30 tott ± 0.18 

Q Q/I 4_ n +0-27 +0.50 
0.0^ IC U.OZ _Q 3g „Q 4Q 

3.90 ± 0.71 ± 0.23 toil 



4.4 Inclusive Cabibbo-suppressed decays 

A recent discussion of the theoretical issues and errors is provided, e.g., in Refs. |159j and 
An independent estimate of the size of the theoretical uncertainties can be found in Ref . |129j . 
The following description is focusing on a technical description of the relevant details of the 
averaging procedure. 



4.4.1 Determination of \Vub\ 

Inclusive determinations of \Vub\ based on charmless semileptonic B decays have been presented 
by all four LEP experiments |161j and, more recently, by the experiments at the T(4S'). In all 
cases, the experimental approach is to measure the charmless semileptonic rate in a restricted 
region of phase space where the background from B XctP is suppressed. Theoretical input 
is used to extrapolate to the full rate and/or to determine \Vub\- The determination of \Vub\ 
from the full rate B{B — > XutP) is accomplished with the following formula: 



\Vub\ = 0.00424 ■ \j ^^^^Q^/""^ ^'^11 X (1.000 ± 0.028pert ± 0.039i/„,3), 

where tb is the average lifetime of B^ and i?"^ mesons and the theoretical error is based on 
Refs. |162l I163j using the moments measurements of the BABAR collaboration |151j . 

Here we combine results of the BABAR, Belle, and CLEO collaborations as they present 
measurements in well-defined phase-space regions and are based on the same theoretical de- 
scription of the decay in terms of an OPE. Due to the large background from Cabibbo-favored 
B — > XciT^ decays, the measurements are sensitive to B — ^ XuiV decays only in restricted 
regions of phase-space. The convergence of the OPE is affected by these restrictions and com- 
plicate the extrapolation to the full rate and thus the determination of \Vub\- 

In the theoretical description, the motion of the b quark inside the B meson is described with 
a "shape function" , usually in the two-parameter "exponential form" . These parameters can be 
determined from a measurement of the photon energy spectrum in the rare decay b —>■ s'j. This 
decay has a smaller branching fraction than the decay B — > XuiT' , which limits the precision 
possible in this approach. Recently, there has been theoretical interest to use moments measured 
in semileptonic B —>■ X^iV decays for the determination of the shape function parameters. 
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In the average \Vub\ presented here, all measurements are scaled to a common set of the 
shape function parameters and their uncertainties. Both are taken from the photon energy 
spectrum in i? — > 57 decays as measured by the Belle collaboration |164j . The (large) errors of 
the shape function parameters are dominated by statistical uncertainties in the measurement 
of the photon energy spectrum. 

For the average, the systematic error is divided into the categories shown in Table 1^ In the 
averaging procedure, the errors belonging to a specific category are taken as 100% correlated 
between the different measurements. The numerical values of the errors for each measurement 
are given in a spreadsheet |165j . 

Table 24: Errors in the average of \Vub\- The numerical values for all measurements can be 
found in a spreadsheet |165j . The entries in the "Name" column correspond to the names in 
the spreadsheet. 



Name 


Explanation 


Uncorr Quadratic sum of statistical and experimental errors 


Statistical 
Exp 


Statistical 

Experimental detector systematics (uncorrelated) 


Corr Quadratic sum of correlated errors 


B2C(bg) 
B2U(eff) 
BF^Vub 
fu(extrapol) 


Modeling of B ^ XctU (Branching fractions, form factors, etc) 
Modeling oi B ^ X^tu (Branching fractions. Ems-, ss popping, etc) 
Theoretical error F — > \ Vub\ {rrib, as) 
Extrapolation from visible range 


Additional errors not provided by the experiments 


SSF 
WA 
QHD 


Subleading shape functions (power corrections) 

Weak annihilation 

(Local) Quark-Hadron Duality 



In Table |211 the largest and most critical uncertainties are the extrapolation "fu" and the 
additional errors not provided by the experiments: 

• fu: This error estimates the uncertainty in the extrapolation from the observed phase- 
space to the full phase-space. For this, all measurements use the triple-differential cal- 
culation of Ref . . The shape function is parametrized with the "exponential form" . 
All measurements are rescaled to correspond to the shape function parameters deter- 
mined from the photon energy spectrum measurement of Belle |164j . In previous update 
(Winter 2004), the data of CLEO |167j had been used for this purpose. The change sub- 
stantially reduces the errors and leads to a shift in the central value. The extrapolation 
errors for each measurement are determined by the uncertainty of the shape function pa- 
rameters |168j . largely dominated by statistics. Recent theoretical work |169j may allow 
the extraction of the shape function parameters from the high-precision measurements of 
moments in semileptonic B X^tP decays. 

• SSF: Uncertainties due to power corrections are summarized into the term labeled "sub- 
leading shape function" effects. They have been calculated for the endpoint measure- 
ments |17()j : the errors of the corrections are used as estimate for the error (without 
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scaling the central values). These corrections strongly depend on the fraction of phase- 
space covered in a measurement. They are much larger for the CLEO and Belle endpoint 
analyses with the higher lepton energy threshold. We use the uncertainty for the case 
> 2.2GeV as estimate for the mx analysis of BABAR. This is roughly consistent with 
theoretical expectations |171j . For the {mx,q^) results, the error is expected to be sig- 
nificantly smaller; here we use half of the mx error. The SSF errors are summarized in 
Table 123 

• WA: The error for weak annihilation estimates the uncertainty from parametrically en- 
hanced nonperturbative effects predominantly expected at high-g^. The estimates in 
Ref. |172j are used for the (mx, analysis and scaled to the other phase-space regions 
as shown in Table 123 Because of the concentration at high-g^, the effects of weak anni- 
hilation are diluted as more phase-space is covered. 

• QHD: The error due to the violation of quark-hadron duality is estimated for the (mx, q^) 
analysis in Ref. |129j . It is found to be substantially smaller than the uncertainties from 
weak annihilation and power corrections. Therefore, we neglect it here. 



Table 25: Scaling of theoretical errors in different regions of phase-space. The mx analysis 
requires for the hadronic system in 5 — >• Xutu decays mx < 1.55 GeV/c^, the (mx, q^) analysis 
selects the phase-space nix < 1-7 GeV/c^, > 8 GeV^. 



Component 


mx 


{mx,q^) 


El > 2.0GeV 


Ee > 2.2GeV 


Ei > 2.3GeV 


SSF 


4% 


2% 


2% 


4% 


6% 


WA 


2% 


4% 


4% 


6% 


8% 



The average of \Vub\ is based on the measurements summarized in Table |2E1 The BABAR 
and Belle collaborations provided these central values. The average B meson lifetime used 
varies by less than 0.5% between the BABAR and Belle collaboration; no rescaling to a common 
lifetime value is applied because of its small effect. No re-evaluation of the endpoint results 
was available from the CLEO collaboration. Therefore their previous result was scaled with 
(2.23/1.72)^/^, the ratio of branching fractions for > 2GeV as determined with the Belle and 
CLEO b ^ s'j photon energy spectra, as tabulated in Ref. |17(ij . The relative extrapolation 
error in Ref. |176j is applied for the CLEO measurement. After this adjustment, the average 
\Vub\ amounts to 

= (4.70 ±0.44) X lO-\ 

This error contains contributions from the perturbative expansion and the uncertainty of 
the b quark mass. The error of this average contains a relative contribution of 4.8% for the 
translation of F(6 — > uiv) to \Vub\, which has been derived using the semileptonic moments 
measurements of the BABAR collaboration |151j . 

The current extrapolation error on \ Vub\ is based on the approach of Ref. |166j for the shape 
function effect, where the shape function parameter errors are based only on the photon energy 
spectrum. Recent reports I18()j studied the radiative corrections to the shape function and 



47 



the prescription of Ref. must be modified. In particular Bosch et al. |169j showed that 
the sensitivity of the extrapolation to the shape function is reduced for the mx and {mx^(f) 
analyses, but becomes larger for the endpoint analyses. Furthermore, it may become possible 
to substantially improve the precision of the shape function parameters (and hence reduce the 
extrapolation errors) when using the moments measurements in semileptonic B decays. 

Table 26: Average of \Vuh\ determined in inclusive measurements on the T(4S'). The errors 
quoted for the measurements correspond to uncorrelated and correlated, respectively. The 
correlated error of the measurements do not include the common and constant error of 2.8% © 
3.9% from the full rate to \Vub\'i however, this error is included in the average \Vub\- 



Experiment 


Method ([GeV^^^]) 


\Vub\ [10~3] (rescaled) 


\Vub\ [10^3] (published) 


CLEO [nni 

Belle [mi 
Belle [T7H] 
BABAR [T7H] 
BABAR [TTTj 
BABAR [TTTj 
BABAR [UHl 
Belle [T7U1 


Endpoint, pn > 2.2 

{mx < 1.70, > 8) 

Endpoint, pe > 2.3 

Endpoint, pe > 2.0 

nix < 1-55 

{mx < 1.70, > 8) 

{Ee,q') 

{mx < 1.70, q^ > 8) 


4.69 ±0.23 ±0.63 
4.75 ± 0.46 ± 0.46 
4.46 ±0.23 ±0.61 
4.40 ±0.15 ±0.44 
5.22 ±0.30 ±0.43 
5.18 ±0.52 ±0.42 
4.99 ±0.34 ±0.51 
5.54 ±0.65 ±0.54 


4.08 ±0.20 ±0.62 
4.66 ±0.45 ±0.61 
3.99 ±0.20 ±0.61 
4.40 ±0.15 ±0.44 
5.22 ±0.30 ±0.43 
5.18 ±0.52 ±0.42 
4.99 ±0.34 ±0.51 
5.54 ±0.65 ±0.54 


Average 




4.70 ± 0.44 


xVdof = 6.7/7 



ALEPH 

4.12±0.67±0.71 
L3 

5.70 + 1.00 ± 1.40 

DELPHI 

4.07 ±0.65 ±0.61 

OPAL 

4.00 ±0.71 ±0.71 



CLEO (endpoint) 

4.69 ± 0.23 ± 0.63 
BELLE sim. ann. (m^, Q") 
4.75 ± 0.46 ± 0.46 
BELLE (endpoint) 

4.46 ±0.23 ±0.61 
BABAR (endpoint) 
4.40 ±0.15 ±0.44 
BABAR mj, 
5.22 ± 0.30 ± 0.43, 
BABAR (mj;, Ql 
5.18 ± 0.52 ±0.4^ 
BABAR (E,, Q ) 
4.99 ±0.34 ±0.51 
BELLE B„„(m^. Q"") 
5.54 ± 0.65 ± 0.54' 
Average 

4.70 ± 0.44 



a) 



HFAG 



J 



X-/dof=6.7/7(CL = 46 5%) 



IV,J [X 10-1 



CLEO(B"^p-l-*v) 

3.24±0J4 ±0.5K 1 1 U 1 1 

BABAR (B -* p"l%) , , _ 

3.59 + 0.35 + 0.66 1 1 U 1 




b) 


3..O±0.3(j,±O.6O 1 1 A 1 1 
CLEO(B -»it-l V) , ■ ^ , ■ 
3.24 + 0.26 + 0.4S 1 1 O 1 1 




CLEO(B -»p"l*v) , , r-, . , 

3.00 ±QJ9± 032 1 1 U 1 1 

BellelB ^nlvl , , ^ 






3.90 ±0.70 ±0.50 1 1 O 






ALEPH 






4,12 + 0.67 ±0.71 1 1 O- 
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Figure 12: (a) Illustration of the average \Vub\ from inclusive semileptonic B decays measured 
at the T{4:S). The measurements from LEP |161j are shown only for illustration, they are not 
included in the average, (b) Illustration of additional exclusive \Vub\ measurements together 
with the average \Vub\ from inclusive semileptonic B decays measured at the T(4S'). 
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4.4.2 Discussion 

In the future, improvements to this average will be studied and implemented. The following is 
a starting point, not an exhaustive list, for possible improvements. 

• The paper by de Fazio and Neubert with a shape function parametrization provides 
a triple differential decay rate for B — *• X^tP decays that is used in all measurements 
over the entire phase-space. Recent theoretical work seems to indicate that it can be 
improved even at its nominal order (l/m^ and Ug) |169| I18U] . 

• Every experiment should provide the result in an unfolded rate so that theoretical im- 
provements can be applied in the future. The idea is to separate the "unfoldings" of 
a rate or distribution from the "extrapolation" (the former is experimental, the latter 
theoretical). 

• All three experiments have a "hybrid" MC model for signal decays. The differential 
spectra should be compared. Its applicability should be understood (probably good for 
unfolding, but questionable for extrapolation). 

• There should be agreement on one specific recipe to evaluate the common errors. This 
can include items like 

— B ^ X(.tP modeling: 

- Variation of inclusive and exclusive branching fractions for B and D decays (full 
list and range) 

- Variation of form factors describing B XciT' decays 

— B — s> XuiV modeling and theoretical error 

- Variation of exclusive branching fraction (full list and range) 

- Comparison with inclusive-only signal model and comparison with exclusive- 
only signal model. 

- ss popping (if applicable in the analysis) 

- Central values of shape function parameters 

- Variation of shape function parameters (this is linked to recent theoretical work) 

A close collaboration between the experiments will be necessary to obtain an even more con- 
sistent treatment of errors not only for the most critical ones (as shown here), but also for the 
remaining smaller errors. 

5 Measurements related to Unitarity Triangle angles 

The charge of the "CP(t) and Unitarity Triangle angles" group is to provide averages of mea- 
surements related (mostly) to the angles of the Unitarity Triangle (UT). To date, most of the 
measurements that can be used to obtain model-independent information on the UT angles 
come from time-dependent CP asymmetry analyses. In cases where considerable theoretical 
input is required to extract the fundamental quantities, no attempt is made to do so at this 
stage. However, straightforward interpretations of the averages are given, where possible. 
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In Sec. 15.11 a brief introduction to the relevant phenomenology is given. In Sec. 15.21 an 
attempt is made to clarify the various different notations in use. In Sec. 15.31 the common 
inputs to which experimental results are rescaled in the averaging procedure are listed. We 
also briefly introduce the treatment of experimental errors. In the remainder of this section, 
the experimental results and their averages are given, divided into subsections based on the 
underlying quark-level decays. 

5.1 Introduction 

The Standard Model Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix V must be 
unitary. A 3 x 3 unitary matrix has four free parameters,^ and these are conventionally written 
by the product of three (complex) rotation matrices |181j . where the rotations are characterized 
by the Euler angles 6*12, ^13 and 6*23, which are the mixing angles between the generations, and 
one overall phase 6, 

(Vud Vus Vub \ I C12C13 S12C13 5136^*"^ \ 

Vcd Vcs Vcb = -S12C23 - Ci2S23Si3e*'' C12C23 - Si2S23Sl3e*'' S23C13 (47) 
Vtd Vts Vtb ) \ S12S23 - Ci2C23Si3e*'^ -C12S23 - Si2C23Si3e*'^ C23C13 / 

where Cij = cos 9ij, Sij = sin 9ij for i < j = 1,2,3. 

Following the observation of a hierarchy between the different matrix elements, the Wolfen- 
stein parameterization |182j is an expansion of V in terms of the four real parameters A (the 
expansion parameter). A, p and 77. Defining to all orders in A |183j 

S12 = A, 

S23 = AX^ (48) 
5136"*'^ = AX%p-tr]), 

and inserting these into the representation of Eq. (|47p. unitarity of the CKM matrix is achieved 
to all orders. A Taylor expansion of V leads to the familiar approximation 

/ 1 - AV2 A AX%p-ir]) \ 

V = [ -X 1 - AV2 AX^ ] + (A^) . (49) 

\ AX^{1- p~iri) -AX"^ 1 / 

The non-zero imaginary part of the CKM matrix, which is the origin of CP violation in the 
Standard Model, is encapsulated in a non-zero value of 77. 

The unitarity relation VW = 1 results in nine expressions which can be written Ylii=u c t ^ij^ik = 
Sjk, where Sjk is the Kronecker symbol. Of the off-diagonal expressions (j 7^ k), three can be 
trivially transformed into the other three (under j ^ k), leaving six relations, in which three 
complex numbers sum to zero, which therefore can be expressed as triangles in the complex 
plane. 

One of these, 

VudV:, + v^v;, + VtdV;, = o, (50) 

is specifically related to B decays. The three terms in Eq. (j5(jp are of the same order [O (A'^)), 
and this relation is commonly known as the Unitarity Triangle. For presentational purposes, it 
is convenient to rescale the triangle by iVcdVX)'^, as shown in Fig. [T51 

^In the general case there are nine free parameters, but five of these are absorbed into unobservable quark 
phases. 
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Two popular naming conventions for the UT angles exist in the literature: 



arg 



P = (pi = arg 



7 = 03 = arg 



In this document the 7) set is used. 

The apex of the Unitarity Triangle is given by the following definition |183j 



p + irj 



(51) 



The sides Ru and Rt of the Unitarity Triangle (the third side being normalized to unity) are 
given by 



Ru 
Rt 



'^udVub 



(52) 
(53) 



5.2 Notations 

Several different notations for CP violation parameters are commonly used. This section reviews 
those found in the experimental literature, in the hope of reducing the potential for confusion, 
and to define the frame that is used for the averages. 

In some cases, when B mesons decay into multibody final states via broad resonances (p, 
K*, etc.), the experiments ignore interference effects in the analyses. This is referred to as the 
quasi-two-body (Q2B) approximation in the following. 



5.2.1 CP asymmetries 

The CP asymmetry is defined as the difference between the rate involving a h quark and that 
involving a h quark, divided by the sum. For example, the partial rate (or charge) asymmetry 
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for a charged B decay would be given as 



Af = 



TiB-^f)-TiB+^f) 



T{B~^f) + T{B+^f) 
5.2.2 Time-dependent CP asymmetries in decays to CP eigenstates 



(54) 



If the amphtudes for 5° and B^ to decay to a final state /, which is a CP eigenstate with 
eigenvalue rjf, are given hj Af and Af, respectively, then the decay distributions for neutral B 
mesons, with known flavour at time At = 0, are given by 



TBO^fiAt) 



-|Ai|/r(BO) 

4r(50) 

-|At|/r(BO) 

4r(E0) 



2Im(A/) 
2 Im(A/ 



1 + |A 



/I 



2 sin(AmAt) — 
sin(AmAt) + 



1 - 




cos(AmAt) 


1 + 




1 - 


lA/P 


cos(AmAt) 


1 + 


|A/P 



(55) 
(56) 



Here A/ = -4^ contains terms related to B^-B^ mixing and to the decay amplitude (the 



, - 1^ 

eigenstates of the effective Hamiltonian in the system are \B±) = p ± q |-B°)). This 

formulation assumes CPT invariance, and neglects possible lifetime differences in the neutral 
B meson system. The time-dependent CP asymmetry, again defined as the difference between 
the rate involving a b quark and that involving a b quark, is then given by 



Af (At) 



Tgo^fiAt) -TBO^fiAt) _ 2Im(A/) 



^BO^fiAt) 



sin(AmAt) 



1 



I \ |2 



1 + |A, 



cos(AmAt). (57) 



BO^fiAt) l + \Xf\ 

While the coefficient of the sin(AmAt) term in Eq. (|37j) is everywhere^ denoted Sfi 

2lm{Xf) 



S 



f 



|2 ' 



1 + IA/i 

different notations are in use for the coefficient of the cos (Am At) term: 



(5J 



C 



f 



-A 



lA 



/ 



/I 



1 + IA; 



|2 ■ 



(59) 



The C notation is used by the BABAR collaboration (see e.g. |192j ). and also in this document. 
The A notation is used by the Belle collaboration (see e.g. |58j). 

Neglecting effects due to CP violation in mixing [i.e., taking \q/p\ = 1), if the decay ampli- 
tude contains terms with a single weak phase then |Aj| = 1 and one finds Sf = — ?7/ sin(0inix + 



Cf = 0, where 



aig^q/p) and 0dec — ^-rg (Af/Af). Note that 



2/3 in the 



Standard Model (in the usual phase convention). If amplitudes with different weak phases 
contribute to the decay, no clean interpretation of Sf is possible. If the decay amplitudes have 
in addition different CP conserving strong phases, then |A/| ^ 1 and no clean interpretation is 
possible. The coefficient of the cosine term becomes non-zero, indicating direct CP violation. 
The sign of A/ as defined above is consistent with that of in Eq. (jS3]). 

■ sin(AmAi) + Aj:" cos(AmAi), or similar. 



^Occasionally one also finds Eq. written as Af {At) 
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5.2.3 Time-dependent CP asymmetries in decays to vector-vector final states 

Consider B decays to states consisting of two vector particles, such as J/ipK*^{-^ /CgTr"), 
D*+]j*~ and p~^p~ , which are eigenstates of charge conjugation but not of parity/ In fact, for 
such a system, there are three possible final states; in the helicity basis these can be written 
ho, /i+i. The ho state is an eigenstate of parity, and hence of CP; however, CP transforms 
/i+i ^ h^i (up to an unobservable phase). In the transversity basis, these states are transformed 
into = (h+i + h^i) /2 and h± = (/i+i — /;,_i)/2. In this basis all three states are CP eigenstates, 
and h± has the opposite CP to the others. 

The amplitudes to these states are usually given by ^o,±,|| (here we use a normalization 
such that l^oP + + = !)• Then the effective CP of the vector-vector state is known 
if |^_Lp is measured. An alternative strategy is to measure just the longitudinally polarized 
component, |^oP (sometimes denoted by /long); which allows a limit to be set on the effective 
CP since < + l^|| P = 1 ~ l^oP- The most complete treatment for neutral B decays 

to vector-vector final states is time-dependent angular analysis (also known as time-dependent 
transversity analysis). In such an analysis, the interference between the CP even and CP odd 
states provides additional sensitivity to the weak and strong phases involved. 



5.2.4 Time-dependent CP asymmetries in decays to non-CP eigenstates 

Consider a non-CP eigenstate /, and its conjugate /. For neutral B decays to these final 
states, there are four amplitudes to consider: those for B^ to decay to / and f {Af and Aj, 
respectively), and the equivalents for B^ [Aj and Aj). If CP is conserved in the decay, then 
Af = Aj and Aj = Aj. 

The time-dependent decay distributions can be written in many different ways. Here, we 

follow Sec. 15.2.21 and define A/ = and Aj = The time-dependent CP asymmetries 

then follow Eq. (lITTjl : 

AfiAt) ^ ^ ^;J,^/ = SfSmiAmAt)-CfCos{AmAt), (60) 



AjiAt) = / ,A,Lr ^A, = SjMAmAt)-CjcosiAmAt), (61) 



with the definitions of the parameters Cj, Sf, Cj and Sj, following Eqs. and (j^ . 
The time-dependent decay rates are given by 

;i + {Afj)) {1 + 5/ sin(AmAt) - Cj cos(AmAt)} , (62) 

;i + (Afj)) {1 - Sf sin(AmAt) + C/ cos(AmAt)} , (63) 

;i - {Afj)) {l + Sjsm{AmAt) - Cj cos(AmAt)} , (64) 

;i - {Afj)) {1 - Sjsm{AmAt) + Cjcos(AmAt)} , (65) 







g-|At|/r(BO) 




./(At) = 


8r(S0) 






g-|At|/r(BO) 




./(At) = 


8r(S0) 






^-\At\/T{B°) 




./(At) = 


8r(50) 






g-|At|/r(i?0) 




./(At) = 


8r(50) 



''This is not true of all vector-vector final states, e.g., D*^p^ is clearly not an eigenstate of charge conjugation. 
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where the time-independent parameter {A^j) represents an overall asymmetry in the production 
of the / and / final states,* 

(^/7) = . ,^,2 ,/,. 'I - (66) 



AfV + \M ) + [\Aj\ +\Aj\ 

Assuming \q/p\ = 1, the parameters Cf and Cj can also be written in terms of the decay 
amplitudes as follows: 

g/= /' ' |2 and /' (67) 

giving asymmetries in the decay amplitudes of B'^ and 5° to the final states / and / respectively. 
In this notation, the direct CP invariance conditions are {^fj) = and Cf = —Cj. Note 
that Cf and Cj are typically non-zero; e.g., for a fiavour-specific final state, Af = Aj = 
{Af = Aj = 0), they take the values Cf = -Cj = 1 {Cf = ~Cj = -1). 

The coefficients of the sine terms contain information about the weak phase. In the case 
that each decay amplitude contains only a single weak phase {i.e., no direct CP violation), 
these terms can be written 

-2 sin(0mix + 0dcc - 5/) -2|Aj| |Zj|sin(0mix + 0dec + 5/) 

\Af\'+\Af\ \Aj\ + \Aj\ 

(68) 

where Sf is the strong phase difference between the decay amplitudes. If there is no CP violation, 
the condition Sf = —Sj holds. If amplitudes with different weak and strong phases contribute, 
no clean interpretation of Sf and Sj is possible. 

Since two of the CP invariance conditions are Cf = —Cj and Sf = —Sj, there is motivation 
for a rotation of the parameters: 

S f -\- S~f S f — S~f C f + cv C f — C~f 

%=^^' ^%=^T^' % = ^^' ^%=^^- (69) 

With these parameters, the CP invariance conditions become Sfj = and Cfj = 0. The 
parameter AC^y gives a measure of the "fiavour-specificity" of the decay: AC^y = ±1 corre- 
sponds to a completely fiavour-specific decay, in which no interference between decays with and 
without mixing can occur, while ACjj = results in maximum sensitivity to mixing-induced 
CP violation. The parameter AS'^-j is related to the strong phase difference between the decay 
amphtudes of to / and to /. We note that the observables of Eq. (jU^ exhibit experi- 
mental correlations (typically of ~ 20%, depending on the tagging purity, and other effects) 
between Sfj and A5'jj, and between Cfj and AC^y. On the other hand, the final state specific 
observables of Eq. ()6()|1 tend to have low correlations. 

Alternatively, if we recall that the CP invariance conditions at the amplitude level are 
Af = Aj and Aj = Af, we are led to consider the parameters |216j 

l^|2 i<i2 l^|2 |.|2 

Al= J 2 . and A-jf= \J\, /' . (70) 



^This parameter is often denoted Af (or Acp), but here we avoid this notation to prevent confusion with 
the time-dependent CP asymmetry. 
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These are sometimes considered more physically intuitive parameters since they characterize 
direct CP violation in decays with particular topologies. For example, in the case of p^-K^^ 
(choosing / = p+vr" and / = p~7i~^), Afj (also denoted A'^~) parameterizes direct CP violation 
in decays in which the produced p meson does not contain the spectator quark, while Ajj (also 
denoted A'p^) parameterizes direct CP violation in decays in which it does. Note that we 
have again followed the sign convention that the asymmetry is the difference between the rate 
involving a h quark and that involving a h quark, cf. Eq. ()54|) . Of course, these parameters are 
not independent of the other sets of parameters given above, and can be written 

{A,j) + C,j + {A,j)AC,j _ -{A,j) + C,j+{A„)AC„ 

l + ACa+{A„)C„ -1+AC„+(A„)C„ ^ 

They usually exhibit strong correlations. 

We now consider the various notations which have been used in experimental studies of 
time-dependent CP asymmetries in decays to non-CP eigenstates. 

The above set of parameters {{A^j), Cj, Sj, Cj, Sj), has been used by both BABAR j212j 
and Belle j2lHj in the D*^D^ system (/ = D*+D~, J = D*~D+). However, slightly different 
names for the parameters are used: BABAR uses {A, 5'+-, C^+, 5'„+); Belle uses (.4, C+, 
5+, C_, S^). In this document, we follow the notation used by BABAR. 

In the p^TT^ system, the {{Ajj), Cj-j, Sj^j, AC^y, AS'jj) set of parameters has been used 
originally by BABAR |221j . and more recently by Belle |223j . in the Q2B approximation; the 
exact names^ used in this case are {A^^p, Cpj,, SpT,, ACp^, AS'p^), and these names are also used 
in this document. 

Since p'^ir^ is reconstructed in the final state vr+vr^Tr'^, the interference between the p res- 
onances can provide additional information about the phases. BABAR |222j has performed a 
time-dependent Dalitz plot analysis, from which the weak phase a is directly extracted. In such 
an analysis, the measured Q2B parameters are also naturally corrected for interference effects. 

Time-dependent CP analyses have also been performed for the final states D'^-k^ , D*^7i^ 
and D^p"^. In these theoretically clean cases, no penguin contributions are possible, so there 
is no direct CP violation. Furthermore, due to the smallness of the ratio of the magnitudes 
of the suppressed (6 — u) and favoured (6 c) amplitudes (denoted Rf), to a very good 
approximation, Cf = —Cj = 1 (using / = D^*'>~h~^, f = D^*^~^h~ h = 7i,p), and the coefficients 
of the sine terms are given by 

Sf = -2Rf sin(0mix + 0dec - Sf) and Sj = -2Rf sin(0mix + 0dec + ^/)- (72) 

Thus weak phase information can be cleanly obtained from measurements of Sf and Sj, al- 
though external information on at least one of Rf or Sf is necessary. (Note that 0mix -l- (pdec = 
2/3 + 7 for the decay modes in question, while Rf and Sf depend on the decay mode.) 

^BaBAR has used the notations A'^p \221\ and Ap-rr \222\ in place of A^'p. 
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Table 27: Conversion between the various notations used for CP violation parameters in the 
D^TT^, Z}*^7r^ and D^p^ systems. The hi terms used by BABAR have been neglected. 







BABM, 


Relle nartial rer 


Relle full rer 


Sd+tt- 




= -(a + q) 


N/A 


2RDn sin (201 + 03 + Son) 


Sd--k+ 


-s+ 


= -(a - Ci) 


N/A 


2Rd^ sin(20i + 03 - 5dt,) 


Sd*+-k- 




= -(a + Ci) 


5+ 


-2Rd*^ sin(20i + 03 + (J^,.^) 


SD*--rr+ 


-s+ 


= -(a - Ci) 




-2RD*n Sin(20i + 03 - Sn^n) 


Sd+p~ 




= -(a + Ci) 


N/A 


N/A 


Sd-p+ 




= -(a - q) 


N/A 


N/A 



Again, different notations have been used in the literature. BABAR \229\ I231j defines the 
time-dependent probability function by 

e-|At|/r 

f^{r], At) = [1 T S^; sin(AmAt) t cos(AmAt)] , (73) 

where the upper (lower) sign corresponds to the tagging meson being a (B^)- [Note here 
that a tagging B^ (B^) corresponds to —5*^ (+5*^).] The parameters r] and ( take the values 
+1 and + (—1 and — ) when the final state is, e.g., D~7r~^ (D+vr"). However, in the fit, the 
substitutions = 1 and 5^ = a^rihi—rjCi are made. [Note that, neglecting h terms, 5*+ = a — c 
and S'_ = a + c, so that a = {S+ + S^)/2, c = {S^ — S+)/2, in analogy to the parameters of 
Eq. (jHUj) .] The subscript i denotes the tagging category. These are motivated by the possibility 
of CP violation on the tag side |233j , which is absent for semileptonic B decays (mostly lepton 
tags). The parameter a is not affected by tag side CP violation. 

The parameters used by Belle in the analysis using partially reconstructed B decays |232j . 
are similar to the S*^ parameters defined above. However, in the Belle convention, a tagging B^ 
corresponds to a + sign in front of the sine coefficient; furthermore the correspondence between 
the super /subscript and the final state is opposite, so that S± (BABAR) = —S"^ (Belle). In 
this analysis, only lepton tags are used, so there is no effect from tag side CP violation. In the 
Belle analysis using fully reconstructed B decays |23Uj . this effect is measured and taken into 
account using D*/z/ decays; in neither Belle analysis are the a, b and c parameters used. In the 
latter case, the measured parameters are 2i?£,(*)^ sin(20i + 03 ± 6d(*).^); the definition is such 
that 5"^ (Belle) = — 2i?£).^ sin(20i + 03 ± 6D*n)- However, the definition includes an angular 
momentum factor (—1)^ |234j . and so for the results in the Dn system, there is an additional 
factor of —1 in the conversion. 

Explicitly, the conversion then reads as given in Table EZl where we have neglected the bi 
terms used by BABAR. For the averages in this document, we use the a and c parameters, and 
give the explicit translations used in Table |2H1 It is to be fervently hoped that the experiments 
will converge on a common notation in future. 

Time- dependent asymmetries in radiative B decays 

As a special case of decays to non-CP eigenstates, let us consider radiative B decays. Here, 
the emitted photon has a distinct helicity, which is in principle observable, but in practice is 



56 



Table 28: Translations used to convert the parameters measured by Belle to the parameters 
used for averaging in this document. The angular momentum factor L is —1 for D*tt and +1 
for Dt[. 



D*TX partial rec. 








full 


rec. 




a -{S+ + S-) 

c -(V - S-) 


i(-l)^+M2i?^(,), 


sin(20i + 1 
sin(20i + ( 


h - 
h - 




) + 2i?^w^ sin(20i + ( 
) - 2i?BM^ sin(20i + , 





not usually measured. Thus the measured time-dependent decay rates are given by |184[ I185j 

Tso^x.m = Tso^^^^m + ^BO^x.^m (74) 

-\At\/riBO) 

= ^^^^0) {1 + + S^) sin(AmAt) - {Cl + Cr) cos(AmAt)} , 
Tso^x^iAt) = TBO^x^^{At) + TBO-.x^^{At) (75) 

-\At\/T(BO) 

= g rno^ + Sr) sin(AmAt) + {Cl + Cr) cos(AmAt)} , 



where in place of the subscripts / and / we have used L and R to indicate the photon helicity. 
In order for interference between decays with and without B^-B^ mixing to occur, the X system 
must not be flavour-specific, e.g., in case of B^ K*^1i the final state must be i^g7r°7. The 
sign of the sine term depends on the C eigenvalue of the X system. The photons from b 
(b q'-f) are predominantly left (right) polarized, with corrections of order of mq/rrib, thus 
interference effects are suppressed. The predicted smallness of the S terms in the Standard 
Model results in sensitivity to new physics contributions. 



5.2.5 Asymmetries in B ^ decays 

CP asymmetries in i? — » _D(*)j'i'(*) decays are sensitive to 7. The neutral D*^*^ meson produced 
is an admixture of Z}(*)° (produced by a 6 — > c transition) and D^*^^ (produced by a colour- 
suppressed h ^ u transition) states. If the final state is chosen so that both and Z}(*)'^ can 
contribute, the two amplitudes interfere, and the resulting observables are sensitive to 7, the 
relative weak phase between the two B decay amplitudes. Various methods have been proposed 
to exploit this interference, including those where the neutral D meson is reconstructed as a 
CP eigenstate (GLW) |186j . in a suppressed final state (ADS) fTH7|, or in a self-conjugate 
three-body final state, such as K^ti^tt~ (Dalitz) jl88j . 

Consider the case of B^ DK^, with D decaying to a final state /, which is accessible to 
both and We can write the decay rates for B" and i?"*" (Tzp), the charge averaged rate 
(r = (r_ + r_|_)/2) and the charge asymmetry {A = (r_ — r_|_)/ (r_ + F^), see Eq. dHU)) as 



Tzp oc + + 2rBrD cos {63 + =F 7) , (76) 
r oc r\ + r'l) + 2rBrD cos [63 + Sd) cos (7) , (77) 
^ ^ ^rBTp sin {Sb + Sp) sin (7) 

^3 + ^0 + '^^bTd cos {5b + 5d) cos (7) , 



57 



where the ratio of B decay amphtudes^" is usually defined to be less than one, 

\A(B- ^ D°K-)\ , , 

and the ratio of D decay amplitudes is correspondingly defined by 

. (80) 

The strong phase differences between the B and D decay amplitudes are given by 5^ and 5/5, 
respectively. The values oivD and bo depend on the final state /: for the GLW analysis, td = \ 
and bo is trivial (either zero or vr), in the Dalitz plot analysis td and br, vary across the Dalitz 
plot, and depend on the D decay model used, for the ADS analysis, the values of rn and bu 
are not trivial. 

In the GLW analysis, the measured quantities are the partial rate asymmetry, and the 
charge averaged rate, which are measured both for CP even and CP odd D decays. For the 
latter, it is experimentally convenient to measure a double ratio, 

Y{B- ^DcpK-) IY{B- ^D'^K') 



T {B- ^ DcpTT-) / T {B~ ^ DH~) ' ' 

that is normalized both to the rate for the favoured A'~7r+ decay, and to the equivalent 

quantities for B~ — > Dn' decays (charge conjugate modes are implicitly included in Eq. (|HTj)). 
In this way the constant of proportionality drops out of Eq. (f77j) . 

For the ADS analysis, using a suppressed D f decay, the measured quantities are again 
the partial rate asymmetry, and the charge averaged rate. In this case it is sufficient to measure 
the rate in a single ratio (normalized to the favoured D f decay) since detection systematics 
cancel naturally. In the ADS analysis, there are an additional two unknowns (r^j and bo) 
compared to the GLW case. However, the value of r/) can be measured using decays of D 
mesons of known fiavour. 

The relations between the measured quantities and the underlying parameters are summa- 
rized in Table QUI Note carefully that the hadronic factors and bs are different, in general, 
for each B decay mode. In the Dalitz plot analysis, once a model is assumed for the D decay, 
which gives the values of ra and bo across the Dalitz plot, the values of (7, r^, bs) can be 
directly extracted from a simultaneous fit to the B" and 5+ data. 



5.3 Common inputs and error treatment 

The common inputs used for rescaling are listed in Table 1^ The B^ lifetime {t{B^)) and mixing 
parameter (Am^) averages are provided by the HFAG Lifetimes and Oscillations subgroup 
(Sec. EI). The fraction of the perpendicularly polarized component (|v4_l|^) in B ^ J/ipK*{892) 
decays, which determines the CP composition, is averaged from results by BABAR |19Uj and 
Belle [TM] . 

^°Note that here we use the notation to denote the ratio of B decay ampUtudes, whereas in Sec. 15. 2.11 we 
used, e.g., Ro-m for a rather similar quantity. The reason is that here we need to be concerned also with D 
decay amplitudes, and so it is convenient to use the subscript to denote the decaying particle. Hopefully, using 
r in place of R will help reduce potential confusion. 
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Table 29: Summary of relations between measured and physical parameters in GLW and ADS 
analyses of 5 ^ D^*'>K^*\ 




GLW analysis 
1 + r"^ ±2rB cos {6b) cos (7) 
±2rB sin {6b) sin (7) / Rcp± 



Rads 

^ADS 



ADS analysis 
^1 + + '^^bTd cos {6b + 5d) cos (7) 
2rBrD sin ((5b + 6d) sin (7) / Rads 



Table 30: Common inputs used in calculating the averages. 



r(fiO) (ps) 1.536 ±0.014 

Arrtd (ps"i) 0.502 ±0.007 
\Ax_\^ {J/ipK*) 0.211 ±0.011 



At present, we only rescale to a common set of input parameters for modes with reasonably 
small statistical errors (h — > ccs and h qqs transitions). Correlated systematic errors are 
taken into account in these modes as well. For all other modes, the effect of such a procedure 
is currently negligible. 

As explained in Sec. we do not apply a rescaling factor on the error of an average that has 
X^/dof > 1 (unhke the procedure currently used by the PDG jS]). We provide a confidence level 
of the fit so that one can know the consistency of the measurements included in the average, 
and attach comments in case some care needs to be taken in the interpretation. Note that, in 
general, results obtained from data samples with low statistics will exhibit some non-Gaussian 
behaviour. For measurements where one error is given, it represents the total error, where 
statistical and systematic uncertainties have been added in quadrature. If two errors are given, 
the first is statistical and the second systematic. If more than two errors are given, the origin 
of the third will be explained in the text. 

Averages are computed by maximizing a log-likelihood function C assuming Gaussian statis- 
tical and systematic errors. When observables are correlated {e.g., sine and cosine coefficients in 
time-dependent CP asymmetries), a combined minimization is performed, taking into account 
the correlations. Asymmetric errors are treated by defining an asymmetric log-likelihood func- 
tion: Ci = {x — Xi)^/{2af), where CTj = crj^+ (uj = (Ti-) if x > Xi {x < Xi), and where Xi is the 
ith measurement of the observable x that is averaged. This example assumes no correlations 
between observables. The correlated case is a straightforward extension to this. 

5.4 Time-dependent CP asymmetries in 6 — ^ ccs transitions 

In the Standard Model, the time-dependent parameters for b ccs transitions are predicted to 
be: Sb^ccs = — ?7sin(2/3), Cb^ccs = to very good accuracy. The averages for Sb^ccs and Cb^ccs 
are provided in Table |^ The averages for Sb^ccs are shown in Fig. averages for Cb^ccs are 
included in Fig. ^1 



Both BABAR and Belle use the r] = -1 modes J/i^K^, ip{2S)Kl, XaK^ and ricK^, as well 
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Table 31: Sb^ccs and Cb^ca 



Experiment 


V^b^ccs 


>ccs 


BABAR una 

Belle IHHl 


0.722 ±0.040 ±0.023 
0.728 ± 0.056 ± 0.023 


0.051 ±0.033 ±0.014 
-0.007 ± 0.041 ± 0.033 


B factory average 

Confidence level 


0.725 ± 0.037 
0.91 


0.031 ± 0.029 
0.30 


ALEPH 

OPAL dnu 

CDF 1951 


0.84 +^;«^±0.16 
3.2;^-^±0.5 

U. 1 J _o.44 




Average 


0.726 ±0.037 


0.031 ±0.029 


Table 32: Avera^ 


res from i?*^ ^ J/ipK*^ transversity analyses. 


Experiment 


sin(2/3) 


cos (2/5) 


BABAR dnni 

Belle [M] 


-0.10 ±0.57 ±0.14 
0.30 ±0.32 ±0.02 


3.32 +^;j^;?± 0.27 
0.31 ±0.91 ±0.11 


Average 


0.21 ±0.28 


1.69 ±0.67 



as J/ipKl, which has t] = ±1 and J/ip K*^ (892), which is found to have rj close to ±1 based on 
the measurement of \A±\ (see Sec. 15. 3^ . ALEPH, OPAL and CDF use only the J/ipK^ final 
state. 

These results give a precise constraint on the {j),f]) plane, in remarkable agreement with 
other constraints from CP conserving quantities, and with CP violation in the kaon system, in 
the form of the parameter ex- Such comparisons have been performed by various phenomeno- 
logical groups, such as CKMfitter |196j and UTFit |197j . Figure ITU displays the constraints 
obtained from these two groups. 

5.5 Time-dependent transversity analysis of — ^ J/il^K*^ 

B meson decays to the vector- vector final state J/ipK*^ are also mediated by the h — > ccs 
transition. When a final state which is not fiavour-specific [K*^ —)■ K^tt^) is used, a time- 
dependent transversity analysis can be performed allowing sensitivity to both sin(2/9) and 
cos(2/3). Such analyses have been performed by both B factory experiments. In principle, the 
strong phases between the transversity amplitudes are not uniquely determined by such an 
analysis, leading to a discrete ambiguity in the sign of cos(2/3). The BABAR |19()j collaboration 
resolves this ambiguity using the known variation |198j of the P-wave phase (fast) relative to 
the S-wave phase (slow) with the invariant mass of the Ktt system in the vicinity of the K*{892) 
resonance. The result is in agreement with the prediction from s quark helicity conservation, 
and corresponds to Solution II defined by Suzuki |199j . We use this phase convention for the 
averages given in Table ESI 

While the statistical errors are large, and exhibit non-Gaussian behaviour, cos(2/3) > is 
preferred by the experimental data in J/ipK*. 
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Figure 14: Standard Model constraints on the (p, 77) plane, from (top) |196j and (bottom) |197j . 
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5.6 Time-dependent CP asymmetries in 6 — > qqs transitions 

The flavour changing neutral current h ^ s penguin can be mediated by any up-type quark in 
the loop, and hence the amplitude can be written as 

= (F.-FJK^K*. + {Ft - F,)Vtf,V* (82) 
0{X^) + 0(A2) 

using the unitarity of the CKM matrix. Therefore, in the Standard Model, this amplitude is 
dominated by VtbV^l, and to within a few degrees {6(3^2° for (3 ^ 23.3°) the time-dependent 
parameters can be written^^ S^^q-qg ~ — ?7sin(2/3), Cb->qqs ~ 0, assuming b ^ s penguin contri- 
butions only (g = u, d, s). 

Due to the large virtual mass scales occurring in the penguin loops, additional diagrams 
from physics beyond the Standard Model, with heavy particles in the loops, may contribute. In 
general, these contributions will affect the values of Sb^qqg and Cb^qqs. A discrepancy between 
the values of Sb^ccs and Sb^qqs can therefore provide a clean indication of new physics. 

However, there is an additional consideration to take into account. The above argument 
assumes only the 6 — > s penguin contributes to the b qqs transition. For q = s this is a 
good assumption, which neglects only rescattering effects. However, for g = m there is a colour- 
suppressed 6 — > u tree diagram (of order 0{X^)), which has a different weak (and possibly 
strong) phase. In the case q = d, any light neutral meson that is formed from dd also has 
a uu component, and so again there is "tree pollution". The decays to iv^Kg and uiK^ 
belong to this category. The mesons /o and r]' are expected to have predominant ss parts, 
which reduces the possible tree pollution. If the inclusive decay K^K~K^ (excluding 

(pK^) is dominated by a non-resonant three-body transition, an OZI-rule suppressed tree-level 
diagram can occur through insertion of an ss pair. The corresponding penguin-type transition 
proceeds via insertion of a uu pair, which is expected to be favored over the ss insertion by 
fragmentation models. Neglecting rescattering, the final state K^K^K^ has no tree pollution. 

The averages for Sb^qqs and Cb^qqs can be found in Table ESI The averages for Sb^qqs are 
shown in Fig. ^1 averages for Cb^qqs are included in Fig. ^| Results from both BABAR and 
Belle are averaged for the modes (pK'^ (both (pK^ and 4>K^ are used), rj'K^, K~^K~Kg, foK^ 
and T^^Kg. In addition, results from Belle are taken for the modes ujK^ and K^K^Kg. Of these 
modes, (pK^ r]'Kl ti^K^ and uK^ have CP eigenvalue t] = -1, while /oK° and K^KlKl 
have ?7 = +1. 

The final state K^K~K^ (contributions from (pK^g are implicitly excluded) is not a CP 
eigenstate. However, the CP composition can be determined using either an isospin argument 
(used by Belle to determine a CP even fraction of 1.03±0.15±0.05 |2()()j ) or a moments analysis 
(used by BABAR who finds a CP even fraction of 0.89 ± 0.08 ± 0.06 [2DS]). The uncertainty in 
the CP even fraction leads to an asymmetric error on Sb^qqs, which is taken to be correlated 
among the experiments. To combine, we rescale the results to the average CP even fraction of 
0.93 ±0.09. 

If we treat the combined error as a Gaussian quantity, we note that the average of —rjSb^qqs 
of all b — > qqs dominated modes (0.41 ±0.07) is more than 5a from zero, and hence CP violation 

^^The presence of a small {0{X^)) weak phase in the dommant amplitude of the s penguin decays intro- 
duces a phase shift given by Sb-^qqs ~ — ?7sin(2/3) • (1 + A). Using the CKMfitter results for the Wolfenstein 
parameters j216| . one finds: A ~ 0.033, which corresponds to a shift of 2/3 of +2.1 degrees. Nonperturbative 
contributions can alter this result. 
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Figure 15: Sb^ccs and Sb 



*qqs- 



in b qqs transitions is established. Furthermore, the averages of —f]Sb^qqs for the modes 
rj'K^ and K~^K~Kg are more than 3a from zero. 

Neglecting theory errors due to suppressed contributions with different weak phases, the 
difference between Sb^ccs and Sb~tqqs can be calculated. We find the confidence level (CL) of 
the joint Sb^ccs and Sb^qqs average to be 0.00012, which corresponds to a 3.8a discrepancy. 
To give an idea of the theoretical uncertainties involved. Fig. shows coarse estimates of 
the theoretical errors associated with the non-charmonium modes. These crude estimates are 
obtained from dimensional arguments only, based on the CKM suppression of the Vub penguin, 
and on the naive contribution from tree diagrams. Including these estimates according to the 
procedure defined in Ref. |21(ij improves the CL of the joint average to about 3a. 



5.7 Time-dependent CP asymmetries in 6 — ^ ccd transitions 



The transition b 
Similarly to Eq. 



ccd can occur via either a 6 — * c tree or a. b ^ d penguin amplitude. 
, the amplitude for the b ^ d penguin can be written 



FuVubV:, + F^^,bV:, + F,v,bV:, 

{F^ - F,)V^bV:, + {F, ~ F,)V,bV*, 
0{X^) + C(A=^). 



(83) 



From this it can be seen that the b d penguin amplitude does not have a dominant weak 
phase. 
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Figure 16: Averages for Sb^qqs compared to Sb~*ccs, with coarse dimensional estimates of 
associated theoretical uncertainties indicated by the open boxes. 



In the above, we have followed Eq. (|H^ by eliminating the Fc term using unitarity. However, 
we could equally well write 

Ab^, = [F^ - F,)VubV:, + {F,-F,)V,bV:,, 

= (F, - F^)V,bV:, + (F, - F^)V,bV:,. ' 

Since the b — > ccd tree amplitude has the weak phase of VcbV*^, either of the above expressions 
allow the penguin to be decomposed into parts with weak phases the same and different to the 
tree amplitude (the relative weak phase can be chosen to be either /3 or 7). However, if the 
tree amplitude dominates, there is little sensitivity to any phase other than that from B^-B^ 
mixing. 

The b ccd transitions can be investigated with studies of various different final states. 
Results are available from both BABAR and Belle using the final states J/ipn^, D*^D*^ and 
£)*±^T- ^]-^g averages of these results are given in Table IH^ The results using the CP eigenstate 
(?7 = +1) J/ipTT^ are shown in Fig. El The vector- vector mode D*^D*~ is found to be dominated 
by the CP even longitudinally polarized component; BABAR measures a CP odd fraction of 
0.063 ± 0.055 ± 0.009 [210] while Belle measures a CP odd fraction of 0.19 ± 0.08 ± 0.01 pTT] 
(here we do not average these fractions and rescale the inputs, however the average is almost 
independent of the treatment). For the non-CP eigenstate mode D*^D^ BABAR uses fully 
reconstructed events while Belle combines both fully and partially reconstructed samples. 

In the absence of the penguin contribution, the time-dependent parameters would be given 
by = -r/sin(2/5), Cb-.c-cd = 0, = sin(2/3 + 6), 5_+ = sin(2/3 - 6), C+_ = 

and = 0, where 6 is the strong phase difference between the D*^D~ and D*^D~^ decay 
amplitudes. In the presence of the penguin contribution, there is no clean interpretation in 
terms of CKM parameters, however direct CP violation may be observed as any of Cb^ccd 7^ 0, 
C+_ or A+_ ^ 0. 
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Table 34: Averages for b ccd modes. 



Experiment 



b~*ccd 



b—fccd 



BABAR 
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Average 

Confidence level 
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Figure 17: Averages of (left) Sh^cai and (right) Cb^ccd for the mode B° Jjipn^. 



The averages for the h — > ccd modes are shown in Fig. Comparisons of the results for 
the h — s> ccd modes to the h — s> ccs and h qqs modes, can be seen in Fig. ^| 



5.8 Time-dependent asymmetries in 6 — ^57 transitions 

The radiative decays b ^ s'-f produce photons which are highly polarized in the Standard 
Model. The decays B° — > F7 and B*^ F'-f produce photons with opposite helicities, and 
since the polarization is, in principle, observable, these final states cannot interfere. The finite 
mass of the s quark introduces small corrections, to the limit of maximum polarization, but 
any large mixing induced CP violation would be a signal for new physics. Since a single weak 
phase dominates the h ^ S'~^ transition in the Standard Model, the cosine term is also expected 
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Figure 18: Averages of (left) Shaded and (right) Cb-,ccd, compared to Sb^ccs 
tively. 
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Table 35: Averages for 6 ^ 57 modes. 



Experiment 



St 



BABAR 

Belle 

Average 

Confidence level 



0.25 ±0.63 ±0.14 

0.46 



-0.58;-^^ ±0.11 



-0.57 ±0.32 ±0.09 
0.03 ±0.34 ±0.11 
-0.29 ±0.38 -0.32 ±0.24 

combined average: 0.31 (l.Ofj) 



to be small. 

The BABAR collaboration has studied time-dependent asymmetries in 6 — > 57 transitions 
using the decay -^*7; with K* — > i^gVr^. The Belle collaboration also uses the decay 

B^ — > K^gTT^'-), but do not restrict themselves to the K* mass region |185j (the invariant mass 
range 0.60 GeV/c^ < M^o^o < 1.80 GeV/c^ is used). These results, and their averages, are 
given in Table ESI 



5.9 Time-dependent CP asymmetries in 6 — > uud transitions 

The h — *• uud transition can be mediated by either a 6 — > u tree amplitude or a 6 — > d penguin 
amplitude. These transitions can be investigated using the time dependence of B^ decays to 
final states containing light mesons. Results are available from both BABAR and Belle for 
the CP eigenstate (r/ = ±1) tt+tt^ final state. BABAR has also performed an analysis on the 
vector- vector final state p~^p^ , which they find to be dominated by the CP even longitudinally 
polarized component (they measure /long = 0.99 ± 0.03 ^Q^i, |220p . 

For the non-CP eigenstate p^vr^. Belle has performed a quasi- two-body analysis, while 
BAByiR performs a time-dependent Dalitz plot (DP) analysis of the 7r+7r~7r° final state |217j : 
such an analysis allows direct measurements of the phases. These results, and averages, are 
listed in Table l36l The averages for 7r~^7T~ are shown in Fig. 1201 

If the penguin contribution is negligible, the time-dependent parameters for B° tt+tt" 
and B° P~^P~ are given by Sb^und = V sin(2a) and C^^uud = 0. With the notation described 
in Sec. 15.21 fEq. (|69p). the time-dependent parameters for B" p^Tr"^ are, neglecting penguin 

contributions, given by SpT, = y^l — (^)2sin(2a) cos((5), ASp^r = y''l — (^)2cos(2a) sin((5) 

and CpT, = A'^p = 0, where 6 = a.Tg{A_^A'^_) is the strong phase difference between the 
p~7i~^ and p"^7r~ decay amplitudes. In the presence of the penguin contribution, there is no 
straightforward interpretation of B° p'^n^ in terms of CKM parameters. However direct 
CP violation may arise, resulting in either or both of Cp^ 7^ and A'^p ^ 0. Equivalently, 
direct CP violation may be seen by either of the decay-type-specific observables Ap~ and A~J', 
defined in Eq. (I70|) . deviating from zero. Results and averages for these parameters are also 
given in Table EHl They exhibit a linear correlation coefficient of ±0.59. The significance of 
observing direct CP violation computed from the difference of the obtained in the nominal 
average, compared to setting Cp^ = A'^p = is found to be 3.4cr in this mode. The confidence 
level contours of A'^~ versus A~J~ are shown in Fig. EU 

Some difference is seen between the BABAR and Belle measurements in the 7r"'"7r~ system. 
The confidence level of the average is 0.0014, which corresponds to a 3.2cr discrepancy. Since 
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Table 36: Averages for b uud modes. 
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Figure 20: Averages of (left) Sb^und and (right) Cb^uud for the mode 5° n^n^ . 
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Figure 21: Direct CP violation in — p^ir^. The no-CP violation hypothesis is excluded at 
the 3.4(7 level. 

there is no evidence of systematic problems in either analysis, we do not rescale the errors of 
the averages. 

The precision of the measured CP violation parameters in 6 — > uud transitions allows 
constraints to be set on the UT angle a. In addition to the value of a from the BABAR 
time-dependent DP analysis, given in Table EEl constraints have been obtained with various 
methods: 

• In the Belle analysis of B^ tt+tt" |219j . constraints on the time-dependent CP violation 
parameters are used to obtain 90° < a < 146° (95.5% CL) following the method proposed 
by Gronau and Rosner |224j . This result includes an assumption on the relative size of 
tree and penguin amplitudes. 

• Using the measured time- dependent CP violation parameters in longitudinally polar- 
ized B^ P~^P~ decays |220j . in combination with the upper limit for the B^ — > p^p^ 
branching fraction |225j . and the measurement of the branching fraction and longitudi- 
nal polarization of B~^ — > p'^p^ |226| I227j . BABAR performs an isospin analysis |228j and 
obtains a = (96 ± 10 ± 4 ± 11)°, where the third error is due to the unknown penguin 
contribution. 

• The CKMfitter group |196j uses the measurements from Belle and BABAR given in Ta- 
ble inni with other branching fractions and CP asymmetries in B ^ irir, pir and pp 
modes, to perform isospin analyses for each system. They then combine the results to 
obtain a = {lOOtilY- 

Note that each method suffers from ambiguities in the solutions. All the above measurements 
correspond to the choice that is in agreement with the global CKM fit. 

At present we make no attempt to provide an HFAG average for a. More details on proce- 
dures to calculate a best fit value for a can be found in Refs. jl96| I197j . 







- 3a 
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Table 37: 


Averages for b cud/ucd modes. 




Experiment 




a 


c 


BABAR (full rec.) 
Belle (full rec.) 
BABAR (partial rec.) 
Belle (partial rec.) 
Average 
Confidence level 


[222] 


D*V 
-0.049 ±0.031 ±0.020 
0.060 ±0.040 ±0.019 
-0.041 ±0.016 ±0.010 
-0.031 ±0.028 ±0.018 
-0.030 ±0.014 
0.19 


0.044 ±0.054 ±0.033 
0.049 ±0.040 ±0.019 
-0.015 ±0.036 ±0.019 
-0.004 ±0.028 ±0.018 
0.010 ±0.021 
0.66 


BABAR (full rec.) 
Belle (full rec.) 
Average 

Confidence level 


[22i 
12201 


-0.032 ±0.031 ±0.020 
-0.062 ±0.037 ±0.018 
-0.045 ±0.027 
0.59 


-0.059 ±0.055 ±0.033 
-0.025 ±0.037 ±0.018 
-0.035 ±0.035 
0.66 


BABAR (full rec.) 


[22i 


-0.005 ±0.044 ±0.021 


-0.147 ±0.074 ±0.035 



5.10 Time-dependent CP asymmetries in 6 — ^ cud/ucd transitions 

Non-CP eigenstates such as D'^tt^, D*^7r^ and D'^p^ can be produced in decays of B^ mesons 
either via Cabibbo favoured (6 — > c) or doubly Cabibbo suppressed (6 — > u) tree amplitudes. 
Since no penguin contribution is possible, these modes are theoretically clean. The ratio of the 
magnitudes of the suppressed and favoured amplitudes, R, is sufficiently small (predicted to be 
about 0.02), that terms of 0{R^) can be neglected, and the sine terms give sensitivity to the 
combination of UT angles 2(3 ± 7. 

As described in Sec. I5.2.4| the averages are given in terms of parameters a and c. CP 
violation would appear as a 7^ 0. Results are available from both BABMi and Belle in the 
modes D'^n^ and D*^7r^; for the latter mode both experiments have used both full and partial 
reconstruction techniques. Results are also available from BABAR using D^p^ . These results, 
and their averages, are listed in Table [23 and are shown in Fig. |221 

5.11 Rates and asymmetries in — > decays 

As explained in Sec. 15.2.51 rates and asymmetries in B^ decays are sensitive to 

7. Various methods using different D^*^ final states exist. 

Results are available from both BABAR and Belle on GLW analyses in the decay modes 
B=F DK^, B^ D*K^ and B^ DK*^. Both experiments use the CP even D decay 
final states K^K~ and 7r"*"7r~ in all three modes; both experiments also use only the D* Dir^ 
decay, which gives CP{D*) = CP{D). For CP odd D decay final states. Belle uses Kgir'^, K^r] 
and Kg(j) in all three analyses, and also use K^uj in DK^ and D*K^ analyses. BABAR uses 
A'gTT^ only for DK^ analysis; for DK*^ analysis they also use Kg(j) and K^uj (and assign an 
asymmetric systematic error due to CP even pollution in these CP odd channels |237j ). The 
results and averages are given in Table IHHl 

For ADS analysis, both BABAR and Belle have studied the mode B^ DK^, and BABAR 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
BABAR (partially reco.) , . , 

-0.041 ±0.016±0. 010 ' " ' 

BABAR (fully reco.) 

-0.04910.031 ±0.020 ' * ' 

a * Belle (partially reco.) , 

" -0.03110.028 + 0.018 ' " 

Belle (fully reco.) 

0.060 ±0.040 ±0.01 9 

Average , . , 
-0.030±0.014 '"*"' 




1 • 1 


BABAR (partially reco.) 

-0.015±0.036±0.019 ' ' 

BABAR (fully reco.) 

0.044 ± 0.054 ±0.033 1 

r. * Belle (partially reco.) , _ 

" -0.004 ±0.028 ±0.01 8 ' ^ 

Belle (fully reco.) 

0.049 ±0.040 ±0.01 9 

Average , 

0.010±0.021 ^ 


1 

■ 1 

1 

1 ■ 1 

■A 1 


BABAR (fully reco.) 

-0.032±0.031 ±0.020 ' * 

a Belle (fully reco.) 

" -0.062±0.037±0.018 ' * ' 

Average 

-0.045±0.027 ' ^ ' 


H 


BABAR (fully reco.) , _ 

-0.059 ±0.055 ±0.033 ' " 

r. Belle (fully reco.) , _ 

" -0.025±0.037±0.018 ' " 

Average , 

-0.035 ±0.035 ' ^ 




a BABAR (fully reco.) , _ 

P -0.005 ± 0.044 ±0.021 ' ■ 


BIhfagI 


r. BABAR (fully reco.) , 


1 FRCP 2004 


P -0.1 47 ±0.074 ±0.035 ' " ' 





-0.4 -0.3 -0.2 -0.1 0.1 



Asymmetry in D**'jc{p) decays 
Figure 22: Averages for b cud/ucd modes. 



Table 38: Averages from GLW analyses of 6 ^ cud/ucd modes. 



Experiment 


Acp+ 


AcP- 


Rcp+ 


Rcp- 






DcpK- 






BABM. PSH] 


0.40 ±0.15 ±0.08 


0.21 ±0.17 ±0.07 


0.87 ±0.14 ±0.06 


0.80 ±0.14 ±0.08 


Belle I22S1 


0.07 ±0.14 ±0.06 


-0.11 ±0.14 ±0.05 


0.98 ±0.18 ±0.10 


1.29 ±0.16 ±0.08 


Average 


0.22 ± 0.11 


0.02 ±0.12 


0.91 ±0.12 


1.02 ±0.12 


BABAR PITTI 


-0.02 ±0.24 ±0.05 


DhpK- 


1.09 ±0.26 




Belle 


-0.27 ±0.25 ±0.04 


0.26 ±0.26 ±0.03 


1.43 ±0.28 ±0.06 


0.94 ±0.28 ±0.06 


Average 


-0.14 ±0.18 


0.26 ±0.26 


1.25 ±0.20 


0.94 ±0.29 






DcpK*- 






BABAR I22H1 


-0.09 ±0.20 ±0.06 


-0.33 ± 0.34 ±0. 10 


1.77 ±0.37 ±0.12 


0.76 ±0.29 ±0.0612:?^ 


Belle !^ 


-0.02 ±0.33 ±0.07 


0.19 ±0.50 ±0.04 ' 






Average 


-0.07 ±0.18 


-0.16 ± 0.29 
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Table 39: Averages from ADS analyses of 6 ^ cud/ucd modes. 

Experiment Aads Raps 

DK-, D K+TT" 

BABAR [2101 0.013 ^SJJ^^ 

Belle 1211] 0.49 ;|]j^± 0.06 0.028 ± 0.010 
Average 0.49 tpjl 0.017 ±0.009 

D*K-, D* Dtt^, D ^ K+Tv- 
BABAR ^ -0.001 

D*K-, D* ^ D-i,D ^ K+TT- 
BABAR 12101 0.011 ^SJJ}^ 



Table 40: Averages from Dalitz plot analyses of 6 ^ cud/ucd modes. 

Experiment 7 (°) 6b (°) rs 

DK~, D fsr^vr+TT- 

BABAR 1212] 70 ±44 ±10 ±10 114±41± 8 ± 10 < 0.19 @ 90% CL 

Belle 1213 64 ±19 ±13 ±11 157 ±19 ±11 ±21 0.21 ± 0.08 ± 0.03 ± 0.04 
Average IN preparation 

D*K~, D* Dir^ or D-f, D K^^tt+ti- 

73±35± 8 ± 10 303 ±34 ±14 ±10 0.16 ± 0.04 ± 0.02 

75 ±57 ±11 ±11 321 ±57 ±11 ±21 0.12 +°:}^ ± 0.02 ± 0.04 

IN PREPARATION 

DK^ and D*K^ combined 
BABAR 1212] 70 ± 26 ± 10 ± 10 
Belle 1213 68lJ^±13±ll 

Average IN preparation 



BABAR 1212] 
Belle I21S] 
Average 



has also analyzed the D*K^ mode [D* Dn^ and D* are studied separately). 

In all cases the suppressed decay D — > K^n" has been used. The results and averages are 
given in Table EOI Note that although no clear signals for these modes have yet been seen, 
the central values are given in Table |201 In B~ D*K~ decays there is an effective shift of 
71 in the strong phase difference between the cases that the D* is reconstructed as Dir^ and 
|189j . As a consequence, the different D* decay modes are treated separately. 

For the Dalitz plot analysis, both BABAR and Belle have studied the mode B^ DK^ . 
Both have also studied the mode B^ — > D*K^] Belle has used only D* —>■ Dir^, while BABAR 
has used both D* decay modes and taken the effective shift in the strong phase difference into 
account. In all cases the decay D A'^tt^tt" has been used. The results are given in Table EOl 
Since the measured values of are positive definite, and since the error on 7 depends on 
the value of tb, some statistical treatment is necessary to correct for bias. Belle has used a 
frequentist treatment, while BABAR used a Bayesian approach. At present, we make no attempt 
to average the results. 



73 



6 Averages of charmless B-decay branching fractions and 
their asymmetries 

The aim of this section is to provide the branching fractions and the partial rate asymmetries 
{Aqp) of rare B decays. The asymmetry is defined as Acp = Ty^^nvf ' where and A''^ are 

number of / B~ and B^/B~^, respectively. Four different B decay categories are considered: 
charmless mesonic, baryonic, radiative and leptonic. Rare mesonic decays with charm are not 
in our scope but results of charmful baryonic decays are included. Measurements supported 
with written documents are accepted in our the averages; written documents could be journal 
papers, conference contributed papers, preprints or conference proceedings. Results from Aqp 
measurements obtained from time dependent analyses are listed and described in Sec. El 

So far all branching fractions assume equal production of charged and neutral B pairs. 
The best measurements to date show that this is still a good approximation (see Sec. 13. Lip . 
For branching fractions, we provide either averages or the most stringent 90% confidence level 
upper limits. If one or more experiments have measurements with >4cr for a decay channel, all 
available central values for that channel are used in the averaging. We also give central values 
and errors for cases where the significance of the average value is at least 3a, even if no single 
measurement is above 4a. For Acp we provide averages in all cases. 

Our averaging is performed by maximizing the likelihood. 



where Vi is the probability density function (PDF) of the ith measurement, and x is the 
branching fraction or Aqp- The PDF is modeled by an asymmetric Gaussian function with the 
measured central value as its mean and the quadratic sum of the statistical and systematic errors 
as the standard deviations. The experimental uncertainties are considered to be uncorrelated 
with each other when the averaging is performed. No error scaling is applied when the fit is 
greater than 1 since we believe that tends to overestimate the errors except in cases of extreme 
disagreement (we have no such cases). 

At present, we have measurements of 230 B decay modes and asymmetry measurements 
for 35 of these decays. These results are reported in 112 separate papers. Because the number 
of references is so large, we do not include them with the tables shown here but the full set 
of references is available quickly from active gifs at the ICHEP04 link on the rare web page: 
http : //www. slac . Stanford. edu/xorg/hfag/rare/ index .html 

6.1 Mesonic charmless decays 
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Table 41: branching fractions (in units of 10 Upper limits are at 90% CL. Values in 
red (blue) are new published (preliminary) result since PDG2004 [as of August 25th, 2004]. 



RPP# 


Mode 


PDG2004 Avg. 


BABAR 


Belle 


CLEO 


CDF 


New Avg. 


117 




18.8 ±2.1 


26.0 ±1.3 ±1.0 


22.0 ±1.9 ±1.1 


IS fi+3.7+2.1 




24.1 ± 1.3 


118 


K+tt" 


12.9 ± 1.2 


12.0 ± 0.7 ± 0.6 


12 ± 1 3+i'n 


12 9+2'4+l-? 




12.1 ± 0.8 


1 1 Q 


Tj I\ 


7s -1- 
/ m 


(D.y ZLL 0.0 III 


/ ^ D ^ y 


on"^-'-'^ _|_ 7 

OU__g III / 




77 

< '■D-4.5 








^ ^ A 


^ on 


■\ 00 




^ 1 A 
<^ 14 


121 


nK+ 


< 6.9 


"^4-1-0^^-1-09 
0.4 ^ u.o ^ u.^ 




Q Q + 2.8 

■^■'^ — 2.2 




2.6 ± 0.5 


122 








99 8+^'^ -1-9 9 
^^.0 3 ^ HI • 


9fi 4+9-6 _|_ 




94 S+3" 




ciQ\i)ou)i\ 1 




^ 9 






^ 9 ^ 






New 


< o.y 








< O.y 


123 




9.2t^j 


4.8 ± 0.8 ± 0.4 


, , 3 

6.5_5^ 2 ± 0.6 


3.2_j g ± 0.8 




5.1 ± 0.7 


124 


ljK*+ 


< 87 


< 7.4 




< 87 




< 7.4 


125 


K''°n+ 




10.5 ±2.0 ±1.4 


9.83 ± 0.90Ii;24 


7.6I30 ± 1.6 




9-7611:22 


126 


K* + TT° 


< 31 






< 31 




< 31 


127 


K+n+n- 


57 lb 4 


61.4 ±2.4 ±4.5 


46.6 ±2.1 ±4.3 






53.5 ± 3.5 


128 


K + TT+TT~{NR) 


< 28 


4.9 ± 0.6 ± 1.4 




< 28 




4.9 ± 1.5 


129 


K+M980) t 


seen 


9.2 ±1.5 ±0.8 


7.55 ± 1.24^^;^^ 






8.49ll:i 


130 




< 12 


5.2 ±1.2 ±0.7 


4.78 ± 0.75tJ;g7 


^ . -1-4 n , ^ 

8.413;^ ± 1.8 




_ Lfl Q1 

■J-^^-o.sg 




/2(1270)ii:+ t 


New 




< 1.3 






< 1.3 




/2(1525)K^ t 


New 




< 2.1 






< 2.1 




-TVq -L'iOU^ TV 1 


New 


QOQ^l Q-Lor: 

oz.o zii -L.y in z.o 








Q9 Q + Q 1 

oz.o ^0.1 


1 "^i 


-£V2 ^J-'iOU^ TV 1 


, - f::Qn 
<^ UoU 




^ 9 






^ z.o 






New 




^ 0. -L 








1 'iO 
±OZ 


ft' — ■7T- + -7r + 
I\ TV Ti 


<^ i.O 


^ 1.0 


^ A ^ 






<^ 1.0 


135 




< 66 






< 66 




< 66 


136 


K°o+ 


< 48 






< 48 




< 48 






New 


1 7 n -1- 2 Q+20 

i- t -yj _i_ ^ .zj si 


\J • VJ _1_ Zj.^ _i_ vj.o 






9.2 ± 2.0 


138 




11 lb 4 


10.6+2'5 ± 2.4 




< 74 




10 6+?? 


139 


K ^ K ^ 


< 71 






< 71 




< 71 


142 




< 2.0 


< 2.4 


< 3.3 


< 3.3 




< 2.4 


143 




< 24 






< 24 




< 24 


144 


R-^KsKs 


13.4 lb 2.4 


10.7 ± 1.2 ± 1.0 


Id. 4 ± 1.9 ± 1.0 






11.5 ± 1.3 


I'^O 


Rs^s'^'^ 












o.z 




i\ i\ TT 


<s. D.O 


<^ D.O 








^ fi 

*^ D.O 


148 




< 1.3 


< 1.3 


< 2.4 






< 1.3 


150 




< 5.3 






< 5.3 




< 5.3 


152 


K'^ K~ 


KjyJ . _1_ . -L 


90 f! + 9 1 -1- 1 R 

iiC'.U _l_ . i _1_ -L . "J 


■jn fi -1- 1 2 + 23 






30.1 ± 1.9 






Q Q _i_ 1 n 
y .0 ^ 1 .u 


1 n n+'~*-^ -1- n 

lU.U Q g KJ.O 


Q fin + n Q9+1-05 
y.ou i u.yz jj 34 


t- c; + 2.1 + n 

0.0 -j^ 3 in u.u 


( .z in 1 .0 in u. 1 


Q 074-0.65 
-0.63 




^ IV 4- + 
12 T 


New 




< i.i 






< 1.1 




</>(1680)i^+t 


New 




< 1.1 






< 1.1 


156 




9.6 lb 3.0 


12.7"'"on lb 1.1 


+2.1+0.7 
' -1.9-1.0 


10.6+5-^+li 

^^—4.9 — 1.6 




9.7 ± 1.5 


159 


§ 






2.6lJ:g ±0.3 






2-6lJ:J 


173 




5.6+?-? 
— 1 . i 


5.8 ±0.6 ±0.4 


5.0 ±1.2 ±0.5 


.g+1. 8+0.6 
■ —1.6 — 0.7 




5.5 ±0.6 


174 




11 ±4 


16.2 ± 2.1 ± 1.3 






16.2 ± 2.5 


175 




8.6 ±2.0 


9.4 ±1.3 ±1.0 


8.0+H ±0.7 


10.4+H ± 2.1 

— 0.4 




9.1 ± 1.3 




p" (1450)77+ 


New 


2.2 ±0.5 ±0.3 




2.2 ±0.6 


176 


7r+/o(980) t 


< 140 










< 140 


177 


/2 (1270)77+ 


< 240 


2.3 ±0.4 ±0.3 








2.3 ±0.5 


180 


p+7r0 


< 43 


10.9 ±1.9 ±1.9 


13.2±2.3li:9 


< 43 




12.0 ±2.0 


182 


p+pO 


26 ± 6 


22.5t5 4 ± 5.8 


31.7±7.ll|:? 






26.416;! 


185 


a;7r+ 




5.5 ±0.9 ±0.5 


5.7ti;3 ±0.6 


11.312:9 ± 




5.9 ±0.8 


186 


ajp+ 


< 61 


12.613 3 ± 1-6 


< 61 




12.6lti 


187 


T)7r+ 


< 5.7 


5.3 ±1.0 ±0.3 


4.8 ±0.7 ±0.3 


1 9-1-2.8 
^■^-1.2 




4.8 ±0.6 


188 


»)'7r+ 


< 7 


4.2 ±1.0 ±0.5 


< 7 


^•'^-1.0 




4.2 ± 1.1 


189 


Vp+ 


< 33 


< 22 




< 33 




< 22 


190 


r;p+ 


< 15 


8.6 ± 2.2 ± 1.1 


8.5l^;| ± 1.0 






8-611:^8 




a^(980)7r+ f 


New 


< 5.8 








< 5.8 


191 


(/)7r+ 


< 0.41 


< 0.41 




< 5 




< 0.41 


192 




< 16 






< 16 




< 16 



fProduct BF - daughter BF taken to be 100%; §M^^ < 2.85 GeV/c^ 
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Table 42: branching fractions (in units of 10 ^). Upper limits are at 90% CL. Values in 



red (blue) are new published (preliminary) result since PDG2004 [as of August 25th, 2004]. 



RPP* 


Mode 


PDG2004 Avg. 


BABAR 


Belle 


CLEO CDF 


New Avg. 


123 


K + TT 


18.5 ± 1.1 


17.9 ± 0.9 ± 0.7 


18.5 ± 1.0 ± 0.7 


„ ri+2.3+1.2 

18.0t2.?+o.9 


18.2 lb 0.8 


124 




9-5-1.9 


11.4 ± 0.9 ± 0.6 


11.7 ± 2.3l;j'3 


1 9 O+4.0+1.7 
^^■"-3.3-1.4 


11.5 di 1.0 


125 




63 ±7 


60.6 ±5.6 ±4.6 


68 ± lOIg 


89Iig ± 9 


65.21^-0 


126 


ri'K*° 


< 24 


< 7.6 


< 20 


< 24 


< 7.6 


127 


•qK " 


1 ,1+6 


18.6 ± 2.3 ± 1.2 


19.8t|J ± 1.4 


13.811' g ± 1.6 


18.7 ± 1.7 


128 




< 9.3 


2.5 ± 0.8 ± 0.1 


O.Sloy ± 0.1 




1.5 ± 0.6 






New 




^+'3 54-2 1 


oo .1+4 1 


— 


a~(980)X+ t 


New 


< 2.1 


< 2.9 




< 2.1 




a{](980)X'' t 


New 


< 7.8 






< 7.8 


129 




< 13 


5.9li;3 ±0.5 


4.0+1 6 ±0.5 


10. ol^:^ ± 1.4 


5.511;? 


131 




< 23 


< 6.0 




< 23 


< 6.0 


132 


K+K~ 


< 0.6 


< 0.6 


< 0.7 


< 0.8 


< 0.6 


133 


RO'ko 


< 3.3 


1.19l°;«±0.13 


< 1.5 


< 3.3 


1 IQ+0'12 
^•^='-0.37 


134 


KsKsKs 


— 1.5 


6.5 ±0.8 ±0.8 


4.2ti;3 =•= 0-8 




5.8 ±1.0 


135 


K+TT--K° 


< 40 


34.9 ±2.1 ±3.9 


36. 6+41 i 3-0 


< 40 


35.6+3-1 

— O.O 


136 


K+p- 


7.3 ± 1.8 


8.6 ± 1.4 ± 1.0 


1+3.4+2.4 
±o.± 2 3 2 6 


16lg ±3 




— 


K+p{U50)- t 


New 


< 3.2 




< 3.2 


— 


ii:+p(1700)- t 


New 


< 1.7 




501^° ± 7 


< 1.7 


137 




47 ±7 


43.7 ±3.8 ±3.4 


45.4 ±5.2 ±5.9 


44.9 ±4.0 




K + -K-TT°(NR) 


New 


< 4.6 


< 9.4 


< 4.6 




is:o(1430) + 7r- t 


New 


11.2 ±1.5 ±3.5 






11.2 ±3.8 


— 




New 


7.9 ±1.5 ±2.7 






7.9 ±3.1 


138 




< 39 


5.1 ± 1.0± 1.2 


< 12.4 


< 39 


5.1 ± 1.6 


139 


KO/o(980) t 


< 36 


6.0 ± 0.9 ± 1.3 


< 14 




6.0 ± 1.6 


140 




16« 


11.9±1.7±1.1 


14 c+4.6+2.8 
14.8_4.4_i.3 


16l^ ± 2 


12.7t\i 


141 




< 3.6 


3.0 ±0.9 ±0.5 




n n+i-3+0.5 


1.7 ±0.8 


142 


is:|(1430) + 7r- 


< 18 


< 13.2 




< 13.2 






New 


< 3.6 






< 3.6 


- 




New 


< 19.4 






< 19.4 


- 


K*(1680)07r0 


New 


< 5.0 






< 5.0 


143 




< 21 




< 18 


< 21 


< 18 


144 


K+K-iT° 


< 19 






< 19 


< 19 


145 


K+K-K° 


28 ±5 


23.8 ± 2.0 ± 1.6 


28.3 ±3.3 ±4.0 




24.7 ± 2.3 


146 


<j>K° 






9.0ti g ± 0.7 


5.4t^;^ ± 0.7 


8.3ll:g 


149 




< 34 




< 2.6 


< 34 


< 2.6 


- 


K*+p- 


New 


< 24 






< 24 


154 


<f>K*° 


10.7 ± 1.1 


9.2 ±0.9 ±0.5 




11 c;+4.5+1.8 
^^■'-'-3.7-1.7 


9.5 ±0.9 


155 




< 22 






< 22 


< 22 


157 


K*+K*- 


< 141 






< 141 


< 141 


176 




4.8 ±0.5 


4.7 ±0.6 ±0.2 


4.4 ±0.6 ±0.3 


. ^+1.4+0.5 
— 1.2— 0.4 


4.6 ±0.4 


177 




1.9 ±0.5 


1.17 ±0.32 ±0.10 


r, qo+0.44+0.22 
— u.^o — U.J-O 


< 4.4 


1.51 ±0.28 


178 




< 2.9 


< 2.5 


< 2.5 


< 2.9 


< 2.5 


179 


rjT} 


< 18 


< 2.8 


< 2.0 


< 18 


< 2.0 


180 




< 5.7 


< 3.7 




< 5.7 


< 3.7 


181 


rj' rj' 


< 47 


< 10 




< 47 


< 10 


182 


r/'q 


< 27 


< 4.6 




< 27 


< 4.6 


183 


Vp" 


< 12 


< 4.3 


< 14 


< 12 


< 4.3 


184 




< 10 


< 1.5 


< 5.5 


< 10 


< 1.5 




7J7r~^ TV 


Ngw 




^"■"—3.2—1.0 




16 6+^-^ 




.-f"F /'Q^^^^'7r=t + 


New 


^ ^ 1 

^ 0.1 


<^ o.o 




^ o.o 


185 


LOT] 


< 12 


< 2.3 




< 12 


< 2.3 


186 


UJT]' 


< 60 


< 2.8 




< 60 


< 2.8 


187 


ujp° 


< 11 


< 3.3 




< 11 


< 3.3 


189 




< 5 


< 1.0 




< 5 


< 1.0 


190 


0»? 


< 9 


< 1.0 




< 9 


< 1.0 


191 


4>ri' 


< 31 


< 4.5 




< 31 


< 4.5 


192 




< 13 






< 13 


< 13 


194 


H 


< 12 


< 1.5 




< 12 


< 1.5 


196 


pOwO 


< 5.3 


< 2.9 


5.1 ± 1.6±0.9 


< 5.5 


< 2.9 


197 




22.8 ±2.5 


22.6 ±1.8 ±2.2 


29.lt4 9 ±4.0 


27.61?-:^ ±4.2 


24.0 ±2.5 


199 


pOpO 


< 2.1 


< 1.1 




<'l8 


< 1.1 


200 




< 490 


42.6 ±4.2 ±4.1 






42.6 ±5.9 


203 


p+p- 


< 2200 


30 ± 4 ± 5 






30 ±6 


205 




< 3 


< 1.2 


< 1.9 


< 5.5 


< 1.2 



fProduct BF - daughter BF taken to be 100% 



6.2 Radiative and leptonic decays 



Table 43: Compilation of semileptonic and radiative branching fractions (in units of 10 ^). 
Upper limits are at 90% CL. values in red (blue) are new published (preliminary) result since 
PDG2004 [as of August 25th, 2004]. 





iViOQC 




Tin Roi. 


oCiiG 




New Avg. 


160 


K (8y2J^7 


38 lb 5 


QO v -1— o o _l_ o a 
oo. t ± 2.5 ± 2.6 


/I O K _1_ Q 1 -1— O A 

4z.o it o.i ± 2.4 


' •D_g 3 it 2.5 


4U.O ± z.b 


161 




< 99 




4.9Q-I-Q4-1-4Q 




49 Q _j_ 1 n 

1^.0 -LU.O 






Q /[ _|_ 1 fl 




o.^ HI u.y ^ u.^ 




Q yl _|_ 1 n 
0.'^: Ill l.U 


1 dQ 

loo 


rr^^ ^, C 

iv ' 7r 7r ' 7 5 






ZO III l.o ± z.z 




zo.U iL Z.o 


164 




201^ 




20lg ± 2 




20l6 


165 


K+p'^-l § 


< 20 




< 20 




< 20 


166 


_fC+7r~7r+7 (N.R.) § 


< 9.2 




< 9.2 




< 9.2 


167 


A'i(1400)+7 


< 50 




< 14.4 




< 14.4 


168 


i(-|(1430) + 7 


< 1400 


14.5 ±4.0 ±1.5 






14.4 ±4.2 


172 


p+7 


< 2.1 


< 1.8 


< 2.2 


< 13 


< 1.8 


207 


pyl7 


New- 




2.16t° g3 ±0.20 




^■-'-"-0.57 


208 


pX«7 


New 




< 3.3 




< 3.3 


- 


TT+l'F 


New 


< 100 






< 100 


226 




0-63l°;l? 


l.OSlg 22 ±0.07 


0.63t° i7 ± 0.03 


< 2.4 


0.80 ±0.15 


227 


K+fi+ii- 




0.07to 11 ±0.02 


0.451:^ ± 0.03 


< 3.68 


0.34 ±0.10 


229 


K+vV 


< 240 


< 52 




< 240 


< 52 


230 


K*{%92)+e+e- 


< 4.6 


0.20j;o;g7 ± 0.28 t 


r, „r, + 1.27+ 0.23 ^ 

^■^^-1.01-0.24 + 




. 20+0-90 

-'-■^^-0.77 


231 


K" (892)+ 11+ fi- 


< 2.2 


3.071 i ll ± 0.42 t 


n cc:+0.69+0.14 ^ 
U.OO_Q 53_Q 15 4. 




n Q2+0-70 
^■^^-0.58 


238 


■K~e+e+ 


< 1.6 






< 1.6 


< 1.6 


239 


n~ ij,+fj,+ 


< 1.4 






< 1.4 


< 1.4 


240 


■K~ e+ IX+ 


< 1.3 






< 1.3 


< 1.3 


241 


p~e+e+ 


< 2.6 






< 2.6 


< 2.6 


242 


p- 11+^1+ 


< 5.0 






< 5.0 


< 5.0 


243 


p~ e+ ^+ 


< 3.3 






< 3.3 


< 3.3 


244 


K'e+e+ 


< 1.0 






< 1.0 


< 1.0 


245 


K-p+fi+ 


< 1.8 






< 1.8 


< 1.8 


246 


K'e+^i+ 


< 2.0 






< 2.0 


< 2.0 


247 


K'-e+e+ 


< 2.8 






< 2.8 


< 2.8 


248 


K*~fj,+^i+ 


< 8.3 






< 8.3 


< 8.3 


249 


K*-e+p+ 


< 4.4 






< 4.4 


< 4.4 



§ Mktvtv < 2.4 GeV/c^ f Central values are not significant. 
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Table 44: Compilation of semileptonic and radiative branching fractions (in units of 10 ^). 
Upper limits are at 90% CL. values in red (blue) are new published (preliminary) result since 
PDG2004 [as of August 25th, 2004]. 



RPP# 


Mode 


PDG2004 Avg. 


BaBar 


Belle 


CLEO 


New Avg. 


162 


if* (892)07 


43 ±4 


39.2 ±2.0 ±2.4 


40.1 ± 2.1 ± 1.7 


45.5tg:g ±3.4 


40.1 ± 2.0 


163 




< 8.3 




< 8.3 


< 8.3 


164 


K^n~^ t 


4.6 ± 1.4 




. f;+1.3 + 0.5 

^■°-1.2-0.7 




4.6 ± 1.4 






New 




24.3 ±3.6 ±3.4 




24.3 ± 5.0 


165 


X*(1410)07 


< 130 




< 130 




< 130 


166 


K+TT-f (N.R.) t 


< 2.6 




< 2.6 




< 2.6 


169 


X| (1430)07 


13 ± 5 


12.2 ±2.5 ±1.0 


13±5±1 




12.4 ±2.4 


173 


P°7 


< 1.2 


< 0.4 


< 0.8 


< 17 


< 0.4 


174 




< 1.0 


< 1.0 


< 0.8 


< 9.2 


< 0.8 


175 


07 


< 3.3 


< 0.94 




< 3.3 


< 0.94 


237 


ROe+e- 


< 0.54 


-0.2lj:°-^g ±0.08 t 


n |-,n+0-20+0.02 ^ 
U.UU_Q j^2-0.05 + 


< 8.45 


-O.O6I0 10 


238 






1.63lo;g±0.14 


0.56lo 2I ± 0.05 


< 6.64 




240 


K*{892fe+e- 


< 2.4 


i.iiig-j? ±0.11 


1 r,n+6.57+0.13 
J^-^M_0.49-0.10 




1 20+°-*^ 


241 


X*(892)0^i+^- 


1.3 ±0.4 


o.86i;;;^l ±0.11 


1.33 j; 55:;^^ ± 0.11 




1 99+0.39 
J^-^^-0.33 



t 1.25 GeV/c2 < A/xtt < 1-6 GeV/c^ t Central values are not significant. 



Table 45: Compilation of B semileptonic and radiative branching fractions (in units of 10 ^). 
Upper limits are at 90% CL. values in red (blue) are new published (preliminary) result since 
PDG2004 [as of August 25th, 2004]. 



RPP# 


Mode 


PDG2004 Avg. 


BaBar 


Belle 


CLEO 


New Avg. 


60 


^3* (1780)7 


< 3000 




< 34 




< 34 


67 


s-y 


330 ± 40 


388 ± 36^45 


355 ±32t|Jlii 


321 ± 43^29 


ot;9 + 30 




S7 with baryons 


New 






< 38 t 


< 38 t 


71 


fn 


< 1.9 


< 1.2 


< 1.4 


< 14 


< 1.2 






New 




f, q+l.7+1.3 
0-^-1.6-1.0 






101 


se+e~ X 


5.0 ±2.6 


6.0 ± 1.7 ± 1.3 


4.04± 1.30to 7g 


< 57 


4.70lli^ 


102 




7 0+30 
' -^-2.6 


5.0 ±2.8 ±1.2 


4.13 ± 1.051:0^1 


< 58 


4 26+1-18 
^-^"-1.16 


103 


se+£- t 




5.6 ± 1.5 ± 1.3 


4.11 ±0.83tg-^^ 


< 42 


4 46+0-98 
^•™-0.96 


104 


Ke+e- 


0.4810:11 


0.74lg;J^ ± 0.05 


n /I t:4+0. 116+0.023 

U.^O'i_Q jg4_Q Q25 




n ^547+0. 098 


105 


X*(892)e+e- 


1.5 ±0.5 


0.98lo;42 ± 0.11 


1.8410-48^0.17 




1 44+0.35 


106 


K/i+tJ.- 


0.48 ±0.12 


OAbto ll ± 0.04 


n cnc+0. 103+0.033 
U.DZD_o Qg4_o 




n 605+0-090 

U.DUO_(j Qg4 


107 


K*{892)ii+/i- 


1 17+0.37 

^•^'-0.33 


1.27lg;^^ ±0.16 


1.8110^0 ±0.11 
0.550to o70 ± 0.027 




1 70+0.30 
''-'-0.27 


108 


Kl+f.- 


5.4 ± 0.8 


0.65to \t ± 0-04 


< 1.7 


n 574+0.071 

U.0^4_o Qgg 


109 


K*(892)e+e- 


1.05 ± 0.20 


O-SSiaS ± 0-10 


1.65l°J3±0.10 


< 3.3 


1.38 ±0.20 


111 




< 1.6 






< 1.6 


< 1.6 


112 




< 3.2 






< 3.2 


< 3.2 


113 




< 1.6 






< 1.6 


< 1.6 


114 




< 6.2 






< 6.2 


< 6.2 



t > 2.0 GeV; {£+£-) > 0.2 GeV/c^ 
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Table 46: Compilation of B leptonic branching fractions (in units of 10 ^). Upper limits are 
at 90% CL. values in red (blue) are new published (preliminary) result since PDG2004 [as of 
August 25th, 2004]. 





Mode 


PDG2004 Ave 


RaRar 


Relle 


CLEO 


CDF 


DO New Ave: 


12 




< 15 




< 5.4 


< 15 




< 5.4 


13 




< 21 


< 6.6 


< 2.0 


< 21 




< 2.0 


14 


r+z/ 


< 570 


< 330 


< 290 


< 840 




< 290 


15 




< 200 




< 22 


< 200 




< 22 


16 




< 52 




< 23 


< 52 




< 23 


235 


e~^e~ 


< 0.19 


< 0.06 


< 0.19 


< 0.83 




< 0.06 


236 




< 0.16 


< 0.08 


< 0.16 


< 0.61 


< 0.15 


< 0.08 


244 




< 0.17 


< 0.18 


< 0.17 


< 1.5 




< 0.17 


247 




< 530 






< 110 




< 110 


248 




< 830 






< 38 




< 38 




vV 


New 


< 220 








< 220 






New 


< 47 








< 47 



Table 47: Compilation of Bs leptonic branching fractions (in units of 10 ^). Upper limits are 
at 90% CL. values in red (blue) are new published (preliminary) result since PDG2004 [as of 
August 25th, 2004]. 



RPP# Mode 


PDG2004 Avg. 


CDF 


DO 


New Avg. 


— /.i^/i" 


New 


< 0.58 


< 0.41 


< 0.41 
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6.3 Baryonic decays 



Table 48: Compilation of B'^ baryonic branching fractions (in units of 10 Upper limits are 
at 90% CL. values in red (blue) are new published (preliminary) result since PDG2004 [as of 
August 25th, 2004]. 



RPP# 


Mode 


PDG2004 Avg. 


BABAR 


Belle 


CLEO 


New Avg. 


201 
204 

206 


ppK*+ 

pA 
AAK+ 
AA-K+ 


< 3.7 

New 

< 1.5 
New 
New 


6.7 ±0.9 ±0.6 


3.06to g2 ± 0.37 
5.66+0 57 ±0.62 

in Qi +3.62+1.34 
lu.o±_2 g5 

< 0.46 

2.9llo.TO ± 0-38 

< 2.8 


< 160 

< 1.5 


3.06+n S 

— 0.72 

6.08lgS 

10-3ll3.22 

< 0.46 

2 Q1+0.98 

< 2.8 


214 


Acp-n-+ 


210 ± 70 




201 ± 15 ± 56 


240 ± eolga 


213 ± 48 


215 
216 


Acpn+ir'^ 

A'^p7T~^7T'^7T~ 


1800 ± 600 
2300 ± 700 






1810 ± 290lgoo 
2250 ± 250lgio 


2250l«g 


218 


Z^(2455)p 


< 80 




36.7lg;g ± 10.2 


< 80 




219 


Z^(2520)p 
X0(3340)7r+ 


< 46 
New 




I2.6I4 g ± 3-5 
38.7tl'l ± 11.0 




12.6l«:g 


220 


r°(2455)p7r0 


440 ± 180 






420 ± 130 ± 170 


420 ± 214 


221 


X^^(2455)p7r+7r- 


440 ±170 






440 ± 120 ± 120 


440 ± 169 


222 


I:;:"(2455)p7r+7r+ 


280 ± 120 






280 ± 90 ± 90 


280 ± 127 


223 


^;r(2593)p7r+ 


< 190 






< 190 


< 190 



Table 49: Compilation of baryonic branching fractions (in units of 10 ^). Upper limits are 
at 90% CL. values in red (blue) are new published (preliminary) result since PDG2004 [as of 
August 25th, 2004]. 



RPP# 


Mode 


PDG2004 Avg. 


BABAR 


Belle 


CLEO 


New Avg. 


212 


pp 


< 1.2 


< 0.27 


< 0.50 


< 1.4 


< 0.27 


214 


ppK° 


< 7.2 




1.88lJ;eQ ± 0.23 




l-»8_0.64 




ppK*° 


New 




< 7.6 




< 7.6 


215 


pA-K~ 






3.97tJ;go ± 0.56 


< 13 


Q q7+l'l5 
•^■^'-0.98 


216 


pAK- 


< 0.82 




< 0.82 




< 0.82 


217 


pZj TT 


< 3.8 




< 3.8 




< 3.8 


218 


AA 


< 1.0 




< 0.79 


< 1.2 


< 0.79 


224 


A^ pVT + TT" 


1300 ± 400 




1030 ± 90 ± 295 


1670 ± 1901:1™ 


1207 ± 262 


225 




22 ±8 




21.9t4 g ± 6.5 


< 90 


21.9tli 


229 


E~''(2520)p7r+ 


160 ± 70 




104 ± 23 ± 30 




104 ± 37 


230 


r"(2520)p7r+ 


< 121 




33 ± 19 ± 10 




33 ±21 


231 


X^(2455)p7r- 


100 ± 80 




97 ± 21 ± 30 


220 ± 60 ± 64 


115 ±33 


232 


I:;:"(2455)p7r+ 


280 ± 90 




115 ±22 ±33 


370 ± 80 ± 113 


134 ± 38 


233 


yi;:(2593)p 


< 110 






< 110 


< 110 



6.4 Charge asymmetries 
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Table 50: Compilation of charmless hadronic CP asymmetries for charged B decays. Values in red (blue) are new 
published (preliminary) result since PDG2004 [as of August 25th, 2004]. 



RPP# 


Mode 


PDG2004 Avg. 


BABAR, 


Belle 


CLEO 


CDF 


New Avg. 


117 




0.03 ±0.08 


-0.087 ±0.046 ±0.010 


0.05 ±0.05 ±0.01 


0.18 ±0.24 ±0.02 




-0.020 ±0.034 


lio 


J\ TV 


— U.iU it U.Uo 


U.UD it U.UD it U.Ui 


U.U4 it U.UO it U.U.^ 


n oo -U n -U n no 
— U.zy it U.Zo it U.Uz 




U.U4 it U.U4 


119 


ri'K+ 


0.009 ±0.035 


0.037 ±0.045 ±0.011 


-0.015 ±0.070 ±0.009 


0.03 ±0.12 ±0.02 




0.022 ± 0.037 


121 


r)K+ 


New 


-0.52 ±0.24 ±0.01 


-0.49 ±0.31 ±0.07 






-0.51 ±0.19 


122 


riK*+ 


New 


0.13 ±0.14 ±0.02 


-0.09to ;tg ±0.01 








123 


LUK + 


-0.21 ± 0.28 ± 0.03 


-0.09 ±0.17 ±0.01 


0.06io:ig ±0.01 






-0.02l»;i3 


127 


K + TT+TT- 


0.01 ±0.07 ±0.03 


0.01 ±0.07 ±0.03 






0.01 ±0.08 




K'°p+ 


New 


-0.14 ±0.17 ±0.4 








-0.14 ±0.43 


138 




0.20j:j;;^g ± 0.04 


0.20j;o2g ± 0.04 










144 


K+KsKs 


New 


-0.04 ±0.11 ±0.02 








-0.04 ±0.11 


152 


K+K-K+ 


0.02 ±0.07 ±0.03 


0.02 ±0.07 ±0.03 








0.02 ± 0.08 


153 


<PK+ 


0.03 ±0.07 


0.054 ±0.056 ±0.012 


0.01 ±0.12 ±0.05 






0.038 ± 0.050 


156 


<i>K*+ 


0.09 ±0.15 


0.16 ±0.17 ±0.03 


-0.13±0.29lEJ°f 






0.09 ±0.15 


173 




0.05 ±0.15 


-0.01 ±0.10 ±0.02 


-0.02 ±0.10 ±0.01 






-0.02 ±0.07 


174 




-0.39 ± 0.33 ± 0.12 


-0.39 ±0.33 ±0.12 








-0.39 ±0.35 


175 




New 


-0.19 ±0.11 ±0.02 








-0.19 ±0.11 


180 


P+ttO 


New 


0.24 ±0.16 ±0.06 


0.06±0.19tE5;o6 






0.16 ±0.13 


182 


p+p° 


-0.09 ± 0.16 


-0.19 ±0.23 ±0.03 


0.00 ± 0.22 ± 0.03 






-0.09 ±0.16 


185 


a;7r+ 


-0.21 ± 0.19 


0.03 ±0.16 ±0.01 


0.50to;2o ±0.02 


-0.34 ±0.25 ±0.02 




0.10 ±0.11 


187 




New 


-0.44 ±0.18 ±0.01 


0.07 ±0.15 ±0.03 






-0.14 ±0.12 


188 




New 


0.24 ±0.19 ±0.01 








0.24 ±0.19 


190 


r)p+ 


New 


0.07 ±0.19 ±0.02 


-0.17i°;2g ±0.02 






01+° " 



Table 51: Compilation of charmless hadronic CP asymmetries for B"^ / admixtures. Values in red 
(blue) are new published (preliminary) result since PDG2004 [as of August 25th, 2004]. 



RPP# 


Mode 


PDG2004 Avg. 


BABAR 


Belle 


CLEO CDF 


New Avg. 


56 




-0.01 ± 0.07 


-0.013 ±0.036 ±0.010 


-0.015 ±0.044 ±0.012 


0.08 ±0.13 ±0.03 


-0.010 ±0.028 


67 




-0.079 ±0.108 ±0.022 


0.025 ±0.050 ±0.015 


0.002 ±0.050 ±0.026 


-0.079 ±0.108 ±0.022 


0.005 ± 0.036 


103 


see 


New 


-0.22 ±0.26 ±0.02 






-0.22 ±0.26 



Table 52: Compilation of charmless hadronic CP asymmetries for neutral B decays. Values in red (blue) are new 
published (preliminary) result since PDG2004 [as of August 25th, 2004]. 



RPP# 


Mode 


PDG2004 Avg. 


BABAR 


Belle 


CLEO 


CDF 


New Avg. 


123 


K + 7T- 


-0.09 ±0.04 


-0.133 ±0.030 ±0.009 


-0.101 ±0.025 ±0.005 


-0.04 ±0.16 ±0.02 


-0.04 ±0.08 ±0.006 


-0.109 ± 0.019 


124 


ROnO t 














127 




New 


0.02 ±0.11 ±0.02 


-OMtoil, ±0.01 






-0.01 ± 0.08 


136 


K+p- 


0.28 ± 0.19 


0.13lg ±0.14 


n Qr,+0.'22+0.06 
'-'•^^-0.23-0.02 






17+0'15 
'^•^'-0.16 




K*(U30) + n- 


New 


-0.07 ±0.12 ±0.08 








-0.07 ±0.14 






New 


-0.34 ±0.15 ±0.11 








-0.34 ±0.19 


140 


K*+n- 


0.26 ±0.35 


-0.04 ±0.13 




n 9f;+0.33+0.10 
'^■^"-0.34-0.08 




-0.00 ±0.12 




K+7T-1T° 


New 




0.07 ±0.11 ±0.01 






0.07 ±0.11 


141 




New 


-0-Ollo.22±0-13 


0.07±0.15j;o;°3 






-0 01+"-^'^ 


154 




0.05 ± 0.10 


0.04 ±0.12 ±0.02 






0.05 ±0.10 


177 




New 


0.12 ±0.56 ±0.06 


0.43 ± O.blto.ll 






0.28 ±0.39 


197 


p+TT" t 













00 



f Measurements of time-dependent CP asymmetries are listed on the Unitarity Triangle home page, 
(http:/ /www. slac.stanford.edu/xorg/hfag/triangle/index. html) 



Table 53: Brief summary of the world averages as of 2004 summer conferences. 



b-hadron lifetimes 

r{Bt) 


1.534 ±0.013 ps 
1.653 ±0.014 ps 
1.469 ± 0.059 ps 
0.45 ±0.12 ps 
1.232 ±0.072 ps 


b-hadron fractions 

in T{AS) decays 
fd = fu at high energy 

at high energy 
/baryon at high energy 


1.026 ±0.032 
0.398 ±0.010 
0.104 ±0.015 
0.100 ±0.017 


and mixing parameters 

Iq/pU 

Arris 


0.502 ±0.006 ps^i 
1.0013 ±0.0034 
> 14.5 ps-i at 95% CL 


Semileptonic B decay parameters 

^ D+i~V) 
b\b ^ XiV) 

\Vcb\ (B° D*+e-v) 

\Va\ (b" D+i-V) 
\Vub\ (inclusive) 


(K QQ _)- n 201% 

(2.13 ±0.20)% 
(10.90 ±0.23)% 
[41.4 ± l.O(exp) ± 1.8(theo)] x 10"^ 
[40.4 ± 3.6(exp) ± 2.3(theo)] x lO'^ 
(4.70 ±0.44) X 10"3 


CP{t) and Unitarity Triangle angles 

sin2/3(</)i) (all charmonium) 
sin2/3(0i)cfT (all 6 — > s penguin) 


0.726 ±0.037 
0.41 ±0.07 


Rare B decays 

Acp{B^ ^ A'+TT-) 


-0.109 ±0.019 (5.7(t) 


Summary 



This article provides the updated world averages for 6-hadron properties as of 2004 summer 
conferences (ICHEP04 and FPCP04). A brief summary of the results described in Sees. OHHlis 
given in Table ESI 

The HFAG provided the first averages at 2003 winter conferences and have given updates 
at summer and winter conferences, as well as annual PDG averages. 

The accuracies of 6-hadron lifetimes and B^ mixing parameters have been considerably 
improved with the asymmetric B factory results compared to the previous LEP working group 
averages j2]. In 2004 summer, the BABAR and Belle collaborations reported the simultaneous 
measurements of B lifetime and Arrid with improved precision using increased data samples. 
Because of correlation between lifetime and Am^ which complicates the average procedures, 
these new results are not included in current averages and there are no substantial changes from 
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previous averages at 2004 winter. It is planned to incorporate them in the next update. In the 
previous update, new results of 6-hadron fractions from CDF are included and they remain the 
same in this update. The average of f~^~/f^^ was newly added to HFAG average items from 
the previous update. With new results from BABAR, the accuracy has been improved to ~ 3% 
and the average value is consistent with 1. In this update, the average of \q/p\d is newly added. 

For |V"cfe| and \Vub\ average values, various new measurements have been available from 
CLEO, BABAR, and Belle. Accordingly, the statistical uncertainties of averages are significantly 
reduced from LEP working group averages j2]. Considerable progress has been also made in 
theoretical side, but the reduction of uncertainty is somewhat slower in some cases [e.g. \Vcb\ in 
exclusive modes). The determination of \Vcb\ from the inclusive semileptonic branching fraction 
in combination with hadronic mass and lepton energy moments is under discussion and foreseen 
in the future update. The average of \ Vub\ using inclusive b — > uiV is provided from the previous 
update. In this update, a substantial change has occurred with the usage of the Belle photon 
spectrum measurement in 6 — S7 decays. This is used for the determination of shape function 
parameters which enter the \Vub\ calculation in inclusive charmless semileptonic B decays. The 
new shape function parameters result in substantially smaller error in \Vub\ than previous one. 

Measurements by BABAR and Belle of the time-dependent CP violation parameter Shades 
in B decays to charmonium and a neutral kaon have established CP violation in B decays, 
and allow a precise extraction of the Unitarity Triangle parameter sin2/3/ sin20i. The B facto- 
ries have also provided various measurements of time-dependent CP asymmetries in hadronic 
6 — * s penguin decays, establishing CP violation in these modes. Intriguingly, the measured 
parameters exhibit deviations from the Standard Model expectation. The significance of this 
effect depends on the treatment of the theoretical error. Results from time-dependent analyses 
with the decays B° — > 7r"'"7r~, p^vr^ and p'^p~ allow, via various methods, constraints on the 
Unitarity Triangle angle a/(f)2- Constraints on the third Unitarity Triangle angle 7/(^3 have 
been obtained by BABAR and Belle, using B^ D^*^K~ decays with Dalitz plot analysis of 
the subsequent D KgH'^7i~ decay. 

For the rare B decays, branching fractions and charge asymmetries of many new decay modes 
have been measured recently, mostly by BABAR and Belle. Since there are several hundred 
measurements in the tables in SeclHl we highlight only the measurement of Acp{B^ — > K~^7r~), 
which provides the observation of direct CP violation (5.7 a) in in the B meson system. 
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